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This  thesis  presents  the  theory  necessary  for  the  concep¬ 
tual  design  of  a  tactical  missile.  The  design  process  begins 
with  the  well  known  linear  aerodynamic  theory  for  initial  sizing 
and  later  includes  nonlinear  effects  to  determine  the  final 
design  of  the  missile.  Where  theory  does  not  apply,  empirical 
methods  are  presented  which  are  known  to  give  accurate  results. 
An  air-to-air  missile  is  designed  for  a  specific  threat  as 
an  example  which  immediately  follows  the  development  of  the 
theory  for  each  section.  Several  small  digital  computer  pro¬ 
grams  are  presented  and  used  for  analysis  of  specific  areas 
of  the  design.  One  large  program  (AER01)  is  used  for  deter¬ 
mining  the  aerodynamic  coefficients  of  the  final  design. 
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I.  INTRODUCTION 


This  study  was  made  to  present  a  method  for  the  conceptual 
design  of  tactical  missiles.  The  starting  point  for  the  design 
was  a  recently  completed  report  by  General  Dynamics,  Convair 
Division  entitled  Rapid  Approach  to  Missile  Synthesis  (RAMS) .  A 
procedure  was  then  developed  more  akin  to  aircraft  design  and 
which  bears  little  resemblance  to  RAMS  and  which  uses  basic 
equations  to  size  components  rather  than  nomograms  and  table 
look-up.  The  procedure  starts  with  a  threat  description  and 
proceeds  with  the  formulation  of  performance  objectives  and 
a  conceptual  design  of  a  tactical  missile  to  counter  the  threat. 
The  design  is  not  unique,  but  as  will  be  shown,  is  a  compromise 
of  parameters  to  give  one  possible  solution  to  the  design  prob¬ 
lem;  therefore,  the  point  design  arrived  at  is  not  necessarily 
the  optimum  design  for  the  presented  threat.  An  attempt  is 
made  to  find  the  optimum  performance  within  specific  areas  of 
the  design  process. 

Throughout  this  study  the  theory  involved  is  explained  and 
specific  examples  are  worked.  A  complete  design  example  is 
worked  out  in  detail.  It  is  an  air-to-air  missile  designed  to 
counter  the  new  Soviet  RAM-K  fighter  aircraft.  This  example 
is  worked  in  each  section  immediately  following  the  develop¬ 
ment  of  the  theory  for  that  specific  area. 
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II.  PROBLEM  DEFINITION 


A.  THREAT  ANALYSIS 

1.  Operational  Requirements 

The  design  of  a  new  missile  is  usually  in  response  to 
an  operational  requirement  which  arises  as  the  result  of  one 
of  the  following:  (1)  A  new  technology  provides  the  means  to 
design  a  more  effective  missile  to  meet  a  current  threat.  An 
example  of  this  might  be  an  advance  in  material  science,  which 
allows  higher  inlet  turbine  temperatures  for  a  turbojet  engine; 
therefore,  allowing  higher  missile  flight  speeds.  (2)  Intelli¬ 
gence  indicates  a  new  threat  for  which  existing  missiles  are 
not  effective.  (3)  Operational  reports  indicate  a  current  mis¬ 
sile  is  inadequate  against  a  current  threat. 

Regardless  of  how  the  operational  requirement  is  derived, 
a  statement  of  the  threat  is  required  before  the  design  process 
can  proceed.  Experience  has  proven  that  one  missile  cannot  be 
designed  to  meet  all  types  of  threats  without  seriously  com¬ 
promising  performance  or  effectiveness.  This  can  be  illustrated 
with  the  design  of  the  warhead.  A  contact  fuze,  shaped  charge 
warhead  designed  to  penetrate  and  kill  hard  targets  such  as 
tanks,  would  not  be  effective  against  a  highly  maneuverable 
aircraft  for  which  the  expected  miss  distance  is  several  feet. 

For  this  reason  the  design  of  a  missile  must  start  with  a  de¬ 
tailed  analysis  of  the  threat.  The  more  detailed  this  analysis 
is,  the  more  effective  the  final  design  can  be. 


2 .  Design  Example  (Operational  Requirement) 

A  design  example  will  be  used  as  a  continuous  thread 
throughout  this  thesis  to  demonstrate  applications  of  the  theory. 
An  air-to-air  missile  will  be  designed  to  counter  the  new  class 
of  Soviet  fighter,  which  is  in  the  advanced  development  stage 
at  the  Ramenskoye  Experimentation  Center.  The  fighter,  as  des¬ 
cribed  in  Aviation  Week(l,23has  been  designated  the  RAM-K.  The 
RAM-K  is  a  twin  engine  fighter  with  variable  geometry  inlets  and 
swing  wings.  The  aircraft  bears  a  resemblance  to  both  the  F-14 
and  F-15.  It  is  expected  to  be  the  recipient  of  a  new  look 
down,  shoot  down  radar  and  the  40  km  range  AA-X-9  missile.  The 
following  unclassified  dimensions  and  performance  data  are  avail¬ 
able  on  the  RAM-K: 

Wing  span  40  ft 

Overall  length  64  ft 

Gross  weight  6Q000  lbs 

Maximum  speed  M=2.5 

Service  ceiling  60,000  ft 

Figure  (2-1)  is  a  drawing  of  the  RAM-K. 

3 .  Scenario 

The  scenario  within  which  the  missile  is  expected  to 
operate  should  also  be  described.  If  the  normal  mode  of  operation 
of  the  threat  is  not  known,  an  attempt  should  be  made  to  define 
the  most  demanding  scenario  that  can  be  expected.  For  a  defen¬ 
sive  weapon  the  most  challenging  incoming  threat  will  normally 
be  a  head-on  encounter.  The  threat  profile  may  vary  from  a  high 
level  attack  with  a  terminal  dive  to  a  low  level  attack  with  a 


terminal  pop-up  maneuver.  In  the  case  of  the  AS-6  (Kingfish) 
two  modes  of  attack  can  be  expected.  In  such  a  case  both  pro¬ 
files  must  be  evaluated  to  determine  the  most  demanding  in 
terms  of  missile  performance  objectives.  For  an  offensive  sys¬ 
tem,  such  as  an  air-to-air  missile  designed  to  intercept  and 
destroy  an  enemy  fighter  before  it  launches  its  weapons,  the 
scenarios  analyzed  should  include  all  possible  encounter 
geometries . 

4 .  Design  Example  (Scenario) 

The  scenario  for  the  above  threat  would  likely  be  an 
intercept  situation  in  defense  of  the  fleet  high  value  unit. 

The  scenario  is  taken  to  be  a  head-on  encounter  with  the  missile 
and  the  target  at  the  same  altitude.  Since  the  combat  specifi¬ 
cations  of  an  aircraft  are  normally  given  at  10,000  ft.,  this  is 
taken  as  the  scenario  altitude. 

B.  HISTORICAL  SURVEY 

Missile  design  is  an  iterative  process,  and  the  first  time 
through  the  design  loop  many  assumptions  have  to  be  made  con¬ 
cerning  component  sizes  and  weights.  One  method  of  approach 
at  the  early  stage  is  to  employ  historical  data  of  existing  mis¬ 
sile  sizes  and  weights;  since  justifications  for  these  parameters 
were  made  duirng  their  design  processes.  An  example  of  the  use 
of  historical  data  in  determining  the  initial  missile  length 
can  be  made  with  the  length  to  diameter  ratio.  The  length  to 
diameter  ratios  of  existing  missiles  of  the  same  class  as  that 
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being  designed  are  collected,  and  an  average  is  computed.  The 
diameter  of  the  design  is  fixed  by  one  of  three  driving  factors 
(propulsion,  warhead,  or  guidance) .  From  the  average  length  to 
diameter  ratio  the  initial  missile  length  is  then  estimated. 

From  this  historical  data,  initial  choices  based  on  the  experi¬ 
ence  of  others  can  be  made  for  many  of  the  missile  parameters. 
These  parameters  define  a  baseline  missile,  which  is  the  initial 
configuration  from  which  design  iterations  and  refinements  can 
be  made. 

Since  missiles  are  designed  for  specific  missions,  speci¬ 
fic  parameters  such  as  length  and  diameter  are  of  little  value 
in  comparing  missiles.  Dimensionless  ratios  such  as  length  to 
diameter,  L/D,  ratio  and  aspect  ratio,  AR,  are  more  meaningful 
when  relating  missiles.  Some  parameters  which  are  useful  in 
defining  the  baseline  missile  are  listed  below: 

L/D  =  Length  to  diameter  ratio 

L  /D  =  Nose  length  to  diameter  ratio 
n 

ARw  =  Aspect  ratio  of  the  wing 

AR^  =  Aspect  ratio  of  the  tail 

W/S  =  Weight  to  lifting  surface  area  ratio 

V,  =  S^l^/fS  ,d  _)  Tail  volume  coefficient 
t  t  t  ref  ref 

WG/Ww^  =  Gross  weight  to  warhead  weight  ratio 

The  tail  volume  coefficient,  Vt,  is  a  dimensionless  para¬ 
meter  used  to  initially  size  the  tail.  For  a  tail  control  mis¬ 
sile,  it  is  a  measure  of  the  relative  control  effectiveness  when 
comparing  missiles.  For  a  wing  control  missile  it  is  a  relative 
measure  of  stability. 
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A  complete  historical  survey  should  not  be  limited  to 
the  parameters  listed  here.  Any  dimensionless  parameter  which 
will  add  information  about  the  proposed  design  should  be  inclu¬ 
ded  for  completeness. 

Table  2-1  is  an  example  of  a  collection  of  such  para¬ 
meters  for  existing  air-to-air  missiles  El,  3] .  in  this  table 
the  subscript,  c,  is  used  to  indicate  a  canard  control  surface. 

C.  LAUNCH  PLATFORMS  AND  PHYSICAL  CONSTRAINTS 

The  problem  definition  must  also  include  a  description  of 
the  intended  launch  platform  for  the  missile.  The  aircraft  or 
shipboard  launcher  from  which  the  missile  will  be  launched  will 
fix  many  design  features  of  the  missile.  For  instance  the  most 
important  consideration  in  the  problem  definition  phase  is  any 
physical  constraints  imposed  by  the  launcher.  Since  it  is  not 
normally  feasible,  economically  to  design  a  launcher  to  fit  the 
missile,  most  new  missiles  must  fit  existing  launchers.  In  the 
case  of  shipboard  launchers  there  will  be  a  maximum  length  and 
diameter  and  a  maximum  launch  weight  which  can  be  accomodated. 

For  the  case  of  an  air  launched  missile,  there  will  be  a 
maximum  weight,  and  the  dimensions  may  also  be  limited  due  to 
the  performance  requirements  of  the  aircraft. 

For  the  air-to-air  missile  design  example  of  this  study, 
the  launch  platforms  will  be  the  F-16  and  F-18A.  Figure  (2-2) , 
which  is  from  Interavia  [4, 5]  shows  the  pylon  weight  limitations 
for  these  aircraft.  From  these  figures  it  can  be  seen  that  the 
maximum  launch  weight  for  this  design  is  limited  to  2500  pounds 

by  the  pylon  limitations  of  the  F-18A. 
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Load  Carrying  Capability  of  F-18A 


Load  Carrying  Capability  of  F-16 


2200/5.5  2500/5.5 


Figure  (2-2).  Weight  Limitations. 


D.  MISSION  PROFILES  AND  PERFORMANCE  OBJECTIVES 

The  mission  profile  of  a  missile  consists  of  dividing  the 
flight  into  fundamental  segments  which  consist  of  a  single  func¬ 
tion,  such  as  boost  to  cruise  speed  and  altitude,  cruise  to 
the  target,  and  terminal  homing  phases.  The  mission  profile 
will  vary  from  missile  to  missile.  For  a  cruise  missile  it 
may  consist  of  a  series  of  pop-up  maneuvers  and  low  level  cruises. 
For  a  short  range  missile,  the  profile  may  be  entirely  terminal 
homing.  The  Mach  number  and  altitude  are  specified  at  the  begin¬ 
ning  and  end  of  each  mission  segment,  as  well  as  the  range 
covered  by  each  segment.  The  range  covered  by  a  segment  can  be 
considered  in  one  of  two  ways.  If  the  segment  is  short,  such 
as  the  terminal  phase,  the  distance  along  the  intended  flight 
path  is  considered.  For  longer  range  missiles,  the  distance  over 
the  ground  is  of  importance.  The  mission  profile  must  be  defined 
during  the  problem  definition  phase  in  order  to  specify  missile 
performance  objectives. 

The  mission  profile  of  a  missile  normally  consists  of  a 
boost,  cruise  (mid-course)  and  terminal  phase.  The  boost  phase 
accelerates  the  missile  to  its  flight  speed.  This  acceleration 
may  be  large  for  a  surface-to-air  missile  which  must  be  acceler¬ 
ated  from  rest  to  a  high  supersonic  speed  or  it  may  be  small  for 
an  air  launched  missile  which  has  the  speed  advantage  of  the  air¬ 
craft  from  which  it  is  launched. 

The  cruise  segment,  or  mid-course  phase,  primarily  is  used 
to  deliver  the  missile  to  a  point  in  space  where  the  seeker  can 
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acquire  the  target.  The  range  and  speed  of  the  cruise  segment 
is  then  a  function  of  the  required  stand-off  distance  for  the 
target. 

The  terminal  phase  is  somewhat  more  difficult  to  analyze. 

If  the  target  maintains  constant  heading  and  velocity  the  flight 
path  of  the  missile  may  be  modeled  by  circular  arc  segments.  This 
method  will  give  approximate  values  of  range.  If  the  target 
maneuvers,  the  missile  must  follow  and  the  range  and  speed  require¬ 
ments  become  more  complicated. 

In  determining  the  range  and  speed  requirements,  all  expected 
encounter  geometries  should  be  analyzed.  The  most  demanding 
encounter  will  then  fix  the  performance  objectives.  The  most 
demanding  speed  requirement,  in  terms  of  maintaining  a  minimum 
stand-off  distance,  will  normally  be  a  head-on  encounter.  Al¬ 
though  the  required  missile  speed  and  range  are  determined  in 
this  section,  the  missile  velocity  may  be  varied  later  in  the 
design  process  due  to  guidance  considerations. 

1.  Design  Example  (Mission  Profile? 

From  the  threat  defined  above,  the  ideal  situation  would 
be  to  obtain  a  fire  control  solution  and  launch  such  that  the 
minimum  separation  distance  between  the  launch  aircraft  and  the 
target  is  40  km.  This  can  be  accomplished  in  one  of  two  ways. 

The  launch  aircraft  can  fire  a  semi-active  homing  missile  at 
such  a  range  and  speed  that  intercept  occurs  before  the  minimum 
range  is  reached,  or  an  active  homing  missile  can  be  fired,  and 
once  missile  lock-on  is  achieved,  the  launch  aircraft  can  maneu¬ 
ver  to  maintain  the  minimum  separation  distance. 
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The  active  radar  homing  missile  would  decrease  the 
launch  range  and  the  range  required  of  the  missile,  which  will 
lessen  the  cons taints  on  the  launch  aircraft.  The  terminal  por¬ 
tion  of  the  engagement  is  a  functicnof  the  guidance  law  and  will 
be  determined  in  Chapter  3. 

Two  cases  are  investigated  to  determine  the  effect  on 
the  range  requirement  of  the  missile  when  a  minimum  separation 
distance  from  the  target  to  the  missile  of  40  1cm  is  maintained. 
The  first  case  is  a  semi-active  homing  missile,  for  which  the 
launch  aircraft  must  maintain  a  closing  course  until  intercept. 
The  second  case  is  an  active  homing  missile  which  has  a  lock-on 
range  of  10  km.  The  launch  aircraft  may  then  maneuver  to  main¬ 
tain  a  separation  distance. 

a.  Case  1:  Semi-active  homing  missile 


< -  Rmt 


*LT 


Rq  =  Range 

«L  "  1'!  - 

“m  *  2’5  " 


a  =  Speed 


V.  = 


V„  = 


M  a  = 

L 

Mja  = 


VM  - 


M  a  = 
M 


at  which  the  missile  is  launched 
Launch  Mach  number 
Missile  Cruise  Mach  number 
of  sound 

1.5a  =  Launch  speed 
2.5a  =  Target  speed 
2.5a  =  Missile  speed 
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The  instantaneous  range  from  the  missile  to  the  target  is 
given  by, 


*MT  “  R0  -  (VH  +  V* 


(1) 


The  instantaneous  range  from  launch  aircraft  to  the  target 
is  given  by,  R^ 


*LT  *  R0  -  (V1  +  V4 


(2) 


If  the  target  does  not  maneuver,  intercept  will  occur  at 
tf,  when  Rj^  =*  0 


0  = 


t„  = 


R0  '  (VM  +  VT)fcf 


R 


_ q___  a 

(v“+vj  5a 


M  T' 


If  the  launch  aircraft  is  at  the  minimum  separation  distance, 
R^t  =  40  km,  when  intercept  occurs. 


40  km  =  R0  -  (VL  +  VT)tf 


substituting  for  VL»VT  and  tf 


40  km  =  Rq  -  (1.5  +  2.5)a(-jr^) 


Solving  for  R 


rq  =  200  km  =  Launch  range 


The  range  required  of  the  missile  is  then,  R^. 


=  V^tf  =  100  km  =  53.96  nmiles 

If  the  missile  speed  is  increased  to  MM  =  3.0,  then 
Rq  =  148 . 15  k  and, 

Rj^  =  80.81  km  =  43.61  nmiles. 

b.  Case  2:  Active  Homing  Missile 

The  lock-on  range  is  a  function  of  the  seeker  in  the 
missile  and  will  be  covered  later  in  this  thesis.  If  it  is 
assumed  that  the  launch  aircraft  must  maintain  its  course  until 
lock-on  occurs  at  a  range  of  R^,  the  problem  can  still  be  solved. 
The  geometry  is  the  same  as  in  Case  1.  Instead  of  following  a 
constant  course  until  intercept,  the  launch  aircraft  must  now 
only  maintain  a  closing  course  until  RMT  =  R^.  Then  from 
equation  (1) 

*MT  -  \0  ~  R0  '  (VM  +  V'fl 

Solving  for  tfl 


ro'rlo 

£1  w 

If  at  the  time  of  target  lock-on,  tfl,  the  target 
and  launch  aircraft  are  at  the  minimum  separation  distance, 
Rlt  =  Rmin*  from  equation  (2) , 


(3) 


For  the  same  geometry  and  relative  speeds  of  the 

first  case,  with  =  3.0, 

M 

t  =  V”lo 

fl  5.5a 

A  reasonable  value  of  lock-on  range  is  10  km.  This 
will  be  shown  later  in  the  guidance  section  of  the  study.  The 
time  to  lock-on  then  becomes , 

R  -10 

t  =  _0 _ 

fl  5.5a 

From  equation  (3),  RQ  then  becomes,  RQ  =  118.52  km  .  The  missile 
range  to  lock-on  is  then,  R^, 

RMl  =  VMfcfl  =  59*19  km  • 
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If  the  target  does  not  maneuver  the  time  from  lock- 


on  to  intercept  becomes,  t^* 


**1,0 

«  <VV 

and  the  missile  range  from  lock-on  to  intercept  becomes, 

1^2  =  VMtf2  ~ 

The  total  missile  range  is  then  the  sum  of  the  two, 
R^  =  +  RM2  =  64.64  km  =  34.88  nmiles  . 

As  can  be  seen  from  the  above  analysis,  both  the 
detection  range  of  the  target  and  the  required  missile  range  are 
decreased  significantly  when  an  active  homing  missile  is  used. 

On  the  other  hand,  it  must  also  be  remembered  that  the  complexity 
and  cost  of  the  missile  will  be  increased  as  a  result  of  choosing 
an  active  radar  seeker.  For  the  design  example  in  this  study, 
an  active  radar  seeker  is  chosen;  therefore,  the  maximum  range 
requirement  will  be  35  nmiles  at  a  speed  of  MM  =  3.0,  however, 
this  missile  velocity  is  tentative  until  a  guidance  analysis 
is  complete. 

From  the  preceding  analysis  the  mission  profile  is 
determined.  It  must  be  kept  in  mind  that  the  mission  profile 
may  be  changed  during  the  design  process  to  meet  other  design 
objectives.  The  following  profile  assumes  both  the  target  and 
launch  aircraft  at  the  same  altitude. 
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BOOST 


CRUISE 


TERMINAL 

HOMING 


SEGMENT 

^egin 

h, 

begin 

M  , 
end 

^end 

RANGE 

(1) 

Boost 

1.5 

10,000  ft 

3.0 

10,000 

ft 

(2) 

Cruise 

3.0 

10,000  ft 

3.0 

10,000 

ft 

29.6  nmiles 

(3) 

Terminal 

3.0 

10,000  ft 

3.0 

10,000 

ft 

5.4  nmiles 

III.  GUIDANCE  LAW  SELECTION 


Although  the  specifics  of  the  guidance  system  is  beyond  the 
scope  of  conceptual  design,  the  selection  of  a  guidance  law  is 
necessary  for  initial  calculations.  The  warhead  design  depends 
on  the  expected  miss  distance  between  the  missile  and  the  target 
and  the  lifting  surface  area  depends  on  the  maneuvering  require¬ 
ments  of  the  missile.  Both  the  miss  distance  and  maximum  acceler¬ 
ation  required  are  functions  of  the  missile  guidance  law. 

The  guidance  law  for  a  missile  is  the  analytical  formulation 
used  by  the  guidance  system  to  convert  sensed  target  information 
into  missile  steering  commands.  Three  general  guidance  laws  are 
used.  Most  others  can  be  forced  to  fit  into  one  of  these  cate¬ 
gories.  These  are: 

1)  Pursuit  Guidance 

2)  Line-of-Sight  Guidance 

3)  Proportional  Guidance 

A.  PURSUIT  GUIDANCE 

A  pursuit  guidance  law  is  illustrated  in  F igure  (3-la)  and 
is  one  in  which  the  missile  velocity  vector  is  always  directed 
toward  the  target.  The  target  and  the  missile  velocity  vectors 
must  therefore  be  sensed;  so  this  type  of  guidance  normally 
assumes  an  on-board  tracker.  The  missile  may  have  a  separate 
mid-course  guidance  package  to  increase  range,  but  target  lock- 
on  initiates  the  pursuit  guidance  for  the  terminal  homing  phase. 


Figure  (3-la).  Pursuit  guidance. 


X 


N 


Figure  (3-lb) .  Beam  rider  guidance. 


Figure  (3-lc)  .  Proportional  guidance 
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For  this  reason  it  has  the  advantage  of  launch-and-forget  at 
lock-on.  Since  the  signal  processing  is  limited  to  looking  and 
pointing,  the  avionics  are  relatively  simple  and  usually  on¬ 
board  the  missile.  An  option  for  this  type  of  guidance  would 
be  to  include  a  lead  angle  to  accomodate  faster  moving  targets. 

B.  LINE-OF-SIGHT  GUIDANCE 

Line-of-sight  guidance  is  used  in  a  beam  rider  missile.  This 
guidance  scheme  is  illustrated  in  Figure  (3- lb)  and  requires  that 
the  missile  remain  on  a  line  (beam)  joining  the  target  and  a 
control  point.  The  target  tracker  is  located  at  the  control 
point;  therefore,  avoiding  the  necessity  of  an  on-board  tracker. 
Because  of  this,  a  dedicated  fire  control  system  is  needed  from 
launch  to  intercept.  The  range  of  this  type  of  guidance  is 
normally  less  than  with  the  other  types.  A  speed  advantage  is 
required  for  line-of-sight  guidance  since  no  lead  angle  is  incor¬ 
porated.  The  main  advantage  of  this  type  of  guidance  is  the 
simple  avionics  required  to  maintain  the  missile  in  the  beam. 

C.  PROPORTIONAL  GUIDANCE 

A  proportional  guidance  law  is  one  in  which  the  rate  of  change 
of  the  missile  heading  is  made  proportional  to  the  rate  of  change 
of  the  line-of-sight  between  the  missile  and  the  target.  This 
is  illustrated  in  Figure  (3-lc).  Since  the  guidance  law  antici¬ 
pates  the  target's  future  position,  it  can  attain  a  higher  degree 
of  responsiveness  than  other  guidance  laws.  In  proportional  gui¬ 
dance  the  rate  of  change  of  the  line-of-sight  must  be  sensed 
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on-board  the  missile.  Because  of  this  requirement,  and  the  need 
to  provide  anticipated  steering  commands,  the  avionics  required 
are  the  most  complex  of  the  three  guidance  systems. 


D.  COMPARISON  OF  GUIDANCE  LAWS 

In  early  design  considerations  two  parameters  of  interest 
are  acceleration  required  of  the  missile  and  the  miss  distance 
attainable.  Of  the  three  guidance  laws  only  the  proportional 
law  can  respond  to  fast  maneuvering  targets.  Since  the  missile 
must  stay  in  the  line-of-sight  for  a  beam  rider  system,  any  tar¬ 
get  maneuver  will  cause  large  excursions  in  the  missile  flight 
path,  resulting  in  large  normal  accelerations.  The  pursuit  gui¬ 
dance  law  causes  similar  large  excursions  near  intercept  due 
to  the  velocity  vector  always  pointing  at  the  target. 

Several  system  parameters  affect  the  miss  distance  attain¬ 
able  with  a  particular  guidance  law.  An  excellent  source  on 
the  effects  of  these  parameters  is  an  article  written  by  Dr. 
Robert  Goodstein  [6  ] .  The  parameters  studied  for  their  effect  on 
miss  distance  were: 


1) 

Sensor 

Bias  Angle 

2) 

Noise 

3) 

Target 

Heading 

4) 

Target 

Acceleration 

5) 

Target 

Speed 

6) 

Wind  Gusts 

The  results  have  been  reproduced  and  are  included  in  Figures  (3-2) 
through  (3-7).  Table  3-1  provides  overall  guidance  in  the  selection 


of  a  guidance  law  and  is  also  reproduced  from  the  above  reference. 
A  first  glance  would  indicate  that  proportional  guidance  is  a 
proper  choice  for  all  cases.  It  must  be  kept  in  mind,  though, 
that  cost  and  simplicity  are  also  driving  factors  in  the  design 
process.  Furthermore,  it  can  be  seen  that,  while  proportional 
guidance  with  a  high  gain  has  good  performance  against  maneuvering 
targets,  any  noise  in  the  system  will  highly  degrade  this  perfor¬ 
mance.  For  this  reason  another  guidance  law  may  be  desired,  or 
a  compromise  in  the  gain  selection  may  have  to  be  made  in  which 
some  performance  is  given  up  in  order  to  deal  with  a  noisy  system. 
A  reasonable  range  of  proportionality  constants  that  gives  good 
performance  against  both  maneuvering  targets  and  noisy  systems 
is  k  =  2  to  k  =  6 . 

Once  a  guidance  law  is  selected,  a  more  detailed  analysis  has 
to  be  performed  to  determine  if  the  maximum  acceleration  required 
of  that  particular  guidance  system  is  within  the  attainable 
maneuverability  limits  for  the  missile.  The  maximum  acceleration 
required,  in  turn,  determines  the  lifting  surface  area  needed. 

A  good  estimate  of  maximum  acceleration,  which  keeps  the  miss 
distance  less  than  50  feet,  is  to  set  it  equal  to  three  times 
the  target  acceleration  plus  ten. 

a  =  3a.  +  10 
m  u 

Figure  (3-8)  shows  the  miss  distance  sensitivity  to  target 
acceleration. 
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Figure  (3-3) 
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Figure  (3-4) 


Figure  (3-5) 


NOISE  SENSITIVITY 
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E.  PURSUIT  GUIDANCE  (DETAILED  ANALYSIS) 

As  stated  previously,  a  pursuit  guidance  law  requires  the 
missile  velocity  vector  to  always  point  at  the  target.  For 
this  reason  the  missile  always  ends  up  in  a  tail  chase  situation, 
with  the  maximum  acceleration  occurring  at  the  end  of  the 
encounter.  From  this  description  the  maximum  acceleration  of  the 
missile  can  be  determined. 


Reference 

Direction 


Figure  (3-9) .  Pursuit  geometry. 


From  Figure  (3-9)  the  time  rate  of  change  of  the  range,  R,  is 


also 


or 


R  =  VTcos  8  -  VM 


8  =  -VT  sin  6/R 


R  at  =  VT  sin  6 
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d8 


sin  8 


VT 

“I  dt 


dt  = 


dR 


VT  cos  8 


-  v, 


M 


Substituting  equation  (2)  into  equation  (1)  gives. 


d8  _  VT  dR  _  VM 

slrTT  “  R  (VT  cos  8  -  VM)  ?  Lettin<?  ~  VT 


(cos  8  -  k)  _  dR 
sin  8  d6  "  R 


Integrating  equation  (3)  yields, 


£nR  =  k£n(tan  j)  -  £n(sin  8)  +  2-nc.^ 


P.nR  =  9.  n 


c^  tan^  8/2 


sin  6 


By  trigonometric  identity. 


(tan  S/2)k  =  -^".6) 
(1+cos  8) 


Therefore, 


S.nR 


In 


cx  (sin  8) 
(1+cos  8)* 


k-l' 


and 

c.  (sin  8)k_1 

R  =  — - Z — 

(1+cos  8)k 
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From  the  initial  condition  8  =  6q  when  R  =  R^  , 


C1  = 


R  (1+COSg  r 

o _ 0 


(sinSQ  ) 


k-1 


R  =  R 


(1+cos Bfl 

1+cos  6 


h  (saj 


k-1 


Substituting  the  above  equation  for  R  into  the  equation  for  B 


6  -  - 


V  sinB/R 


8  = 


Vt  /l+cosB 


R  y 1+cos B 


-j 

O' 


(sinB) 


2-k 


(sinBQ  ) 


The  missile  acceleration  can  be  expressed  as  a  normal  and  a 
tangential  component. 


a  =  V  8  n  +  v  t 
M  M 


The  normal  component  is  am,  where, 


m 


=  VM  8 


V  V 
M  T 

am  "  R 


1+cos  B 
1+cos  ^ 


) 


k 


sin  6 

sinB0 


2-k 

1-k 


The  terminal  acceleration  for  a  pursuit  guidance  law  will  occur 
at  the  end  of  the  encounter  (B  *  0)  .  From  the  above  expression 
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the  terminal  acceleration  can  be  evaluated 


a  =  0 ,  for  1  <  k  <  2 

in 

a  =  — ^  T  (  .  2  Q  \  sinS  for  k  =  2 

m  Rq  \  1+cosBq  j  M0 

a  =  °°  for  k  >  2 
m 

Since  pursuit  guidance  always  ends  up  in  a  tail  chase  situation, 
the  missile  will  never  intercept  if  k  <  1  .  Thus,  for  pursuit 
guidance  operating  against  a  non-maneuvering  target,  the  velocity 
ratio  should  be  between  one  and  two.  These  results  indicate 
this  guidance  system  would  not  be  effective  against  air  targets; 
therefore,  results  for  a  maneuvering  target  were  not  pursued. 

F.  LINE-OF-SIGHT  GUIDANCE  (DETAILED  ANALYSIS) 


Figure  (3-10).  Line-of-sight  geometry. 


Figure  (3-10)  illustrates  the  geometry  used  to  derive 
the  beam  rider  equations  of  motion.  The  basic  concept  of  beam 
rider  guidance  is  that  the  missile  is  maintained  in  the  line-of- 
sight  of  the  target  and  a  control  point.  This  can  be  expressed 
in  an  equation  as  follows: 


4>  = 


VT  sin  afc 


V„  sin  a 
M  m 


‘m 


(1) 


where,  r  =  range  from  point  0  to  the  target 
rm  =  range  from  point  0  to  the  missile 

From  equation  (1) 


r.Vw  sin  a  =  r  V_  sin  a. 
t  M  m  m  T  t 


(2) 


As  in  the  case  of  pursuit  guidance,  the  missile  and  target  acceler 
ations  can  be  divided  into  normal  and  tangential  components.  If  j 

the  target  is  limited  to  contant  g  turns,  and  the  normal  component 
of  missile  acceleration  is  of  interest;  then. 


Differentiating  equation  (2)  with  respect  to  time  yields. 


r. VM  sin  a  +  r.  V„  a  cos 
t  M  m  t  M  m 


m 


=  r  Vm 
m  T 


sm  + 


v„ 


afc  cos 


Solving  for  am  , 


a  =  — T7-~ -  [r  V_  sm  a.  +  r  V_  a.  cos  a.  -  r.  V,,  sm  a  ] 

m  t.  V,.cosa  lmT  t  mTt  t  tM  m 

t  M  m 
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From  the  original  figure 


r 


t 


r 


m 


VT  cos  at 

V„  cos  a 
M  m 


Also 


m 


«t  +  ♦ 
%  *  * 


afc  +  $ 

.  * 

a  +  <£ 
m  Y 


The  target  and  missile  accelerations  (normal  components)  become. 


at  =  VT  9t 


a  =  VM  0m 
m  Mm 


Collecting  equations; 


9t  =  at/vt 


♦  =  VT  sin  at/rt 


“t  "  9t  "  ♦ 


rt  =  VT  COS  at 


rm  = 


V„  cos  a 
m  m 


am  = 


r, V  coso  [rm  VT  Sln  a 
t  M  in  c 


t  +  rm  VT  “t  cos  “t 


-  r.  V„  sin  a  1 
t  M  ml 


9m  - 


%  +  ♦ 


a  =  v  0m 
m  Mm 
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The  above  equations  are  the  equations  of  motion  which  describe 
the  target  and  missile  trajectories.  These  equations  cannot  be 
solved  analytically  except  for  highly  specialized  cases.  The  com¬ 
plete  set  of  equations  can  be  solved  using  a  numerical  integra¬ 
tion  technique.  If  Euler's  one  step  method  is  used,  the  algorithm 
is  as  follows ; 

9 1 ( i+1)  =  0 1 (i)  +  At  8 1 ( i ) 

<j>(i+l)  =  <j)(i)  +  At  <j>(i) 


a^,  (i+1) 

rt(i+l) 

r  (i+l) 
m 

ct  (i+l) 
m 

0  (i+l) 
m 


afc(i)  +  At  at(i) 

rfc (i)  +  At  rfc(i) 

r  (i)  +  At  r  (i) 
m  m 

a  (i)  +  At  a_(ij 
m  m 

0  (i)  +  At  0  (t) 
m  m 


With  initial  conditions; 


rt(°)  -  r„ 


r  (0)  *  0 
m 


<M0)  =  <j>0 

o  (0)  =  0 
m 


0.(0)  =  0. 

^0 


0(0) 

rn 


o.(0)  =  a 
c  c0 
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The  target  and  missile  positions  can  be  expressed  as  follows; 


X  (i+D 

m 

-  *m(i) 

+ 

O 

ft 

VM 

COS 

eu) 

m 

ym(i+i> 

m 

-  *m(i> 

+ 

At 

VM 

sin 

9m(i) 

xt (i+l) 

=  xfc(i) 

+ 

At 

VT 

cos 

9fc(i) 

yt(i+D 

=  yt(i) 

+ 

At 

VT 

sin 

0t(i) 

Where 

x  (0)  =  y  (0)  =  0 
m  m 

yt(0)  =  r0  COS  ^0 

yt(o>  ■  r0  sin 

The  above  equations  have  been  programmed  on  the  HP  9830  com¬ 
puter.  Table  3-II  is  a  listing  of  this  program.  The  program 
asks  the  user  for  the  initial  conditions  and  the  target  acceler¬ 
ation.  It  also  asks  for  the  integration  step  increment.  At. 

It  should  be  kept  in  mind  when  using  the  program  that  the  error 
involved  in  integrating  is  of  order  At.  The  output  is  a  plot  of 
missile  and  target  trajectories  as  well  as  the  missile  maximum 
acceleration  and  time  of  flight.  Three  examples  follow  which 
demonstrate  possible  uses  of  the  program.  (Note:  All  angles 
are  input  in  radians.) 
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10  PRINT 
20  PRINT 
30  PRINT 
40  PRINT 


TABLE  3- I I 


"THIS  PROGRAM  DETERMINES  THE  ACCELERATION  OF" 
"A  BEAM  RIDER  MISSILE  AND  PLOTS  THE  TP  A SECTOR 

"INPUT  TIME  INCREMENT  FOR  INTEGRATION" 


50  INPUT  D 1 

60  D I M  NC  250  ] , VC  250  1 , UC 250 ] » VC  250 1 > A 
70  PRINT  "INPUT  INITIAL  TARGET  RANGE 

89  INPUT  R 1 

90  PRINT  “INPUT  TARGET  SPEED" 

100  INPUT  VI 

110  PRINT  "INPUT  MISSILE  SPEED" 

120  INPUT  V2 
ISO  PRINT  "INPUT  INITIAL  LINE  OF  SIGHT  ANGLE' 


At  250  1 5  RC 250 


130  PRINT  "INPUT  INITIAL  LINE  OF  SIGH1 
140  INPUT  PI 

150  PRINT  "INPUT  MISSILE  ALPHA" 

160  INPUT  A2 

178  PRINT  "INPUT  TARGET  ALPHA” 

180  INPUT  A 1 

190  PRINT  "INPUT  MISSILE  THETA" 

200  INPUT  T 2 

210  PRINT  "INPUT  TARGET  THETA" 

220  INPUT  T 1 

230  PRINT  "INPUT  TARGET  ACCELERATION" 

240  INPUT  G1 

250  1=1 

260  R2=0 

265  RC13=R1 

270  At  1  ]=0 

280  MC I  3=0 

290  YC I  3=0 

300  UC  I  3=R1*C0SCP1 > 

310  VC  I  3=Rl*SIN(F'l  > 

320  AC  I  3=8 

330  PRINT  "  XM  YM 

395  PRINT  "DO  YOU  WANT  A  PRINT  OF  THE 

396  INPUT  T9 

397  IF  T9=l  THEN  410 

400  PR  I  NT  XC  I  3j  YC  I  3>  IJC  I  3  5  VC  I  3 

410  1=1+1 

420  T3=G1/V1 

430  P 3 = V 1  * S I N ( A 1 >  •  •  R 1 

440  A3=T3-P3 

450  R3=V 1 *COS (A  1  > 

460  R4=V2*C0S  C  A 2  > 

461  D  6 = R  4  *  V 1  *  S I N  <  A 1 > 

462  D7=R2*V 1* A3* COS C H 1 > 

4 6 3  D 8 = - R 3 * V 2*S I N < A 2  > 

464  D  9 = R 1  *  V  2  *  C  0  S  <  A  2 ) 


THE  OUTPUT 


NT 

0=YESj  1: 


TABLE  3-II  (cont) 


465  fi  4  — B  6 + 0  7  +  B  y )  B  9 

466  T4=A4+F‘3 
470  T1=T1+D1*T3 
480  P 1 =P 1 +B 1 *P3 
490  A1=A1+D1*A3 

495  R 1 =R 1 +B 1 *R3 

496  R2=R2+B 1 *R4 
500  A2=R2+D1*A4 
510  T2=T2+D1*T4 
520  AC  I ]=V2*T4 

530  IF  PBS ( RC I  ] > < ABS < flC 1  —  13 )  THEN  600 
540  M=AC I  ] 

600  T2=fl2+Pl 
610  RC I  ]=R1-R2 

6  2  0  X  C  I  ] = X  C  I  - 1  ] + B 1  *  V  2  *  C  0  S  ( T  2  > 

630  YC I  3*YC I - 1  ]  +  B 1 *V2*S I N ( T2  > 

640  IJC  I  ]=IJ[  I  - 1  ]+B  1  *V  1  *COS  <  T 1  > 

650  VC  I  ]=VC  I - 1 3+D 1  *V  1  *S  I N  <  T 1 ) 

660  IF  RCI]>0  THEN  397 
665  T5=B 1  * (  I  - 1 ) 

678  PRINT  "INPUT  MINIMUM 
680  INPUT  X6 

690  PRINT  "INPUT  MAXIMUM 
7@0  INPUT  X7 

710  PRINT  "INPUT  MINIMUM 
720  INPUT  Y6 

730  PRINT  "INPUT  MAXIMUM 
748  INPUT  Y 7 
8 0 0  8  C ALE  X 6 » X  7  * Y  6 ? Y 7 
810  PRINT  "HAS  AXIS  BEEN 
820  INPUT  FI 
830  IF  F 1=0  THEN  860 
840  XAXIS  0 .*  X7  -10 
850  YAK  IS  0  j  Y  7  1 0 
860  PEN 

870  FOR  W= 1  TO  I 
880  PLOT  XCW3»YCW3 
890  NEXT  W 
900  PEN 

910  FOR  S= 1  TO  I 
920  PLOT  IJC  S  VC  S  3 
930  NEXT  S 
940  PEN 

950  PRINT 

951  PRINT 

952  IF  ACI3>RCI-13  THEN  954 

953  M=AC 1—13 

954  PRINT  "THE  MAXIMUM  ACCELERATION  I 
960  PRINT  "THE  TIME  TO  INTERCEPT  IS 
1000  STOP 


X  VALUE" 
X  VALUE" 

Y  VALUE" 

Y  VALUE" 


BRAUN »  0=YES 


5' 


Example  I  (Non-maneuyering  Crossing  Target) 


\ 

\ 

\ 


\ 

\ - — 

Figure  (3-11) .  Non-maneuvering  crossing  geometry. 


ro 

‘ 

359  meters 

*0 

= 

109.57°  = 

1.9124 

rad 

e. 

17.57°  = 

.3067  rad 

0 

— 

109.57°  = 

1.9124 

rad 

mo 

= 

-91.0°  = 

01. 5882 

rad 

V„  =  373  m/sec 
VT  =  221  m/sec 

Table  3-III  is  the  computer  output.  As  indicated  the  missile 
maximum  acceleration  is, 


a  =  -459.25  m/sec/sec  =  -46.86  g's 


TABLE  3-III 


INPUT  TIME  INCREMENT  FOR  INTEGRATION 

INPUT  INITIAL 

TARGET  RANGE 

INPUT  TARGET  6 

:PEED 

INPUT  MISSILE 

SPEED 

INPUT  INITIAL 

LINE  OF  SIGHT 

ANGLE 

INPUT  MISSILE 

ALPHA 

INPUT  TARGET  ALPHA 

INPUT  MISSILE 

THETA 

INPUT  TARGET  THETA 

INPUT  TARGET  ACCELERATION 

XM 

YM 

XT 

DO  YOU  WANT  A 

P  R I N  T  0  F  T  H  E  0  U  T  P  U  T  ?  0 = Y  E  3  ?  1  = 

0 

0 

-120. 2644734 

-5. 154984716 

17.92341018 

-109. 7301208 

-9. 197499767 

36. 1308181 

— 99.  1 9 576 a  1 y 

-12. 11631881 

54.55019713 

-88. 66141559 

-13. 90888092 

73. 1 1385036 

-78. 127063 

-14.58120543 

91 . 75172791 

-67. 5927104 

-14. 14756982 

1 1 U .  3  9  6  6  y  5  9 

-57.05835782 

-12. 62997166 

1 2  y .  9  y  4  y  y  8  o 

-46. 5240 0  5  2  3 

-18. 0574058 

147.456557 

- 3 5 . 98965264 

-6 . 464996448 

165.7572973 

-25.45530005 

-1 . 893028052 

183.8382156 

-14.92094745 

3.614079320 

201.6565844 

- 4 . 3 8  6594  3  6 3 

10. 08880633 

219,  1760013 

6. 147757728 

17.24154285 

2  3  6 . 3  6  6  4  0  6  2 

16,6821 1032 

25.26133519 

253. 2039305 

27. 21646291 

34. 01774632 

269. 6706819 

37.7508155 

43 . 4596 1 099 

285.7539377 

48. 28516809 

53. 53767559 

301 . 4464555 

58.81952068 

64.2041 194 

316.745131 

69.35387327 

75.41315987 

33 1 . 65083 1 8 

-?  q  _  p;  q  2  2  5  8  7 

87. 12134900 

346. 1677463 

90 . 42257846 

99 . 28777 1 53 

360 . 3028325 

100. 956931 

111.8741571 

374. 0652959 

1 1 1 . 4912836 

124.8449199 

3 8  7 . 4  6  6 1 0  9  7 

122. 0256362 

138. 1671393 

400.5175336 

132.5599888 

151.8104937 

413.232983 

143.0943414 

INPUT  MINIMUM 

X  VALUE 

INPUT  MAXIMUM 

X  VALUE 

INPUT  MINIMUM 

Y  VALUE 

INPUT  MAXIMUM 

Y  VALUE 

HAS  AXIS  BEEN 

DRAWN?  0=YES? 

i=NO 

THE  MAXIMUM  A I 

8CELERATION  IS 

-459. 2526072 

THE  TIME  TO 

INTERCEPT  IS 

1.3  SEC 

256495 
5926479 
9  2  8  8  0  0  7 
2649534 
601 1062 
937259 
2734118 
6895646 
9457 174 
2818702 
618023 
9541758 
2903286 
6  2  6'  4  y  1 4 
S’  £•  y  o  4  d 

k  !  z<  4  9  o  9  y 
97 1 0926 
3  ti  7  2  4  5  4 
6433982 
979551 
3157038 
6518566 
9  y  y  0i  0i  9  4 
3241622 
6  6  0i  3 1 5 


METERS/ SEC 
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y (meters) 


target. 


The  large  acceleration  is  typical  of  line-of-sight  due  to  the 
missile  requirement  to  stay  in  the  beam.  The  trajectories  are 
plotted  in  Figure  (3-12) . 


2.  Example  II  (Effect  of  VM) 

This  example  is  presented  to  study  the  effect  of  missile 
velocity  on  a  crossing  target  (non-maneuvering) . 


Figure  (3-13).  Crossing  target  geometry. 


The  initial  conditions  are  as  follows; 


Tq  =  4000  m 

4>o  =  450  =  .7854  rad 

9  =  45°  =  .7854  rad 

m0 

a.  =  135°  =  2.3562  rad 
t0 

0.  =  180°  =  3.1416  rad 

c0 


VT  =  200  m/sec 


TABLE  3- IV 


THIS  PROGRAM  DETERMINES  THE  ACCELARRTIOH  OF 
A  BEAM  RIDER  MISSILE  AND  PLOTS  THE  TRAJECTORY 

INPUT  TIME  INCREMENT  FOR  INTEGRATION 

INPUT  INITIAL  TARGET  RANGE 

INPUT  TARGET  SPEED 

INPUT  MISSILE  SPEED  «  m/SSC, 

INPUT  INITIAL  LINE  OF  SIGHT  ANGLE 

INPUT  MISSILE  ALPHA 

INPUT  TARGET  ALPHA 

INPUT  MISSILE  THETA 

INPUT  TARGET  THETA 

INPUT  TARGET  ACCELERATION 

KM  YM  KT 

D 0  Y 0 U  W ANT  A  P R  I  N T  OF  THE  0 IJ T P U T  <  0  =: Y ES  »  1  =  N 0 
INPUT  MINIMUM  K  VALUE 

INPUT  MAXIMUM  X  VALUE 

INPUT  MINIMUM  Y  VALUE 

INPUT  MAXIMUM  Y  VALUE 

HAS  AXIS  BEEN  DRAWN ,  0=YES>  l=NO 


THE  M A X I M LI M  ACCELERATI 0 N  I S  5 S .  Q 3 1 2 9 5 4 8  METE R S  S E C /' S E C 

THE  TIME  TO  INTERCEPT  IS  8  SEC 


THIS  PROGRAM  DETERMINES  THE  ACCEL AR AT  I  ON  OF 
A  BEAM  RIDER  MISSILE  AND  PLOTS  THE  TRAJECTORY 

INPUT  TIME  INCREMENT  FOR  INTEGRATION 

INPUT  INITIAL  TARGET  RANGE 

INPUT  TARGET  SPEED 

INPUT  MISSILE  SPEED  n/S^c 

INPUT  INITIAL  LINE  OF  SIGHT  ANGLE 

INPUT  MISSILE  ALPHA 

INPUT  TARGET  ALPHA 

INPUT  MISSILE  THETA 

INPUT  TARGET  THETA 

INPUT  TARGET  ACCELERATION 

KM  YM  KT  YT 

D 0  Y 0 U  M ANT  A  P R I N T  OF  THE  0 IJ T P U T j  O = Y ES*  1=N 0 
INPUT  MINIMUM  X  VALUE 
INPUT  MAXIMUM  X  VALUE 
INPUT  MINIMUM  V  VALUE 
INPUT  MAXIMUM  Y  VALUE 
HAS  AXIS  BEEN  DRAWN >  0=YES>  l=NO 


THE  MAXIMUM  ACCELERATION  IS  72.78774168  METER S/SEC- SEC 

THE  TIME  TO  INTERCEPT  IS  5.6  SEC 
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TABLE  3- IV  (cont) 


THIS  PROGRAM  DETERMINES  THE  ACCEL BRAT  I  ON  OF 
A  BEAM  RIDER  MISSILE  AND  PLOTS  THE  TRAJECTORY 

INPUT  TIME  INCREMENT  FOR  INTEGRATION 
INPUT  INITIAL  TARGET  RANGE 
INPUT  TARGET  SPEED 

INPUT  MISSILE  SPEED  =  SOO  rr>/s£C 

INPUT  INITIAL  LINE  OF  SIGHT  ANGLE 

INPUT  MISSILE  ALPHA 

INPUT  TARGET  ALPHA 

INPUT  MISSILE  THETA 

INPUT  TARGET  THETA 

INPUT  TARGET  ACCELERATION 

XM  YM  NT 

D 0  V 0 U  M ANT  A  P R I N T  OF  THE  0 U T P U T ,  0 = V E 3 » 1 = N 0 
INPUT  MINIMUM  X  VALUE 
INPUT  MAXIMUM  X  VALUE 
INPUT  MINIMUM  Y  VALUE 
INPUT  MAXIMUM  Y  VALUE 
HAS  AXIS  BEEN  DRAWN?  0=YES?  i=NO 


THE  MAXIMUM  ACCELERATION  IS  87.22950371 
THE  TIME  TO  INTERCEPT  IS  4,4  SEC 


METERS-  -SEC/  SEC 
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Figure  (3-14).  Effect  of  V  on  the 


The  program  was  run  three  times  for  a  missile  velocities  of 
400,  600,  and  800  m/sec.  Table  3-IV  contains  the  program  out¬ 
puts  and  Figure  (3-14)  is  the  plot  of  the  trajectories.  As 
can  be  seen  from  the  output,  the  missile  maximum  acceleration 
increases  with  increasing  missile  speed.  For  example, 

VM/VT  =  2  gives  the  smallest  acceleration,  although  the  maximum 

acceleration  for  V  /V_  =  4  is  not  exceedingly  large  for  this 

M  t  ... 

scenario. 


Figure  (3-15).  Maneuvering  target. 


TABLE  3-V 


THIS  PROGRAM  DETERMINES  THE  fiCCEL.HP.HT  I  ON  OF 
fl  BEAM  RIDER  MISSILE  AND  PLOTS  THE  TRAJECTORY 

INPUT  TIME  INCREMENT  FOR  INTEGRATION 

INPUT  INITIAL  TARGET  RANGE 

INPUT  TARGET  SPEED 

INPUT  MISSILE  SPEED 

INPUT  INITIAL  LINE  OF  SIGHT  ANGLE 

INPUT  MISSILE  ALPHA 

INPUT  TARGET  ALPHA 

INPUT  MISSILE  THETA 

INPUT  TARGET  THETA 

INPUT  TARGET  ACCELERATION 

XM  YM  NT 

D  0  V  0  U  W ANT  A  P  R I N  T  0  F  THE  0  U  T  P U  T  ?  0 = Y E  3  ? 1 = N 0 
INPUT  MINIMUM  X  VALUE 
INPUT  MAXIMUM  X  VALUE 
INPUT  MINIMUM  Y  VALUE 
INPUT  MAXIMUM  V  VALUE 
HAS  AXIS  BEEN  DRAWN?  @=YE3?  i=NO 


THE  MAXIMUM  ACCELERATION  -IS  -355.8215593 
THE  TIME  TO  INTERCEPT  IS  5.8  SEC 


YT 


METE  R  S  •••■'  S  E  C  S  E  C 
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(3-16) .  Effect  of  a  7g  maneuver  upon  the  encounter. 


The  initial  conditions  are; 


rQ  =  10000  m  VT  =  821.436  m/sec 

4>0  =  90°  =  1.5708  rad  VM  =  985.723  m/s ec 

9  =  90°  =  1.5708  rad 

% 


9.  =  -90° 

0 


a^  -  68.60  m/sec/sec 

From  the  output  (Table  3-V) ,  notice  the  large  missile  acceler¬ 
ation  required  (211.55  m/sec/sec)  to  intercept  a  maneuvering 
target.  The  trajectories  are  plotted  in  Figure  (3-16) . 

G.  PROPORTIONAL  GUIDANCE  (DETAILED  ANALYSIS) 


Figure  (3-17) .  Proportional  guidance  geometry. 
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Proportional  guidance  automatically  establishes  a  lead  angle 
and  reacts  to  a  changing  line-of-sight.  The  basic  guidance  law 
used  equates  the  rate  of  change  of  the  missile  heading  to  a 
constant  times  the  rate  of  change  of  the  line-of-sight.  From 
the  above  figure  this  law  can  be  expressed  as. 


0  =  k  a 

m 

From  the  figure,  the  rate  of  change  of  the  line-of-sight,  a,  is 
given  by, 

V_  sin  0  -  V  sin  6 

•  _  T  t  M  m 


As  in  the  case  of  pursuit  and  line-of-sight  guidance,  the  para¬ 
meter  of  interest  here  is  the  normal  acceleration?  therefore,  the 
missile  and  target  tangential  accelerations  are  assumed  to  be 
zero.  In  this  case. 


a  =  VM  8 
m  Mm 


at  *  VT  6t 


From  the  guidance  law. 


and 


a  =  VM  k  c 
m  M 


amr  =  VM  k(VT  sln  8t  -  VM  sin  Bm> 
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Since  vm  =  vt  *  0/  the  time  derivative  of  this  equation  is, 


r  a  +  r  a  =  k  V  f(V  0  cos  0.  -  V  8  cos  8  )  1  (1) 

m  m  mtt  tmm  m 


From  the  original  figure. 


6=6+0 
m  m 


m 


6  =  9  -  0=  9  -  , 

mm  m  k 


Also, 


8  =  0m(l~) 

m  m  k 


et  =  0t  -  a 


ft  -  a _ £ 

6t  -  9t  -  k 


Making  these  substitutions,  equation  (1)  becomes, 


.  e 


r  a  =  -a  r  +  k,  V  f  V,,,  (0.  -J£>  cos  8^  -  VM0  (l-±)cos  6m  ] 
m  m  'M  Ttk  t  Mmk  m 


Since, 


am  -  VM  8m 


at  "  VT  6t 


r  a  = 
m 


-a„  r  +  k  VM  a.  cos  6.  -  k  VM  am  cos 
m  -  M  t  t  Mm  m 


-  am  (VT  cos  8t  -  VM  cos  8m> 
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From  Figure  (3-17) , 


i  ■  VT  cos  6t  '  VM  cos  6m 

•  • 

r  a  =  -2  r  +  k  VM  a.  cos  0.  -  k  V„  a  cos  0 
m  m  M  t  t  M  m  m 


k  v  a  cos  0  a 

a  - - r - f  [2r  +  k  VMcos0m  ] 


m 


m 


Collecting  equations; 


r  =  VT  cos  St  -  VM  cos  Sm 


•  •  • 

=  0fc  -  or 


*  *  » 

6m  -  em  -  ° 

J  .  VTslnBt  -  Vin6m 


*  at 

=  ^ 


9  =  k  a 

m 


The  above  equations  are  the  equations  of  motion  for  a  missile 
using  proportional  navigation  assuming  constant  missile  and 
target  speeds.  As  with  the  line-of-sight  equations,  the  motion 
is  quite  complex.  The  equations  cannot  be  solved  analytically 
except  for  special  cases.  One  such  case  will  be  investigated 
here.  That  is  for  a  non-maneuvering  crossing  target. 
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Example  IV  (Non-Maneuvering  Crossing  Target) 


If 


Figure  (3-18).  Non-maneuvering  target. 


For  this  case,  r  = 
a 


-VM  COS  B 
M  m 


a  =  [2r 

m  r 


-  kr] 


=  £  (it  —  2 ) 


m 


£na  =  (k-2)  Unr  +  Jhc. 
m  l 


a„  =  a.  at  r  =  r. 
m  0  0 


a  =  a  (— ) 
m  0  r0 


k-2 


From  this  equation, 


from  equation  (1) 


if  k  >  2 


a 


0  as  r  0 


if  k  =  2 


constant 


a  =  a.  = 
m  0 


if  k  <  2  a„  — *  «  as  r  0 
m 


From  the  above  example  it  can  be  realized  that  the  pro 
portionality  constant,  k,  must  be  greater  than  two.  A  more 
general  analysis  of  the  equations  of  motion  can  be  obtained  by 
solving  the  equations  numerically.  The  same  Euler's  one  step 
method  is  used  here  with  initial  conditions,  at  t  =  0, 


a  (0) 
m 

r  (0) 

et(o) 

Sm(0) 

o  (0) 

9 1  ( 0 ) 

0(0) 

m 


The  algorithm  used  is  as  follows: 

r(i+l)  =  r(i)  +  At  r(i) 

a_(i+l)  =  a  (i)  +  At  a  (i) 
in  in  m 

cr(i+l)  =  a(i)  +  At  a(i) 

0  (i+1)  =  0_ (i)  +  At  0  (i) 
mm  m 


9 t (i+1)  =  et(i)  +  At  et (i) 

S_ (i+1)  =  8_  (i)  +  At  B  (i) 
mm  m 

Bt(i+1)  =  Bt(i)  +  At  Bt(i) 

The  trajectory  of  the  missile  and  target  can  be  deter¬ 
mined  assuming  the  missile  is  at  the  origin  at  t  =  0.  From  the 
original  figure  the  missile  and  target  positions  are  given  by, 

x  (i+1)  =  x  (i)  +  At  V  cos  6(i) 
m  m  Mm 

y  (i+1)  =  y  (i)  +  At  V  sin  0  (i) 
m  m  M  m 

xfc (i+1)  =  yfc (i)  +  At  VT  cos  6fc(i) 

yfc (i+1)  =  yfc(i)  +  At  vT  sin  9t(i) 

With  initial  conditions, 

x  (0)  =  y  (0)  =  0 
m  m 

xfc(0)  =  rQ  cos  a0 

yt(0)  =  rQ  sin  aQ 

Evaluating  the  missile  initial  acceleration  can  be  more  complicated 

One  procedure  which  is  both  realistic  and  of  interest  is  to  have 

the  missile  and  target  on  a  constant  bearing  -  decreasing  range 

course  (afc  =  0)  when  the  target  initiates  a  constant  g  turn  at 

t  «  0.  In  this  case  a  (0)  =0  and  the  subsequent  motion  can  be 

m 

found. 
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TABLE  3 -VI 


10  PRINT  "THIS  PROGRAM  FINDS  THE  MAXIMUM " 

20  PRINT  "MISSILE  ACCELERATION  FOR  A" 

3  0  P  PINT  "  P  R  0  P  0  R  T 1 0  N  A  L.  N  A  V I G  A  T  1 0  N  S  Y  S  T  EM" 

33  PRINT 

3 6  D I M  X C  250  3 j Y [ 2 5 0  3  > U C 2 5 0  3 » V [250  1 
38  DIM  AC  250  3 ?  TC  250  3  ?  RC 259  1 
40  PRINT  "INPUT  TIME  INCREMENT" 

50  INPUT  D1 

60  PRINT  "INPUT  NAVIGATION  CONSTANT" 

70  INPUT  K 

80  PRINT  “INPUT  MISSILE  VELOCITY" 

90  INPUT  VI 

100  PRINT  "INPUT  TARGET  VELOCITY" 

110  INPUT  V2 

120  PRINT  "INPUT  TARGET  ACCELERATION" 

130  INPUT  A 1 

140  PRINT  "INPUT  INITIAL  MISSILE  ACCELERATION" 

150  INPUT  AC  1  1 

160  PRINT  "INPUT  INITIAL  RANGE" 

17@  INPUT  R[ 13 

180  PRINT  "INPUT  BETA  TARGET" 

200  INPUT  B 1 

210  PRINT  "INPUT  BETA  MISSILE- 
220  INPUT  B2 

230  PRINT  "INPUT  THETA  TARGET- 
240  INPUT  T 1 

250.  PRINT  "INPUT  THETA  MISSILE" 

260  INPUT  T2 

270  PRINT  "INPUT  SIGMA" 

280  INPUT  SI 
290  1=1 

295  PRINT  "  MISSILE  POSIT  TARGET  POSIT 

296  PRINT 

297  PRINT  "  XI  Y 1  X2 

300  XC I  3=0 

310  YC I  3=0 

320  IJC  I  3=RC  1  3*COS<Sl> 

330  VC  I ]=R[  1  3+S I H  <  S 1  > 

340  PRINT  "DO  YOU  WANT  A  PRINT  OF  THE  OUTPUT >  @=YESO l=NO" 

345  INPUT  Q1 

350  IF  Q 1 = 1  THEN  365 

360  PRINT  KC  I  If  VC  I  3,  UC  I  3i  VC  I  1 

365  1=1+1 

370  XC I ]*XC I - 1 3+V l*COS C T2 ) *D 1 
380  YC I ]=YC I - 1 3+V 1*D1*SI N <  T2 > 

390  UC I 3=UC 1-1 3+V2+D 1 *COS < T 1 > 

400  VC  I  3= VC  I - 1 3+D 1 +V2+S I N  < T1 > 

410  R  2 = V  2  *  C 0  S  < B 1 > - V 1 * C  0  S  < B  2 > 

420  A  8 = K *  V 1  *  A 1  *  C  0  S  <  B 1 > / R  C I - 1 3 
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TABLE  3 -VI  (cont) 


430  A9=AC  I  - 1  ]*  ( 2*R2+K*V  1  *C0S «.  B2  > ) /RE  1  —  13 

440  A3=AS-R9 

445  TC I  ]=( 1-1 >*B1 

450  S  £ = < V £  *  S I N < B 1 > - V 1  *  S I N ( B  £  >  > - '  R  C I - 1  3 
460  T8=Hl  ■-•"■,'£ 

470  T9=K*S2 
4S0  B3=T8-S2 
490  B4=T9-S£ 

500  RC  I  3*RC  1-1  3+Dl*R2 

518  RC  I  JafiC  1-1  3+DUA3 

5£0  S1=S1+D1*S2 

530  T1=T1+B1*T8 

540  T2=T£+B1*T9 

550  6 1=B 1+01*63 

560  B2=6£+B1*64 

570  IF  fit  I KfiC 1  —  13  THEN  608 

580  Z 1 =ftC I  3 

600  IF  RC I  3 > 8  THEN  358 

601  IF  RES (AC I  ]><RBS<flC 1-1  3>  THEN  693 
692  Z 1 =HC 1-13 

603  Z2=TC I  1 

4  nn  r  lit 


618 

PRINT 

619 

PRINT 

628 

PRINT 

“MAXIMUM  MISSILE  ACC 

IS" 

’21" 

1  METERS/SE 

640 

PRINT 

"MISSILE  TIME  OF  FLIGHT 

IS 

"22"  SEC 

700 

PRINT 

"INPUT  MINIMUM 

VALUE 

OF 

y  " 

719 

INPUT 

X6 

720 

PRINT 

"INPUT  MAXIMUM 

VALUE 

OF 

730 

INPUT 

i  i 

r 

740 

PRINT 

"INPUT  MINIMUM 

VALUE 

OF 

V " 

759 

INPUT 

Y6 

760 

PRINT 

"INPUT  MAXIMUM 

VALUE 

OF 

Y " 

770 

INPUT 

Y  7 

780 

SCALE 

N6 ? N7  ?  Y6 ’ Y7 

790 

PRINT 

"HAVE  THE  AXIS 

BEEN 

DRAWN  j 

0--YE3<  l=NO 

791 

INPUT 

FI 

792 

IF  F 1  = 

■0  THEN  838 

800 

NANIS 

0*  X7.-10 

810 

YANIS 

0  5 Y7/10 

830 

PEN 

840 

FOR  J  = 

=  1  TO  I 

850 

PLOT  > 

:c  j  3  >  yc  j  i 

860  NEXT  J 

879  PEN 

880  FOR  H=1  TO  I 
899  PLOT  UCM3»VCW3 
909  NEXT  U 

919  PEN 
1000  STOP 
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m 


The  above  algorithm  was  programmed  on  the  HP  9830  computer. 
The  listing  is  included  in  Table  3-VI.  The  inputs  and  output 
of  the  program  are  the  same  as  the  line-of-sight  guidance  pro¬ 
gram.  Several  examples  follow  which  demonstrates  the  use  of  the 
program. 


Figure  (3-19) .  Crossing  maneuvering  target. 

VM  =  373  m/sec  9.  =  17.57°  =  0.3067  rad 

M  tQ 

V  =  221  m/sec  8.  =  -91.0°  =  -1.5882  rad 

c0 

a.  =  156.8  m/sec/sec  8  =  -38.0°  =  -0.6632  rad 

t  mQ 

r.  =  359  m  a  =  109.57°  =  1.9124  rad 

0  o 

0  =  71.57°  *  1.2491  rad 

m0 

The  output  is  listed  in  Table  3-VII.  An  important  aspect  of  this 
problem  is  the  maximum  acceleration  required  (202.92  m/sec/sec). 
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TABLE  3- VI I 


THIS  PROGRAM  FINDS  THE  MAXIMUM 
MISSILE  ACCELERATION  FOR  A 
PROPORTIONAL  NAVIGATION  SYSTEM 

INPUT  TIME  INCREMENT 

INPUT  NAVIGATION  CONSTANT 

INPUT  MISSILE  VELOCITY 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 


INPUT  SIGMA 

MISSILE  F'O 

SIT 

TARGET  F'O: 

I:  IT 

XI 

Y 1 

Y  £. 

DO  YOU  WANT  A 

0 

P R I N T  OF  THE  0 U T P U T  ?  0 = Y ES  j 1 = N 0 
0  -i 20. 2644734 

338. 256495 

5 . 8  9  6  6  U  8  8  0  2 

1 7 . 69326330 

-109. 7301208 

341 . 5926479 

1 1 . 72925942 

y  5  >  4  y  7  7*  &  S  y  fr 

-99.32072172 

345.3003306 

17. 50036490 

5 3 . 14238894 

-89. 04937488 

y  4  y  a  y  r'  4  y  \  ( i 

23.21064114 

i'  ti .  y  9  b  b  9 1 1  u 

-78 . 92900524 

353. 81 1 1620 

2  8 . 8  5  8  6  0 4  7  3 

88.67091421 

- 6 8 . 9 7 234776 

358. 603601 

34. 44056825 

1  y  6  ■ 4659749 

-59. 19193143 

3 63. 7461642 

39 . 95057982 

124.2834458 

-49 , 60006345 

369. 2323885 

45.38038997 

142. 1255214 

-40.20881377 

375.0553462 

50.71944940 

159. 994963 

-31 . 02999988 

331.2077340 

55 . 95494 1 06 

177.8950225 

-22.07517195 

387 . 68 1 8020 

6 1 . 07 1 84968 

195.82934 

-13.3555983 

394. 4694  0 3  6 

66. 053071 1 

213.8018172 

-4.882251242 

401 . 5619976 

70. 87956373 

231 . 8164631 

3 . 3  3 4  2  U  6794 

408 . 9506598 

75.53054295 

249. 8772187 

1  1 . 2834366 

4 1 1*  ■  b  y  b  y  y  y  y 

79.98371892 

267 . 9877 60S 

18.95543527 

424.578633 

84.21557618 

286. 1512922 

26.34054872 

432.7982802 

88 . 20 1 69238 

304.370331 

33.  42948387 

441.2746885 

91 . 91709127 

322. 64649  8 9 

40.21 332036 

449.9971918 

95 . 33662288 

340 . 9803284 

46.68352174 

458.9548139 

98.43536132 

359 . 37 1 8956 

52.8319462 

468. 136283 

101 . 1890086 

377 . 8 166893 

58 . 65085689 

477.5300458 

1 0  3  •  5  7  4  d  y  ‘5  ‘5 

396.3135251 

64. 13293155 

487. 1242815 

105.5693244 

414.8565113 

69.27127182 

496. 906917 

107. 1539546 

433. 4390689 

74.05941 184 

5  0  6 . 8  6  5  6  4  2  6 

108.3100167 

452. 0532039 

7 8 . 49132646 

516. 9879266 

109.0215346 

470. 6  8 9626 3 

82. 5614  3 y  7  8 

527. 2610316 

109.2748269 

489 . 3379062 

86. 26462715 

537. 672  0  3  0  4 

109.0585142 

507.9866517 

89.59623167 

548 . 2078223 

108.3634271 

526 . 6236943 

92. 55  2 0  6 

558.8551496 

107. 1824200 

545. 2362632 

95. 12839268 

569, 6  0 06142 

105.5101054 

563.811135 

97.32198776 

5  8  0 . 4  3  0  6946 

103.3425418 

582.3347461 

99. 13008494 

591.3317626 

100.6769169 

600. 7932657 

1 00 . 550409 

602. 290 1 0 1 0 
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TABLE  3- VI I  (cont) 


MAXIMUM  MISSILE  ACC  IS  202.9212389  METER 

MISSILE  TIME  OF  FLIGHT  IS  1.7  SEC 

INPUT  MINIMUM  VALUE  OF  X 

INPUT  MAXIMUM  VALUE  OF  X 

INPUT  MINIMUM  VALUE  OF  Y 

INPUT  MAXIMUM  VALUE  OF  Y 

HAVE  THE  A X I S  BEE N  D R A U N >  0 = Y ES»  1  =  N 0 
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Figure  (3-20) .  Trajectories  of  crossing,  maneuvering  intercept. 


The  same  problem  was  run  for  line-of-sight  guidance  with  no 
target  acceleration  in  Example  I.  The  maximum  acceleration 
was  459.25  m/sec/sec.  This  points  out  the  advantage  of  pro¬ 
portional  navigation  over  line-of-sight  guidance.  Figure  (3-20) 
is  a  plot  of  the  trajectories. 

3.  Example  VI  (Effect  of  k) 

This  example  demonstrates  the  effect  of  varying  the  pro¬ 
portionality  constant,  k.  The  scenario  is  as  follows. 


/  ^ 


Figure  (3-21) .  Initial  geometry. 


Vfc  =  208m/sec 


V  =  413  m/sec 


afc  =  68.60  m/sec/sec 

8.  =  -106°  =  -1.85  rad 

t0 

8  =  -29°  =  -.5061  rad 

rao 
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TABLE  3-VIII 


INPUT  TINE  INCREMENT 

I N P U T  N A V I G R T 1 0 N  C 0 N 3 TRNT  =  3 

INPUT  MISSILE  VELOCITY 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE  POSIT  TARGET  POSIT 

XI  Y 1  X2 

D 0  Y 0 U  W ANT  A  P R I H T  0 F  T H E  0 U T P IJ T<  O = Y E S  ?  1=N 0 


MAXIMUM  MISSILE  ACC  IS  83.70521471  METERS/ SEC/ SEC 

MISSILE  TIME  OF  FLIGHT  IS  3.2  SEC 

INPUT  MINIMUM  VALUE  OF  X 

INPUT  MAXIMUM  VALUE  OF  X 

INPUT  MINIMUM  VALUE  OF  Y 

INPUT  MAXIMUM  VALUE  OF  Y 

H A V E  THE  A X I S  BEE N  D R A M N  ,  0 = Y ES» 1= N 0 


THIS  PROGRAM  FINDS  THE  MAXIMUM 
MISSILE  ACCELERATION  FOR  A 
P R 0 P 0 R T 1 0 N A L  HA V I G A T 1 0 N  S Y S T E M 

INPUT  TIME  INCREMENT 

INPUT  NAVIGATION  CONSTANT  *  V 

INPUT  MISSILE  VELOCITY 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE  POSIT  TARGET  POSIT 

XI  Y1  X2 

D O  Y O U  M ANT  A  P R I N T  OF  THE  0 U T P U J,  0 = Y ES> 1= N 0 


M A X I M U M  M I S S I L E  A C C  IS  71.1  1 3 3 3 197  METER S / S E C / S E C 
MISSILE  TIME  OF  FLIGHT  IS  3.2  SEC 
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TABLE  3-VIII  (cont) 


THIS  PROGRAM  FINDS  THE  MAXIMUM 
MISSILE  ACCELERATION  FOR  A 
P R 0 P 0 R T 1 0 H A L  HA V I G A T 1 0 N  S V S T E M 

INPUT  TIME  INCREMENT 
I N P U T  N A V I  GAT  10 h  C 0 N S T  A  N  T  “  S’ 

INPUT  MISSILE  VELOCITY 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE  POSIT  TARGET  POSIT 

XI  Y 1  X2 

D 0  Y 0 U  U ANT  A  P R I N T  OF  THE  0 U T P IJ T  >  O = Y E S »  1  =  N 0 


MAXIMUM  MISSILE  ACC  IS  65.38577248  METERS/SEC/SEC 

MISSILE  TIME  OF  FLIGHT  IS  3.2  SEC 

INPUT  MINIMUM  VALUE  OF  X 

INPUT  MAXIMUM  VALUE  OF  X 

INPUT  MINIMUM  VALUE  OF  Y 

INPUT  MAXIMUM  VALUE  OF  Y 

HAVE  THE  A  X I S  B E  E  N  D  R A  W N >  0 = Y  E  S  ? 1 = N  0 
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meters) 1200 


Figure  (3-22).  Effect  of  varying  k. 


.8483  4ad 


8  =  48.6°  = 
m0 

9  =  -28.4°  =  -.4956  rad 
=  77.6°  =  1.3544  rad 

The  problem  was  run  for  k  =  3,4,  and  5. 

Table  3-VIII  is  the  output.  It  can  be  seen  that  the 
effect  of  increasing  the  proportionality  constant  is  to  decrease 
the  maximum  acceleration  required.  Figure  (3-22)  is  a  plot  of 
the  trajectories. 

H.  DESIGN  EXAMPLE  (GUIDANCE  LAW  SELECTION) 

From  the  examples  given  in  this  chapter,  it  can  be  seen  that 
a  missile  designed  to  encounter  a  highly  maneuverable  target, 
such  as  a  fighter,  requires  a  proportional  guidance  law  to  limit 
the  maximum  acceleration  required  of  the  missile.  To  select 
the  proportionality  constant,  it  was  assumed  that  the  threat 
could  maintain  a  constant  7g  turn  at  =  1.5  and  an  altitude 
of  10,000  feet. 

Three  cases  were  investigated,  (1)  A  head-on  encounter  with 
the  target  initiating  a  turn  at  10,000  meters  range,  (2)  A 
crossing  encounter  in  which  the  target  turns  into  the  missile  at 
10,000  meters  range,  and  (3)  An  oblique,  closing  encounter  in 
which  the  target  turns  into  the  missile.  The  scenario  and 
computer  outputs  are  shown  in  Figures  (3-23),  (3-24),  and  (3-25). 
From  this  analysis  the  crossing  encounter  requires  the  largest 
acceleration  (263.2  m/sec/sec). 

For  the  crossing  scenario  then,  the  missile  speed  was  varied 
and  the  results  indicated  that  as  the  speed  increased  the  maximum 
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7  g  turn 


7  =  821  m/sec  (M  =  2.5) 

m 

7fc  *  492  m/sec  (M  =  1.5) 
rQ  =  10,000  m 


180° 


THI*  PROGRAM  FINDS  THE  MAXIMUM 
MISSILE  ACCELERATION  FOR  A 
PR  0  P  0  R  T 1 0  N  A  L  N  A  V I G  A  T 1 0  N  S  V  S  T  E  M 


INPUT  TIME  INCREMENT 

I N  P IJ  T  N  A  V I G  A  T 1 0  N  C  0  N  S  T  A  N  T 

INPUT  MISSILE  VELOCITY 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE  POSIT  TARGET  POSIT 

X 1  V 1  !::!  2 

D 0  Y 0 U  U ANT  A  PR  I N T  OF  THE  0 U T P IJ T ,  0 = Y E S »  1  =  N 0 


MAXIMUM  MISSILE  ACC  IS  126.1911994  METERS  ••'SEC 
MISSILE  TIME  OF  FLIGHT  IS  9  SEC 

INPUT  MINIMUM  VALUE  OF  X 


Figure  (3-23).  Head-on  scenario 


THIS  PROGRAM  FINDS  THE  MAXIMUM 
MISSILE  ACCELERATION  FOR  A 
PROPORT I ONAL  NA V I  GAT  I  OH  SYSTEM 


INPUT  TIME  INCREMENT 

INPUT  NAVIGATION  CONSTANT 

INPUT  MISSILE  VELOCITY 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE  POSIT  TARGET  POSIT 

X 1  Y 1  K  ^ 

D 0  Y 0 U  U ANT  A  P R I N T  OF  THE  0 U T P IJ T »  Q = Y E S  t  1  =  N 0 


MAXIMUM  MISSILE  ACC  IS  263.1995227  METERS--' SEC-" 

MISSILE  TIME  OF  FLIGHT  IS  3.3  SEC 
INPUT  MINIMUM  VALUE  OF  X- 

Figure  (3-24).  Crossing  scenario. 
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V 


492  m/sec 

821  m/sec 
135° 

53.9° 

10,000  m 


0  =  135° 

0 

0  =  53.9° 

m0 

a0  "*  0 
k  =  4 


THIS  PROGRAM  FINDS  THE  MAXIMUM 
MISSILE  ACCELERATION  FOR  A 
PROP 0 R T 1 0 N A L  HA V I G A T 1 0 N  S V S T E M 


INPUT  TIME  INCREMENT 

INPUT  NAVIGATION  CONSTANT 

INPUT  MISSILE  VELOCITY 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE  POSIT  TARGET  POSIT 

XI  ■  Y 1  X2 

D 0  Y 0 U  U ANT  A  P R I N T  OF  THE  0 U T P IJ T  >  Q  =  Y ES»  1=N 0 


M A X I M U M  MISSILE  A C C  IS- 156. 97959 6 2  METER S S E C / S E C 

MISSILE  TIME  OF  FLIGHT  IS  3.8  SEC 
INPUT  MINIMUM  VALUE  OF  X 


Figure  (3-25).  Oblique  scenario. 
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TABLE  3 -IX 


INPUT  TIME  INCREMENT 
INPUT  NAVIGATION  CONSTANT 

I N  P IJ  T  M I S  S I L  E  V  E  L  0  C I T  V - -  *  2 .  O 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE  POSIT  TARGET  POSIT 

XI  VI  X2 

B 0  Y 0 IJ  N ANT  A  P R I H T  OF  THE  0 U T P U T »  O = ' < E S »  1=N 0 


MAXIMUM  MISSILE  ACC  IS  275.5734145  METERS/SEC/SEC 

MISSILE  TIME  OF  FLIGHT  IS  10.3  SEC 

INPUT  MINIMUM  VALUE  OF  X 


THIS  PROGRAM  FINDS  THE  MAXIMUM 
MISSILE  ACCELERATION  FOR  A 
P  R  0  P  0  R  T 1 0  N  A  L  N  A  V I G  A  T 1 0  N  S  Y  S  T  E  M 

INPUT  TIME  INCREMENT 

INPUT  NAVIGATION  CONSTANT  ^  _  _ 

INPUT  MISSILE  VELOCITY  - >  -  3.0 

INPUT  TARGET  VELOCITY 

INPUT  TARGET  ACCELERATION 

INPUT  INITIAL  MISSILE  ACCELERATION 

INPUT  INITIAL  RANGE 

INPUT  BETA  TARGET 

INPUT  BETA  MISSILE 

INPUT  THETA  TARGET 

INPUT  THETA  MISSILE 

INPUT  SIGMA 

MISSILE  POSIT  TARGET  POSIT 

XI  Y 1  X2 

DO  YOU  WANT  A  PRINT  OF  THE  OUTPUT ?  0=YESj1=NO 


MAXIMUM  MISSILE  ACC  IS  248.233583  METERS- "SEC. -SEC 

MISSILE  TIME  OF  FLIGHT  IS  7.8  SEC 
INPUT  MINIMUM  VALUE  OF  X 
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Figure  (3-26) .  Selection  of  Navigation  Constant. 


acceleration  decreased.  The  results  for  =  2.0  and  ^  =  3.0 
are  shown  in  Table  3-IX.  This  result  indicates  the  desirablity 
of  retaining  the  missile  Mach  number  originally  selected. 

For  the  crossing  case  and  a  missile  speed  of  =  3.0,  the 
proportionality  constant  was  then  varied  from  k  =  2  to  k  =  6. 

The  results  are  plotted  in  Figure  (3-26) .  If  the  maximum 

acceleration  is  limited  to  31  g's  (3afc  +  10),  the  required  pro¬ 

portionality  constant  is  k  =  3.75.  This  is  well  within  the 
desirable  range  of  2  -  6  indicated  earlier. 

From  this  analysis  the  required  performance  objectives  are: 

M  =  3.0 
m 

k  =  3.75 


(a) 

m 


=  31  g's 


max 


IV.  SIZING  THE  DIAMETER 


The  missile  diameter  is  determined  by  one  of  three  driving 
factors.  For  relatively  short  range  missiles  the  diameter  will 
be  fixed  by  either  the  warhead  or  the  seeker  requirements.  As 
might  be  expected,  for  longer  range  missiles  the  diameter  will 
more  likely  be  fixed  by  either  the  warhead  or  the  seeker  require¬ 
ments.  As  might  be  expected,  for  longer  range  missiles  the 
diameter  will  more  likely  be  fixed  by  the  propulsion  require¬ 
ments  in  order  to  prevent  excessive  propulsion  system  lengths. 

An  initial  estimate  of  the  missile  diameter  must  be  made  at 
this  point  in  order  to  proceed  with  the  design.  The  initial 
seeker  requirement  can  be  determined  from  a  knowledge  of  the 
lock-on  range  requirement  found  in  Chapter  2.  The  warhead 
necessary  to  inflict  a  "kill"  can  also  be  estimated  from  infor¬ 
mation  about  the  target  and  characteristic  explosives.  The 
propulsion  requirement  cannot  be  determined  because  of  the  lack 
of  any  aerodynamic  drag  or  weight  information  at  this  point.  For 
this  reason  the  missile  diameter  will  be  now  sized  for  seeker 
or  warhead  requirements.  The  missile  propulsion  requirements 
will  be  determined  later  in  the  design  process,  and  it  may  be 
necessary  at  that  time  to  resize  the  missile  to  meet  these  pro¬ 
pulsion  requirements. 

Selection  of  the  type  of  seeker  depends  upon  the  operational 
arena  of  the  missile.  The  seeker  of  a  shoulder  fired,  battle¬ 
field  missile  would  not  be  the  optimum  seeker  of  a  shipboard 
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missile  where  antenna  and  component  sizes  are  not  limiting 
factors.  All  types  of  guidance  use  some  portion  of  the  electro¬ 
magnetic  spectrum.  The  three  primary  areas  of  use  are  the 
electro-optical,  infrared,  and  radio  frequencies.  The  milli¬ 
meter  wave  section  of  the  spectrum  is  also  of  current  interest 
in  the  design  of  missiles  due  to  small  component  size  and  will 
also  be  discussed.  The  following  table  lists  some  of  the  major 
advantages  and  disadvantages  of  the  three. 


Optical 


Infrared 


RF 


Advantages 

Target  resolution  (de¬ 
tail) 

Real  time  information 
Three  dimension  effect 

Improved  resolution 
over  RF 

Longest  range 

Least  absorption  and 

attenuation 


Disadvantages 

Bad  weather  degrades 
Night  use  degrades 

Attenuation  due  to 
aerosols  and  atmosphere 

Larger  components 


A.  THE  RADAR  RANGE  EQUATION 

An  omnidirectional  antenna  is  one  that  radiates  power  in  all 
directions  equally.  If  the  power  radiated  by  an  antenna  is  Pfc, 
the  power  density  at  a  distance  from  the  source  is  given  by, 

2 

Power  density  =  Pt/(4ir  Rfc  ) 

2 

4tt  Rfc  =  area  of  a  sphere  of  radius  Rfc 


Since  antennas  are  normally  directive  instead  of  omnidirectional, 
most  of  the  power  is  radiated  in  a  particular  direction.  The 


gain,  Gfc,  is  a  measure  of  the  increased  power  from  a  directive 
antenna  as  compared  to  an  omnidirectional  antenna.  Therefore  the 
power  density  from  a  directive  antenna  can  be  expressed  as. 


Power  density  = 


P  G 
t  t 


4  it  R, 


This  is  the  power  density  which  arrives  at  the  target.  The  tar¬ 
get  intercepts  a  portion  of  this  energy  and  reradiates  it  in  the 
opposite  direction.  The  radar  cross  section,  a,  is  a  measure  of 
the  effective  area  of  the  target.  The  power  radiated  by  the 
target  is  PechQ, 

Pt  Gt  a 

Pecho  4*  Rt2 

This  energy  propagates  as  if  it  were  radiated  by  an  omni¬ 
directional  antenna.  Therefore  if  the  receiving  antenna  is  a 
distance,  Rr,  away  the  power  density  at  the  receiver  is 

Pfc  G  a 

(Power  density)  =  - 5 - 

r  (4ir  Rt  j  ( 4 tt  Rr*) 


If  the  energy  is  intercepted  by  the  receiving  antenna,  which  has 
an  effective  area  as  seen  by  the  returned  energy  of  Ar;  then  the 
power  received  by  the  radar  ,  Pr,  is 
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(1) 


P.  G  0  A 
t  t  r 

(4ir  Rt2)  (4tt  Rr2 


) 


This  is  the  simplest  form  of  the  radar  equation  and  can  be  used 
to  determine  the  size  of  antenna  required. 

B.  ACTIVE  RADAR  HOMING 

Active  homing  is  the  method  of  missile  guidance  in  which  the 
radar  transmitter  and  receiver  are  located  on-board  the  missile. 
In  this  case  the  same  antenna  is  used  for  both  transmitting  and 
receiving.  The  radar  equation  then  becomes, 


P 

r 


pt  Gt  °  At 
(4ir  R2)2 


where,  R,  =  R  =  R  and  A„  ■  A.  . 
t  r  r  t 

The  minimum  power  for  which  the  target  can  be  detected,  pm^n*  is 
a  function  of  many  variables.  A  full  development  of  this  term 
can  be  found  in  reference  (8) . 


s 

P  .  =  k  T-  B  F  (=2.) 

mm  0  n  n  N.  . 

0  min 


A  *)  _ 

Boltzmans  constant  is  k  =  1.38  x  10  joule/  k.  The  value  of 

-  21 

kTQ  at  room  temperature  is  4  x  10  watt/cps  of  bandwidth.  The 

bandwidth,  B  ,  noise  figure,  F  ,  and  minimum  signal  to  noise 
n  n 

ratio,  (Sq/N  }  are  all  functions  of  the  receiver.  Typical 
0  0  min 
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values  are  listed  below. 


B  *  1  MHz 
n 

Fn  -  7.5  db  =  2.37  (for  crystal  mixer) 

(Sq/N.)  =  14.7  db  =  5.43  (for  probability  of  detection, 

min 

PD  =  0.9  and  probability  of  false  alarm, 

P^a  =  1/15  minutes) 

The  above  values  give  P  .  =  5.15  x  lO*"14  watts.  This  is  the 

mm 

value  which  will  be  used  throughout  this  section. 

From  antenna  theory  the  gain  is  related  to  the  effective 
antenna  area  by. 


G  = 


4  fr  A 


The  maximum  radar  range  can  then  be  shown  to  be 


R 


max 


V4 


Pt  J  At 

4tt  X2  P  , 
mm 


(1) 


Equation  (1)  for  radar  range  does  not  include  any  system  losses 
It  also  does  not  include  the  statistical  nature  of  several  of 
the  parameters.  Because  of  these  assumptions  the  actual  range 
of  a  radar  may  be  as  small  as  one-half  of  what  the  radar  range 
equation  predicts  for  laboratory  conditions.  For  this  reason, 
twice  the  required  range  should  be  used  when  using  the  above 
equation. 
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1.  Example 


The  AN/APQ-153  is  the  airborne  attack  radar  system  used 
on  the  F-5E  aircraft.  The  following  parameters  apply  to  this 
radar. 


fQ  =  8-10  GH2 

X  =  c/f  =  m/sec  =  >0333  m 

u  9  x  10  /sec 

c  =  speed  of  electromagnetic  propagation 
=  3  x  108  m/sec 
P  =  80  kW 

Antenna  =  Parabolic  dish  30.5  x  40.6  cm 
Afc  =  0.12383  m2 

for  a  target  of  1  square  meter  of  radar  cross-section, 


R  =  36.16  km  -  19.5  nmiles 
max 

For  the  case  of  an  active  homing  radar,  the  size  of  the 
transmitting  and  receiving  antenna  is  the  parameter  of  interest. 
The  antenna  size  may  very  well  drive  the  design  diameter  of  the 
missile.  The  antenna  diameter  can  be  expressed  as  follows. 


[4tt  X4  P  . 
_ mi 

Pt 


2  „  R  41 

min  max 


1/2 


Let 


then 


'4lt  x2  pmln  1 

'  l  Pt  a  J 


7Tdt  2 

A,  =  =  C  R 

t  4  max 


1/2 


(2) 


4C  R. 


max 


7T 


therefore,  d.  =  2  R 


ma  xw  it 
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RANGE  (NMILES) 


Equation  (2)  has  been  plotted  in  Figure  (4-1)  for  various 
values  of  the  transmitter  frequency.  From  this  plot  and  a  know¬ 
ledge  of  the  maximum  lock-on  range  required,  the  antenna  size 
can  be  determined. 

From  equation  (2)  and  Figure  (4-1)  there  are  two  obvious 
ways  to  decrease  the  antenna  size  required.  (1)  Increasing  fre¬ 
quency  is  the  best  way  to  reduce  antenna  and  electronic  component 
sizes.  The  current  trend  is  toward  higher  frequencies.  (Milli¬ 
meter  waves.)  One  problem  is  that  the  equations  developed  in 
this  section  do  not  include  atmospheric  attenuation.  For  fre¬ 
quencies  above  about  30  GHz  the  absoprtion  due  to  atmospheric 
gases  increases.  This  is  shown  on  Figure  (4-2) .  As  indicated 
on  this  figure  there  are  "windows"  where  the  attenuation  is  less. 
These  "windows"  occur  at  frequencies  of  34  GHz,  94  GHz,  140  GHz, 
and  220  GHz.  These  are  the  frequencies  where  most  of  the  cur¬ 
rent  research  and  development  is  going  on.  As  the  frequency  in¬ 
creases,  the  wavelength  approaches  the  size  of  rain  droplets. 

For  this  reason,  radar  performance  is  greatly  reduced  in  inclement 
weather.  (2)  Increasing  transmitter  power  will  also  decrease 
the  size  of  antenna  necessary.  The  limiting  factor  in  this  area 
is  the  lack  of  high  power  sources.  In  the  millimeter  range  the 
available  power  from  current  traveling  wave  tubes  is  50-100 
watts.  Increasing  power  is  obviously  confined  to  size  and  weight 
limitations  of  missile  components. 


1 


1 
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C.  SEMI-ACTIVE  HOMING 

The  advantage  of  semi-active  homing  is  obvious  when  the  radar 
range  equation  is  investigated.  From  equation  (1) , 


_  2 
Rt  Rr 


P,  G.  a  A 
t  t  r 

(  4  7T  )  2P  . 

min 


In  the  above  equation  the  missile  range  from  the  target  is  R^. 
The  transmitting  and  receiving  antennas  are  at  different  ranges 
and  have  different  characteristics  in  this  case.  As  before, 

4-rr  a. 


therefore. 


R  2R  2  a  PtAtAi- 
t  r  ,2 


4ir  A  P 


min 


The  main  advantage  is  in  the  transmitter  characteristics.  Since 
the  transmitter  is  not  located  in  the  missile,  it  is  not  normally 
limited  in  size  and  weight  requirements.  In  the  above  equation 
if  a  transmitter  power  and  standoff  range,  Rfc,  is  chosen  the 
receiving  antenna  can  be  sized  for  a  maximum  homing  range,  R  , 

lUcLX 

of  the  missile. 

1.  Example 

Pt  =  100  km 
Rfc  =  100  nmiles 
fQ  =  10  GHz 
Afc  =  4  m 

P  .  =  5.15  x  10"14  W 

min 
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R  (nmiles) 

max 

dr  (in) 

10 

U1 

00 

20 

11.6 

50 

29.1 

It  can  be  seen  from  comparing  these  numbers  to  those  of  Figure 
(4-1)  that  the  required  antenna  size  is  less  than  one  half  of 
that  required  for  an  active  homing  radar  of  the  same  frequency. 

D.  DESIGN  EXAMPLE  (ANTENNA  SIZING) 

An  active  radar  was  assumed  in  Chapter  2  to  decrease  the 
required  missile  range.  The  lock-on  range  was  10,000  m  or 
5.4  nmiles.  As  stated  earlier,  twice  this  number  should  be  used 
for  determining  antenna  size.  From  Figure  (1)  for  a  range  of 
10.8  nmiles,  and  a  transmitter  power  of  10  kW  at  fQ  =  20  GHz, 
the  required  antenna  size  is  d  *  10  inches. 


E .  INFRARED  SEEKERS 

In  the  design  of  missile  seekers  two  parameters  of  primary 
importance  are  range  and  size.  The  idealized  range  for  an  infra¬ 
red  tracker  relates  these  two  factors.  The  idealized  range  is 
the  range  at  which  the  signal-to-noise  ratio  is  unity  and  is 
given  by, 


”d*t  Ttt>  d  2  j 
a  IR  a 


4  v/A  f  A  • 


1/2 


(1) 
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Specific  detectivity 

Transmission  through 
atmosphere 

Transmission  through 
IR  optics 

Aperture  diameter 

Radiant  intensity 

Receiver  bandwidth 

Sensitive  area  of 
detector 


D* 


Typical  value 


10 


10 


T  0-1.0 

a 

tir  °-1-0 

D 

/  103 

Af  103 

A.  10-6  -  10-1  cm2 

a 


The  derivation  of  equation  (1)  and  its  use  are  the  subjects  of 
this  section.  Some  references  (9)  may  give  the  above  equation 
in  terms  of  the  Noise  Equivalent  Intensity,  NEI . 


where  \J Af  A. 

NEI  =  D*  A  T  TIR 

ol  cl 


1.  Planck's  Law 

The  radiant  emittance  of  a  body  is  a  measure  of  the  radiant 
power  per  unit  area  emitted  from  the  surface. 

W  =  — ^  watt/cm2  (2) 

The  spectral  radiant  emittance  is  the  radiant  emittance  per  unit 
wavelength  interval. 
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«TWP 


Wx  *  |y  watt/cm2  y 
y  =  micron  =  10~^  meters. 

Planck's  law  gives  the  blackbody  spectral  radiant  emittance  as 
a  function  of  wavelength  and  temperature. 


2 

x  _  2  TT  h  c _ 1 _  n 

X  BB  "  x5  exp (hc/XkT) -1  v 

-34 

h  =  Planck's  constant  =  6.6238  x  10  Joule-sec 

g 

c  =  speed  of  light  =  3  x  10  m/sec 
X  =  wavelength 

23  o 

k  =  Boltzmann's  constant  =  1.38  x  10~  Joule/  K 
T  =  Absolute  temperature,  °K 


Figure  (4-3)  shows  equation  (3)  for  various  absolute  tempera¬ 
tures.  As  can  be  seen  the  wavelength  at  which  maximum  radiant 
emittance  occurs  varies  with  temperature.  This  maximum  occurs 
at  a  wavelength  given  by  Wien's  displacement  law,  ^max 

X  T  =  2897. 8y  °K 
max 


2.  Emissivity 

Actual  bodies  do  not  emit  radiation  according  to  Planck's 
law.  A  more  typical  plot  of  radiant  emittance  is  shown  in 
Figure  (4-4).  Spectral  emissivity,  ex,  is  defined  as  the  ratio 
of  the  actual  spectral  radiant  emittance  to  the  blackbody  spectral 
radiant  emittance, 


W, 


eX  = 


(Wx) 


BB 
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Figure  (4-4) .  Actual  spectral  radiant  emittance. 


As  shown  in  the  above  figure,  e^,  may  vary  with  wavelength.  A 
grey  body  is  defined  as  one  which  has  constant  spectral  emissivity, 

e  =  *  constant 

IR  systems  normally  use  filters  to  limit  the  accepted 
radiation  to  a  specific  wavelength  band.  The  radiant  emittance 
of  a  body  between  wavelengths  X^  and  X2  becomes, 

x2 

W  =  /  e(W, )  dX 

Xx  A  BB 


In  Figure  (4-5)  the  surface  at  the  orgin  emits  a  total 
energy  WA  into  a  hemisphere  normal  to  A. 

The  radiance  is  defined  as  the  radiant  power  per  unit 
solid  angle  per  unit  projected  area, 


Radiance  =  N 


32P  _1 _ 3 

cos63A3ft  ~  cos9  3ft 


3P 

3A 
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Figure  (4-5) .  Radiant  emittance  of  a  body. 

From  equation  (2) 

„  1  3W  .  .  .  .  . .  2 

N  =  ■  -Q  watts/steradian  cm 

COS0  Oio 

The  radiant  intensity  is  defined  as  the  radiant  power  per  unit 
solid  angle  from  a  point  source. 

9P 

Radiant  intensity  =  J  =  watt/steradian 

The  following  is  a  summary  of  the  definitions  of  radiant  energy 
quantities , 

2 

W  =  radiant  emittance,  watts/cm 

J  =  radiant  intensity,  watts/steradian 

2 

N  =  radiance,  watts/steradian-cm 


From  the  definitions  the  relationship  between  W,  J  and  N  are. 


W  =  ttN  =  it  J/A 

From  the  above  definitions  and  Figure  (4-5)  the  energy  into  solid 
angle,  ft.  is, 

NA  ft.  cos  6  =  —ft,  A  cos  0 

1  7T  1 

The  energy  into  solid  angle  is  given  by, 

NAft,  =  §  ft-  A 
2  TT  2 


3 .  Energy  into  a  Hemisphere 

From  the  definition  of  radiance  and  the  radiant  emittance, 
the  radiance  in  terms  of  the  radiant  emittance  can  be  found. 


N 

W 

N 


cos03A3ft 


3P 

3A 


Ji _  3W 

cos0  TH" 


Figure  (4-6)  . 


z 


Energy  emitted  into  a  hemisphere. 
108 


From  Figure  (4-6)  the  incremental  solid  angle  is 


r  sin  6  d  tp  r  dd 

- 2 - 

r 


sin  8  d  9  d  $ 


If  the  above  surface,  A^,  is  considered  a  Lambertian  surface, 
the  radiance,  N,  is  independent  of  the  direction  of  radiation. 


dW  =  cos  8  N  dQ 


dW  =*  cos  8  N  sin  8  d  0  d  (ji 

The  total  radiant  emittance  into  a  hemisphere  above  the  surface 
is  then, 

2tt  -tt/2 

w  /  /  N  cos  8  sin  8  d  8  d  4> 

o  o 

fl  .  2  1 

W  *  2tt  N  sin4  9J 


W  =  it  N 

4.  Targets 

Infrared  targets  include  a  wide  variety  of  radiation 
sources.  The  radiance  of  most  bodies  can  be  divided  into  that 
due  to  self-emission  and  that  due  to  reflection  of  incident 
radiation. 


N  =  N  +  N 
e  r 

The  relative  magnitude  of  these  contributions  depends  on  a  number 
of  factors  and  varies  from  target  to  target  and  operating 


environment. 


a.  Self-emission. 


Self-emission,  also  referred  to  as  thermal  emission, 
depends  primarily  on  the  temperature  of  the  body  and  the  emis- 
sivity.  The  most  often  used  appraoch  is  to  consider  the  body  as 
a  grey  body  which  emits  radiation  according  to  the  Stefan- 
Boltzmann  law. 

4 

Ne  =  ~ tt -  watts/steradian-cni 

where  <i  =  5.67  x  10  ^  watts/cm^  (°K)  ^ 

This  term  is  the  total  radiance  (over  all  wavelengths)  and  is 
not  the  same  as  used  previously, 
b.  Reflection. 

Radiance  due  to  reflectance  depends  on  the  illumi¬ 
nating  source.  This  source  may  be  the  sun,  active,  or  semi¬ 
active  sources.  It  is  obvious  that  at  night  for  a  passive 
infrared  system,  the  radiance  due  to  reflection  is  not  a  con¬ 
tributing  factor.  For  this  reason  only  the  radiance  due  to  self¬ 
emission  is  considered  in  this  section. 

5.  Target  Temperature 

Since  the  self-emittance  of  a  target  depends  on  the 
temperature  of  the  target,  a  method  for  determining  this  tem¬ 
perature  is  needed.  The  temperature  of  an  aerial  target  varies 
depending  on  the  aspect  of  the  target.  The  propulsion  system 
has  hot  surfaces  such  as  the  nozzle  and  exhaust  plumes.  There 
may  also  be  hot  surfaces  due  to  aerodynamic  heating  and/or 
solar  radiation. 


Air  breathing  engines  normally  have  exhaust  plumes 
ranging  from  600  to  1000°K.  Rockets  typically  have  much  hotter 
plumes.  The  flame  temperatures  for  liquid  propellants  range 
from  2500  to  7500°K.  Solid  propellants  flame  temperatures  range 
from  1700  to  3500°K.  The  plume  temperature  can  be  estimated 


from  the  relation. 


T 

flame 

T 

plume 


1  +  Me2 


where 


T^  =  stagnation  temperature 
Tq  =  static  temperature 

M  =  Mach  number  at  the  nozzle  exit 
e 


a.  Example  I 

For  a  flame  temperature  of  2700°K  and  an  exit  Mach 
number  of  3.0  the  plume  temperature  can  be  found. 


Tplume 


T flame 


l+Ir«e 


1270. 6°K 


For  many  missile  encounters  the  exhaust  plume  may 
be  shielded  from  the  infrared  sensor.  For  a  head-on  encounter 
the  temperature  of  interest  is  the  skin  temperature  of  the  target. 
This  temperature  is  due  to  aerodynamic  heating  and  is  a  function 
of  the  target  speed  and  the  target  material.  One  approach  to 
finding  this  temperature  is  through  the  use  of  the  recovery 
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factor,  which  requires  some  knowledge  of  the  material  of  the 
target.  The  recovery  factor,  r,  of  a  material  is  defined  as 
follows : 


T  -  T 

surface  ambient 

T  .  -  T 

stagnation  ambient 


The  skin  temperature  of  the  target  then  becomes. 


T  =  T  +  r  (T  —  T) 

surface  ambient  stag  amb 


The  stagnation  temperature  is  found  from  the  relation. 


T  . 
stag 


T  . 
amb 


(1  + 


M2) 


The  Mach  number,  M,  is  that  of  the  target  and  the  specific  heat 
ratio,  y,  is  for  air. 

b.  Example  II 

A  target  flying  at  M  =  2.5  where  the  ambient  temper¬ 
ature  is  300°K  has  a  recovery  factor  of  0.75. 


T  = 

stag 

Tsurface 


575°K 


=  581°K 


This  is  the  temperature  used,  along  with  the  emissivity  of  the 
target,  to  find  the  radiant  emittance  of  the  target  from 
equation  (4) . 


LU 


Figure  (4-7) .  Simple  IR  system. 

Figure  (4-7) ,  above,  shows  a  simple  IR  system  and  a 
target  at  a  range  R.  If  the  system  is  sensitive  to  radiation 
in  the  3  to  5  micron  region  the  radiant  emittance  becomes , 


The  above  integral  is  best  evaluated  on  the  computer,  if  the 
ambient  temperature  and  target  speed  is  known,  the  skin  temper 
ature  of  the  target  can  be  determined. 


113 


The  radiance  from  the  target  then  becomes , 


If  Figure  (4-7) ,  the  solid  angle  of  the  aperture  as  seen  from 

A 

the  target  is,  ^  _a 

R2 

The  power  seen  at  the  detector  surface  is  then, 

IR  Power  =  N  Q 


The  above  formula  assumes  no  attenuation  by  the  atmosphere  or 
the  IR  system  optics.  This  attenuation  is  significant  in  actual 
IR  systems.  These  factors  are  normally  accounted  for  through 
the  use  of  atmospheric  and  IR  optics  transmission  coefficients. 


T  =  transmission  of  the  atmosphere 

cl 


T 


IR 


transmission  of  the  IR  optics 


The  total  power  at  the  detector  then  becomes. 


P 


T 

a 


tir  At 


irR 


A  W 
a 


(5) 


7.  Detectors 

Detectors  are  devices  which  are  radiation  transducers. 
It's  purpose  is  to  change  the  incoming  radiant  power  to  an 
electrical  signal,  which  can  then  be  amplified.  Detectors  can 
be  divided  into  two  main  categories.  (1)  Thermal  detectors  - 
The  responsive  element  of  a  thermal  detector  is  sensitive  to 
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temperature  changes  brought  about  by  the  incident  radiation. 

(2)  Photodectectors  -  Responsive  elements  of  photodetctors  are 
sensitive  to  the  number  of  incident  photons. 

Detectors  also  are  made  up  of  windows,  apertures  and 
Dewar  flasks.  The  window  restricts  the  bandwidth  to  which  the 
detector  is  sensitive.  The  aperture  may  limit  the  field  of  vie- 
in  order  to  limit  photon  noise.  The  Dewar  flask  cools  the 
detector  which  improves  the  detectivity. 

Detectivity  of  a  detector  is  defined  as, 

D  _  signal/noise  _  S/N  ,g, 

input  power  ~  P 

The  specific  detectivity  is 

D*  -  D  [ Af  Ad]  1/2  (7) 

Af  =  Bandwidth 

Ad  =  Sensitive  area  of  detector 

For  a  tracking  system  the  bandwidth  is  that  of  the  preamplifier 
in  Figure  (4-7) .  The  input  to  the  preamplifier  is  proportional 
to  the  incoming  IR  energy,  which  has  been  modulated  to  give  tar¬ 
get  resolution  from  the  background  and  provide  line-of-sight 
information. 

A  simple  chopper  is  shown  in  Figure  (4-8)  .  It  consists 
of  an  opaque  material  which  has  a  wedge  cut-out  of  angle  a.  The 
rotation  causes  the  input  from  a  point  source  to  be  modulated. 
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Figure  (4-8) .  Simple  Chopper. 


while  that  of  the  background  is  not.  The  input  signal  to  the 
preamplifier  would  look  like  Figure  (4-9) .  The  frequency  content 
of  the  signal  in  Figure  (4-9)  can  be  found  from  a  Fourier  Analysis 
If  the  pulses  are  assumed  to  be  sinusoidal  of  period  T,  the 
optimum  bandwidth  is,  , 

“  -  I 


T  = 


W 


TTCtW_ 


- - - » 


Figure  (4-9) .  Preamp  input. 
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The  specific  detectivity  is  characteristic  of  the 
detector  used.  Reference  [9]  is  an  excellent  source  of  infor¬ 
mation  on  operational  detectors. 

8 .  Idealized  Range 
Equation  (5)  is, 


P  = 


Ta  TIR  At  Aa  W 
7TR2 


From  equations  (6) and  (7)  this  becomes. 


D* 


S/Njr  R 


[Af  Ad] 1/2  Ta  TIR  At  Aa  W 


The  idealized  range,  where  the  signal  to  noise  ratio  is  unity,  is 

At  ; 

T7T 


0  D*  T  Ttt,  A  W 
2  a  IR  t  a 


[Af  Ad]  '  it 


To  simplify  this  equation  the  radiant  intensity  is  given  by. 


J  * 


At  W 


Since  the  parameter  of  interest  in  missile  design  is  the  aperture 

"D  a2 

diameters  Aa,  it  is  replaced  with  — so  that. 


R0  * 


D  = 
a 


D*  Ta  TIR  *  Da  J 


1/2 


4  v/Af  A, 


V  4/^a7 

D*  "T  TIR  it J 


n  1/2 
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a.  Example  (Idealized  Range) 

From  example  II  a  target  flying  at  M  =  2.5  has  a  skin 

temperature  of  581°K.  If  this  target  has  a  presented  area  of 
,  2 

1  m  ,  the  detector  size  needed  to  detect  the  target  at  a  range 
of  10,000  m  can  be  determined.  From  Wien's  displacement  law  the 
maximum  radiation  occurs  at, 


X 


max 


4.99  y 


If  the  system  is  designed  to  accept  radiation  from  3  to  5  microns, 
and  the  emissivity  of  the  target  is  0.7. 

5  2 

W  =  e  2*hc  1-  -  dX 

,5  exp (hc/XkT) -1 
3 

If  this  equation  is  integrated  on  the  computer,  the  radiant 
emittance  becomes, 

W  *  1410  watts/m 

The  radiant  intensity  becomes, 

A  W 

J  =  - =448.82  watts 

TT 

Typical  values  of  the  parameters  in  the  idealized  range  equation 
are  (9) 

T  =0.75 
a 

TIR  "  °*95 
Af  =  1000  cps 

A^  =  1  cm2 

D*  =  1  x  10^°  (cps)1^2  cm/watt 
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Substituting  these  values  into  the  idealized  range  equation  gives, 

D  =  0.0355  m  =  1.4  in. 
a 

As  can  be  seen  the  size  of  the  IR  seeker  is  relatively  small 
compared  to  other  seekers.  The  above  analysis  is  for  the 
"idealized"  range.  Attenuation  of  IR  radiation  can  be  quite 
high  thereby  increasing  the  seeker  size  required. 

F.  WARHEAD  SIZING 

The  conditional  kill  probability  of  a  missile  is  the  prob¬ 
ability  that  the  target  is  destroyed  given  that  the  warhead  is 
delivered  to  a  point  in  space  and  the  fuze  detonates  the  war¬ 
head  at  a  miss  distance  r. 


Figure  (4-10) .  Encounter  geometry. 
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The  fragment  distribution,  D  (<{>),  is  the  number  of  fragments  per 
unit  solid  angle,  ft.  The  total  number  of  fragments  within  the 
cone  is , 


N  =  D(<j»)ft 

The  fragment  density,  p,  is  the  number  of  fragments  per  unit  of 
area  normal  to  the  path. 

_  N  _  D(4>)ft 


however, 


and 


a  -  S 

r 


p  -  p(|)  _  DU)  sin2<)> 

r  rm 

m 


(8) 


1.  Target  Vulnerability 

The  vulnerability  of  aircraft  or  missile  components  is 
normally  determined  experimentally.  Fragments  of  a  specified 
size  are  fired  at  the  component,  and  the  damage  is  assessed  to 
determine  if  the  fragment  would  cause  a  kill.  If  a  large  num¬ 
ber  of  fragments  are  fired,  the  ratio  of  killing  fragments  to 
hits  can  be  determined.  This  ratio  is  defined  as  the  probability 
that  given  a  hit  a  kill  will  result,  =  \  #  Assuming  the 

K/n  Njj 

distribution  of  hits  is  uniform  over  the  target, 


Ap  =  presented  area  of  target 


=  vulnerable  area  of  target 

As  would  be  expected,  P^^,  depends  greatly  upon  the 
encounter  geometry.  It  will  be  dependent  upon  the  aspect  of 
the  aircraft,  and  also  depends  upon  the  type  of  kill  specified. 
If  the  target  is  assumed  to  be  a  spherical  target,  the  prob¬ 
ability  of  kill  given  a  hit  can  be  assumed  constant. 

2 .  Conditional  Kill  Probability 

The  fragment  density  is  given  by  equation  (8).  From 
this  density  expression,  the  average  number  of  hits,  a,  on  the 
vulnerable  area  of  the  target  is, 

2 

_  _  D  (<j))  sin  <j>  * 

ci  —  A  A 

^7 

m 

It  is  customary  to  assume  that  the  distribution  of  hits  on  the 
presented  area  follows  a  Poisson  distribution.  The  conditional 
kill  probability  then  becomes, 


-a 


PD  =  1  -  e 


PD  =  1  -  exp 


2 

-D  (<j>)  sin  4>  Av 


m 


(10) 


The  above  expression  depends  upon  the  fragment  distribution, 
D(<f>).  If  the  warhead  casing  is  scored  such  that  it  produces  N 
fragments  of  uniform  size  and  mass,  m,  the  problem  is  simplified 
by  formulating  an  alternate  expression  for  a. 


Figure  (4-11) .  Static  encounter. 


The  warhead- target  encounter  geometry  is  shown  in 
Figure  (4-11) .  As  shown  the  area  of  the  fragment  ring  depends 
on  the  miss  distance,  r  .  The  area  of  the  ring  is  given  by,  A 

A,  =  2 it  r  (i  +  2r  tan  6) 
f  m  w  m 

If  the  N  fragments  are  distributed  evenly  in  A^,  the  number  of 
fragments  per  unit  area  is 

_N  _  _  N _ 

p  ”  A*  2ir  r  +  2r  tan  8) 
f  m  w  m 

The  average  number  of  hits  on  a  vulnerable  component  is  then, 


From  equation  (9),  Av  =  ?K/H  Ap 


Therefore,  a  =  p  Ap 


N  PK/H  Ap 

a  “  2ir  r  71  +  2r  tan  6) 


m  w 


m 


If  in  Figure  (4-11) ,  the  average  target  width  into  the  paper  is 
given  by  W,  the  presented  area  is, 

A_  =  W  ( Jl  +  2r  tan  3) 
p  w  m 


The  average  number  of  hits  is  then. 


al  = 


N  PK/H  w 
27r  rm 


(11) 


The  distance  at  which  the  target  just  fills  the  fragment  ring 
is  the  critical  miss  distance,  and  the  maximum  distance  for 
which  applies.  This  is  found  by  setting 

L  =  H  +  2r  tan  3 
w  c 

From  which 

L  -  Z 

r  =  _ w 

c  2  tan  3 


If  the  miss  distance  in  Figure  (4-11)  is  such  that  the  entire 
target  is  always  presented  to  the  fragment  ring;  i.e.,  rm  >  rc, 
the  presented  area  becomes. 


A  =  W  L 
P 
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The  average  number  of  hits  then  becomes,  a 


The  parameter  of  interest  in  this  chapter  is  the  diameter 
of  the  warhead  required  to  achieve  a  specified  kill  probability. 
The  development  thus  far  is  the  conditional  probability  of  kill, 
PD.  It  has  been  assumed  that  the  guidance  system  delivers  the 
warhead  to  the  point  of  interest,  and  the  fuze  detonates  the 
warhead  at  this  point.  Since  only  the  conditional  probability 
of  kill  is  determined,  the  purpose  of  this  section  will  be  to 
maximize  PD  or  to  find  the  warhead  diameter  which  sets  P  =  1. 

From  the  threat  to  be  encountered,  the  target  presented 
area,  Ap,  and  the  vulnerability,  can  be  determined.  Also 

from  an  alaysis  of  the  threat,  the  siz  e  and  impact  velocity  of 
the  fragments  necessary  to  kill  the  target  can  be  determined. 

The  initial  velocity  required  to  obtain  the  impact  velocity  is 
a  function  of  the  explosive  used  and  the  charge  to  mass  ratio, 
C/M. 


C  Mass  of  explosive/unit  length 


M  Mass  of  warhead  casing/unit  lengt 


The  initial  velocity  is  given  by  Gurney's  equation, 


Gurney's  constant,  ,/2E for  various  explosives  are  given  below  {11 ] . 


f 


Explosive 

Density,  kg/m 

V2E,  m/s  ec 

TNT 

1590 

2316.5 

RDX 

1650 

2834.6 

HMX 

1840 

3118.1 

PETN 

1730 

2834.6 

Tetryl 

1620 

2500.0 

Composition  B 

1680 

2682.2 

Octol 

1800 

2895.6 

From  equation  (13)  , 

the  charge  to  mass 

ratio  necessary  to  attain 

the  specified  initial  velocity  can  be  determined. 


c  \2'2* 

M  1-Vi2/2(2E) 


Figure  (4-12).  Warhead. 
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With  the  casing  thickness  fixed  as  a  function  of  warhead  radius 
the  number  of  fragments,  N,  can  be  determined.  In  order  to 
achieve  the  desired  velocities  of  fragments,  it  is  essential 
to  have  enough  length  for  the  diameter.  An  acceptable  length 
to  diameter  ratio  for  a  cylindrical  warhead  is  from  two  to  three 
Most  air-to-air  missiles  have  a  length  to  diameter  ratio  of  2.5. 
For  this  analysis  a  value  of  2.5  is  used. 

The  warhead  casing  volume  is  then  given  by,  V  , 


V 

c 


2irr  ti 
w 


w 


lOirtr  2 
w 


The  mass  of  the  casing  is,  mc, 

m  =  10-rrtr  2p 
c  w  c 


The  total  number  of  fragments,  N,  is  obtained  by  dividing  the 
case  mass  by  the  individual  fragment  mass. 


N 


2 

m  lOrtr  p 

c  _  _ w  c 

m  m 


From  the  equation  (11)  or  (12)  the  average  number  of  hits  can 
be  determined  as  a  function  of  warhead  radius 


ai  ” 


5ttw  °cPK/HW 


mr 


m 


a 


II 


5trw  PcPK/HWL 

mr™  U +  2r  tang) 
m  w  m 
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mm 


The  conditional  probability  of  kill  can  be  determined  now  for  a 
given  warhead  radius  using  equation  (10) . 

a.  Design  Example  (Effect  of  radius  on  Pp) 

If  the  above  equations  are  programmed  for  the  con¬ 
ditions  listed  below,  a  plot  of  the  conditional  probability  of 
kill  versus  warhead  radius  can  be  obtained. 

Threat:  RAM-K 

L  =  19.51  m 
w  =  6.10  m 

PK/H  =  -10 
Fragments : 

m  =  105  grains  =  0.0068  kg 
=  2133.6  m/sec 
Pc  =  7000  kg/m 
8  =  20  degrees 


Explosive:  Composition  B 

p  =  1680  kg/m 


Miss  distance:  r  =  50  ft  =  15.24  m 

m 


Figure  (4-13)  is  a  plot  of  the  output.  From  this 
figure  it  can  be  seen  that  a  warhead  radius  of  r^  =  .06  m  is 
required  to  achieve  a  conditional  kill  probability,  PD  =  1.0. 
Therefore  the  missile  diameter  required  for  warhead  considera¬ 
tions  is,  d  =  4.72  inches. 

The  warhead  radius  also  varies  with  the  required 
initial  velocity.  From  the  equations  for  charge  to  mass  ratio 
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.06 

(meters) 


and  casing  thickness  the  radius  required  to  achieve  a  specified 
kill  probability  can  be  determined. 


Since, 


t  = 


1 


and 


C 

M 


Vi2/2E 

1  -  Vi2/2(2E) 


the  average  number  of  hits  becomes. 


al  - 


5pc  PK/H  W  rw3 


v7V2E 


a  -  ViV2(2E )/  pe 


+  1 


mr 


m 


The  above  equation  assumes  r  <  r 

me 


Letting 


b  = 


5  c  PK/H  w 

‘  /  2  \ 

.  /Vi_/2E  ) 

pc 

—  +  3 

.  \l  -  Vi2/2(2E)/ 

T  X 

pe 

mr 


m 


the  probability  of  kill  becomes, 

2 


-br 


PD  =  l-e 


w 


If  the  conditional  probability  of  kill  is  selected  as 
equation  (15)  can  be  solved  for  the  warhead  radius. 

130 


(15) 

=  0.999, 


^6.9078  ^ 1/3 


(16) 


The  initial  velocity,  ,  required  to  achieve  a  target  kill  is 
a  function  of  the  miss  distance,  fragment  size,  expected  en¬ 
counter  altitude  and  target  characteristics.  For  this  simple 
analysis  the  effect  of  initial  velocity  on  warhead  radius  will 
be  studied.  The  effect  of  varying  the  initial  velocity  in 
equation  (16)  is  plotted  in  Figure  (4-14) .  The  results  were 
determined  for  various  explosives  to  show  their  effect  on  the 
warhead  radius.  If  an  initial  velocity  of  7000  ft/sec  is  chosen 
with  Composition  B  as  the  explosive  the  required  warhead  radius 
is  2.6  inches. 

From  this  analysis  the  warhead  weight  can  be  found.  The 
casing  thickness  is  given  by. 


t 

SL 

w 


.115  r 

w 


5  r 


w 


(from  required  C/M) 


Therefore,  the  explosive  weight  is 


W  =  it  ( r  -  t)  l  p 
e  w  w  e 


=  5.9530  kg 
e 


The  casing  weight  becomes, 


2  2 

W  =  ( rrr  -  trr  )  i  p 
c  w  e  w  Kc 


W  =  6.8651  kg 
c 
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The  total  warhead  weight  becomes , 

WWH  =  12.8181  kg  =  28.26  lb 

The  required  diameter  for  the  radar  antenna  was  10  inches 
and  exceeds  that  required  for  the  warhead.  If  the  warhead  is 
made  hollow  and  kept  at  the  same  weight,  equation  (14)  can  be 
reformulated  to  give 


k  -  fc»2  - 

21 
ri  , 

pe 

[2v  - t2 . 

o 

Q. 

The  hollow  portion  of  the  warhead  is  then  found  from 


) 


lw 


Figure  (4-15)  is  the  resulting  warhead. 


Figure  (4-15) .  Hollow  warhead. 
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V.  BASELINE  DEFINITION 


A.  CONTROL  CONCEPTS 

The  lifting  and  control  surfaces  of  a  missile  may  be  of 
monowing,  triwing  or  cruciform  configuration.  Figure  (5-1)  shows 
these  three  arrangements . 


Monowing  Triwing  Cruciform 

Figure  (5-1).  Control  configurations. 

The  monowing  arrangement  is  typical  of  most  cruise  missiles, 
which  require  long  range  and  low  drag.  For  this  type  of  arrange¬ 
ment  the  missile  must  bank  to  orient  the  lift  vector  for  a  maneuver. 
Because  of  this,  the  monowing  missile  is  not  as  rapid  in  maneu¬ 
vering  as  the  cruciform  configuration,  which  can  produce  lift 
in  any  direction  instantaneously.  The  cruciform  control  also 
has  identical  pitch  and  yaw  characteristics  which  results  in 
a  simpler  control  system.  The  triwing  configuration  is  used 
very  seldom  for  conventional  missiles.  It  can  be  shown  that 
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the  triwing  missile  requires  larger  wing  size;  therefore,  there 
is  little  drag  savings  even  though  there  is  one  less  wing 
than  with  the  cruciform  missile  [10].  From  this  discussion  it 
can  be  seen  why  the  cruciform  configuration  is  the  most  com¬ 
monly  used  for  tactical  missiles  of  short  or  medium  range. 

The  positions  of  the  lifting  surfaces,  on  the  missile  body 
depends  on  the  method  of  control  used  for  the  missile.  There 
are  three  conventional  methods  of  control  for  tactical  missiles. 
These  are;  1)  Wing  control;  2)  Tail  control;  3)  Canard 
control. 

1 .  Wing  Control 

Wing  control  missiles  normally  have  large  movable  wings 
located  slightly  behind  the  missile  center  of  gravity.  Because 
of  the  small  moment  arm  the  wing  surface  must  be  relatively 
large  to  provide  control  effectiveness.  As  would  be  expected, 
larger  activators  are  required  for  moving  the  wings.  A  positive 
deflection  of  the  wings  causes  a  positive  normal  force;  there¬ 
fore,  the  missile  reacts  almost  instantaneously;  thus  making 
wing  control  the  fastest  reacting  method  of  control.  Because  of 
the  smaller  moment  arm,  the  resulting  smaller  pitching  moment 
and  instantaneous  lift  result  in  smaller  angles  of  attack.  This 
feature  makes  wing  control  missiles  attractive  for  applications 
where  theincidence  angle  must  be  kept  small.  Air  breathing 
applications  often  use  wing  control  because  of  inlet  perfor¬ 
mance  degradation  at  higher  angles  of  attack.  This  type  of 
control  is  also  good  for  fixed  seekers. 


4*1 
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2.  Canard  Control 


Canard  control  is  normally  the  method  of  control  where 
the  movable  surfaces  are  placed  as  far  forward  as  possible. 

Because  of  the  resulting  larger  moment  arm,  smaller  surfaces  are 
required  to  provide  control  effectiveness.  Lift  on  the  missile 
is  still  developed  primarily  by  the  aft  (wing)  lifting  surface. 
Response  is  slower  than  wing  control  because  of  the  need  to 
pitch  the  missile  to  an  angle  of  attack  before  lift  is  developed. 
Higher  angles  of  attack  are  needed  to  generate  the  required  lift. 

One  advantage  of  this  type  of  control  is  convenience  of  packaging. 
Since  the  controls  and  avionics  are  forward  of  the  propulsion 
system,  the  need  for  connectors  is  elminated.  For  stability 
reasons  the  wing  of  a  canard  control  missile  must  be  located 
farther  aft  than  a  conventional  wing- tail  configuration.  The 
zero  lift  drag  is  normally  lower  than  wing  control  missiles 
due  to  smaller  surfaces. 

3.  Tail  Control 

Like  the  case  of  canard  control,  a  tail  control  missile 
has  the  movable  surface  as  far  from  the  center  of  gravity  as 
possible.  This  also  results  in  a  larger  moment  arm  and  there¬ 
fore  smaller  surface  required.  Control  deflection  is  the  oppo¬ 
site  of  that  for  canard  or  wing  control  since  a  negative  control 
deflection  results  in  a  pitching  moment  that  pitches  the  nose 
up  and  therefore  a  positive  lift  on  the  main  lifting  surface  (wing) . 
Tail  control  is  normally  the  slowest  method  of  control.  One 
advantage  is  that  the  flight  controls  are  at  the  end  of  the 
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missile  requiring  the  propulsion  system  to  be  located  forward 
of  other  types  of  control.  This  results  in  less  center  of  gravity 
movement  as  the  propellant  grain  burns.  This  becomes  a  definite 
advantage  for  longer  range  missiles. 

B.  GROSS  WEIGHT  AND  CENTER  OF  GRAVITY 

Since  the  component  weights  and  their  precise  locations  can¬ 
not  be  determined  at  this  point  in  the  design  process,  some 
method  of  estimating  the  total  weight  and  center  of  gravity  is 
needed.  There  are  two  approaches  commonly  used  to  find  the  gross 
weight  of  missiles  in  the  conceptual  design  phase.  One  is  through 
the  use  of  the  historical  data  discussed  in  Chapter  2.  Since 
the  warhead  weight  is  known,  an  estimate  of  the  gross  weight 
is  now. 


Another  method,  which  is  used  extensively  in  the  design  of 
aircraft,  is  the  use  of  regression  formulas  to  find  gross  weight 
or  component  weights  in  terms  of  parameters  that  are  known  early 
in  the  design.  Reference  (7)  has  derived  such  a  formula  for  the 
gross  weight  of  a  missile.  It  is, 

wQ  =  kg(l)2’13  (d)1,14 

Kg  =  Constant  to  be  determined  ^ 

L  =  Total  missile  length  (inches) 

D  =  Missile  diameter  (inches) 
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The  constant  K  in  equation  (1)  is  derived  from  a  baseline 
(generic)  missile.  In  this  case 


(WG) 


Ks  - 


baseline 


(L) 


2.13 

baseline 


(D) 


1.14 

baseline 


The  accuracy  of  a  regression  formula  such  as  equation  (1)  depends 
upon  how  close  the  synthesized  missile  is  to  the  baseline  for 
which  Kg  was  determined.  If  the  parameters  used  vary  more  than 
20  percent,  the  accuracy  of  the  regression  equation  decreases 
rapidly;  therefore,  if  the  parameters  length  and  diameter  vary 
significantly  from  those  of  the  baseline,  equation  (1)  may  give 
an  erroneous  estimate  of  the  gross  weight.  As  an  aid  in  deter¬ 
mining  the  gross  weight,  the  following  values  of  K_  were  derived 
from  data  given  in  reference  [7]. 


Missile 

SRAAM 

MRAAM 

LRAAM 

SAM 


Figure 

(5-2) 

(5-3) 

(5-4) 

(5-5) 


K, 


.00128 

.00118 

.00093 

.00108 


Figures  (5-2) ,  (5-3)  ,  (5-4) ,  and  (5-5)  show  the  generic  missiles 
from  which  these  values  were  derived. 

The  center  of  gravity  of  the  baseline  missile  for  this  section 
is  taken  to  be  60  percent  of  the  total  length.  At  this  point 
sufficient  information  has  been  developed  to  define  the  baseline 
missile  from  which  design  iterations  will  be  made.  The  lifting 
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Lgure  (5-2).  Generic  SRAAM  configuration  [7]. 


Figure  (5-3) .  Generic  MRAAM  configuration  [7] . 


Figure  (5-5) .  Generic  SAM  configuration  [7] . 


surface  planform  for  the  baseline  is  taken  as  a  delta  planform. 

As  will  be  shown  later,  this  wing  planform  will  be  very  close  to 
an  optimum  wing. 

C.  DESIGN  EXAMPLE  (BASELINE  DEFINITION) 

The  threat  for  which  this  missile  is  designed  is  a  highly 
maneuverable  fighter.  For  this  threat  a  canard  control,  cruci¬ 
form  configuration  is  chosen.  The  rationale  is  that  this  configur¬ 
ation  will  provide  the  fast  response  'necessary  at  minimum  drag. 

The  diameter  was  fixed  at  10  inches  due  to  antenna  considerations 
in  Chapter  3.  From  historical  data  in  Chapter  2, 

L  =  (£)  D  =  (15.89)  10  ■  158.9  inches 
L 

Ln  =  (— gO  D  =  22.3  inches 

The  MRAAM  of  Figure  (5-3)  is  used  as  a  generic  missile  for 
selecting  K_.  Inserting  length  and  diameter  into  equation  (1), 

Cj 

the  first  estimate  of  gross  weight  becomes, 

W_  =  .00118  (158. 9)2*13  (10)1,14  =  794.83  lb 

The  total  lifting  surface  required  is  then,  from  historical 
average , 

3  '  ‘^AVG  <V  ’  <^H794-83) 

S  =  9.02  ft2 
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Since  the  canard  to  wing  area  ratio  is  known  from  historical 


data, 

VSW  =0.20 

S„  =7.52  ft2 
w 

S  =  S  -  Sr  =  1.50  ft2 
c  w 

Lifting  surfaces  (delta  planform) .  From  historical  data, 

ARw  =  1.61 

AR  =3.74 
c 


AR 


s 


Figure  (5-6).  Lifting  surface. 


Wing:  (b)w  =  =  3.48  ft 

(C  )  =  2  S  /b  =  4.32  ft 

W 


Canard: 

(b)  =  /AR  i  =  2.37  ft. 

c  c  c 

(Cr)  =  1.27  ft. 

c 
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The  baseline  missile  is  now  defined.  The  canards  are  placed 
as  far  forward  as  possible.  The  wings  are  placed  as  far  aft 
as  possible  to  ensure  the  center  of  pressure  is  behind  the 
center  of  gravity.  The  exact  location  of  the  wings  will  be 
modified  in  the  next  chapter.  Figure  (5-7)  is  a  drawing  of 
the  baseline  missile. 
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VI.  LINEAR  AERODYNAMICS 


The  total  drag  on  the  missile  is  used  to  size  the  propul¬ 
sion  needed  for  the  cruise  segment  of  flight.  The  wing  size 
depends  upon  the  maximum  lift  required  by  the  missile.  The 
wing  and  tail  are  normally  placed  to  provide  minimum  drag  during 
cruise  or  a  certain  stability  margin  at  launch.  To  make  any 
of  the  above  calculations,  values  of  the  aerodynamic  coefficients 
are  needed.  This  chapter  presents  the  background  aerodynamic 
theory  necessary  for  these  initial  calculations.  The  theory 
used  is  linear  aerodynamic  theory  and  slender  body  theory,  from 
which  simple  calculations  can  be  made.  Where  linear  theory  did 
not  apply,  an  attempt  was  made  to  find  existing  empirical  expre- 
sions,  which  yield  results  accurate  enough  for  initial  calcula¬ 
tions.  The  full  nonlinear  theory  will  be  presented  in  Chapter  9. 
The  reference  area  for  all  coefficients  in  this  report  is  the 
missile  maximum  cross  sectional  area.  The  reference  length  is 
the  maximum  missile  diameter. 


A.  MISSILE  DRAG 

The  total  drag  of  a  missile  consists  of  zero  lift  drag,  Cj 
and  induced  drag,  Cn 


C0  -  CD0  *  S, 


(1) 


The  zero  lift  drag  can  be  found  from  a  component  build  up  method 
in  which  the  contributions  due  to  the  nose,  body  and  lifting 
surfaces  are  added  together  to  obtain  the  total  zero  lift  drag. 
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Care  must  be  taken  to  reference  the  appropriate  areas  when 
using  this  method.  The  total  is  then  multiplied  by  1.25  to 
account  for  interference  effects  and  variations  in  skin  roughness 
[12]  . 


1.25 


;v„  +  s 


)  +  (Cn  )  +  (c  ) 


B 


W 


(2) 


The  method  used  to  find  the  components  drag  depends  upon  the 

speed  regime  in  which  the  missile  is  operating.  Since  this 

report  is  concerned  primarily  with  supersonic  tactical  missiles, 

supersonic  zero  lift  drag  will  be  discussed  here.  The  component 

supersonic  zero  lift  drag  can  be  divided  into  skin  friction, 

C  ,  and  wave  drag,  C  . 

Uf  DW 


W 


V  +  <v 

1  w  w  w 


1.  Supersonic  Skin  Friction 

The  flow  over  a  body  traveling  at  supersonic  speeds  is 
likely  to  be  turbulent;  so  the  incompressible  skin  friction  co¬ 
efficient  is  given  by. 


C 


f . 
x 


.455 

(logiQ  Re> 


2.58 


(3) 


The  Reynolds  number  in  equation  (3)  is  based  upon  the  cruise 
altitude  and  speed  and  upon  the  characteristic  length  for  the 


component  being  determined.  The  Reynolds  number  is  given  by . 


where  x  =  L^,  the  length  of  the  nose 

L^b,  the  length  of  the  body  without  the  nose 

(c) ,  the  mean  aerodynamic  chord 


The  compressibility  correction  to  equation  (3)  is, 

-0.58 

Cf  =  Cfi  (1  +  °-15  MM  > 


(4) 


From  equation  (4)  the  skin  friction  drag  coefficient  for  each 
component  can  be  found  when  referenced  to  the  appropriate  area. 


(c  )  =  (cj  ( wetted^ 

f  N  «  ref 


^Swetted^ 


(CDf}  "  ^f*  S 

f  B  B  ref 


B 


The  lifting  surface  skin  friction  drag  is  determined  in  a  similar 
manner.  Care  must  be  taken  to  include  all  surfaces  in  the 
wetted  area  calculation  for  the  lifting  surfaces. 

2 .  Supersonic  Wave  Drag 

The  supersonic  wave  drag  consists  of  components  contri¬ 
buted  by  the  nose  and  lifting  surface.  Nose  wave  drag  depends 
on  the  shape  of  the  nose,  and  the  most  common  nose  shapes  are 


% 
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conical,  ogival  and  hemispherical.  Reference  [10]  lists 
empirical  formulas  for  finding  the  form (wave)  drag  of  various 
nose  shapes  at  zero  angle  of  attack. 

(1)  Conical;  C  =  (0.083  +  0 . 097/M  2) (a/10) 1* 69 
DW 

a  =  tan  =  Nose  semi-vertex  angle  (Degs) 

2LN 

The  center  of  pressure  for  a  conical  nose  is  at  the  centroid  of 
the  nose  planform  or  two  thirds  the  length  of  the  nose. 

2[196(Ln/D)2  -  16]  ) 

28  (M  +  18)  (Ln/D)2  ) 


(2)  Ogival;  C 


=  P<  1  - 


w 


P  =  (C  )  for  conical  nose 
DW 

The  center  of  pressure  for  an  ogive  noise  is, 

fp  _  1  [50 (M  +  18)  +  7M2P  (5M  -  18) 
LN  2  40 (M  +  18)  +  7M2P  (4M  -  3) 


The  semi-vertix  angle  for  an  ogive  is  twice  the  equivalent  cone 
angle. 

(3)  Hemispherical;  The  drag  on  a  hemispherical  nose  is  extremely 
high  compared  to  other  nose  shapes,  and  is  difficult  to 
estimate.  An  initial  estimate  of  the  wave  drag  can  be 
found  from  Figure  (6-1) . 
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Figure  (6-1) .  Hemispherical  nose  wave  drag 


The  wave  drag  due  to  the  lifting  surfaces  can  be  found 
using  the  methods  of  reference  [12],  For  a  double  wedge  air¬ 
foil  with  sharp  leading  edges  as  shown  in  Figure  (6-2) ,  the 
wave  drag  is  given  by  the  following  formulas: 


(1)  Supersonic  leading  edge. 


(2)  Subsonic  leading  edge. 


B  cot  ale  (l) 


'w 


ref 


Figure  (6-2) .  Double  wedge  wing  [12] . 


where  n  _  c/xt 

B  "  1  -  xfc/c 

n  r 


A  =  leading  edge  sweep 

LiE 


S„  =  planform  area 


3.  Base  Pressure  Drag 

The  drag  contribution  of  the  blunt  base  for  non-boattailed 
bodies  is  given  in  reference  [4]  as. 


This  term  assumes  no  jet  thrust  from  the  base  of  the  missile, 
or  that  the  missile  is  operating  in  the  power  off  condition. 
This  term  is  not  included  in  equation  (2) ,  which  is  not  a  bad 
assumption  for  powered  flight  where  the  nozzle  exit  area  is 
approximately  equal  to  the  base  area  of  the  missile.  For  the 
case  where  the  nozzle  exit  area  is  much  less  than  the  base 
area  as  in  Figure  (3),  the  base  pressure  contribution  should  be 
included. 


Figure  (6-3).  Base  pressure  areas  [10]. 
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In  this  case  the  base  drag  is. 


=  shaded  area  of  Figure  (6-3) 


4 .  Induced  Drag 

The  induced  drag  on  a  missile  is  the  drag  due  to  lift. 
This  drag  is  caused  by  the  component  of  the  lift  vector  in  the 
drag  direction.  For  supersonic  flow  the  induced  drag  is  given 
by. 


C 


D 


i 


2 


where;  w 

The  lift  curve  slope,  CNa,  will  be  developed  later  in  this  chapter. 

B.  DESIGN  EXAMPLE  (ZERO-LIFT  DRAG  CALCULATION) 

The  thickness  to  chord  ratio  of  the  wing  and  tail  have  not 
yet  been  determined;  however,  it  is  desirable  to  construct  the 
lifting  surface  as  thin  as  structurally  possible  to  minimize  the 
wave  drag.  Since  structures  have  not  been  covered  the  minimum 
thickness  to  chord  ratio  is  estimated  at  3  percent.  The  flight 
and geometric  conditions  determined  thus  far  are; 
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h  =  10,000  ft 


“m  58  3-° 


s 


ref 


(t/c)w  = 


0.5454 
( t/c ) c*  0 . 0  3 


Sw  =  7.52  ft 
Sc  =  1.50  ft 


L  =  1.8583 
N 

Lcb  =  11.3833 
(c)w  =  2.8816 
(c)c  =  0.8467 


From  these  conditions  the  following  table  can  be  constructed. 


Nose 

Afterbody 

Wings 

Canards 

X 

1.8583 

11.3833 

2.8816 

0.8467 

wetted 

3.3394 

29.8014 

15.0400 

3.0000 

R 

e 

2.984 

18.282 

4.6280 

1.360 

Cfi 

0.0025 

0.0020 

0.0024 

0.0029 

C£ 

0.0015 

0.0012 

0.0015 

0.0017 

C'>£ 

0.0095 

0.0652 

0.0827 

0.0187 

C^w 

0.1545 

— — 

0.0400 

0.0050 

CD_ 

0.1640 

0.0652 

0.1227 

0.0237 

The  zero  lift  drag  for  the  wings  and  canards  in  the  above  table 
take  into  account  that  there  are  two  sets  of  wings  (4  panels) . 
The  total  missile  zero  lift  drag  then  becomes, 


1.25 


[(V  + 

u  N 


(c_  )  +  (c„  )  +  (c 


AB 


W 


D  >  1 

0  CJ 


C  =  0.4695 
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C.  MISSILE  LIFT  CURVE  SLOPE 


T(B){ 


>L 
'  B(W) 


Jw(B) 


Figure  (6-4) .  Wing-body- tail  lift  [13] . 


For  the  purposes  of  this  section  lift  on  a  wing-body-tail 
combination  can  be  taken  as  the  sum  of  the  components  in  Figure 
(6-4) .  These  consist  of, 

1^  =  lift  on  the  nose 

L  =  lift  on  the  wing  in  the  presence  of  the  body 
W(B) 

Lfi  =  additional  lift  on  the  body  due  to  the  presence 
a (W)  of  the  wing 

L  -  lift  on  the  tail  in  the  presence  of  the  body 
T(B) 

L_  =  additional  lift  on  the  body  due  to  the  presence 
b(T)  of  the  tail 

The  lift  of  only  the  wing-body  combination  can  be  defined  as,  L^, 


where. 


L  56 


The  lift  of  the  wing  alone,  1^,  is  that  obtained  from  thin 
airfoil  theory  or  experiment  considering  the  exposed  area  only. 
The  constant  in  equation  (5)  is  defined  as 


+ 


KW 


(B) 


From  this  equation,  K  is  the  ratio  of  the  additional 

B(W) 

body  lift  in  the  presence  of  the  wing  to  the  wing  alone  lift  for 
zero  control  deflection,  6=0. 


K 


(W) 


(CL  ) 
La 
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6  =  0 


6=  0 


W 


The  interference  factors,  K  and  K..  have  been  deter- 

B(W)  ™(B) 

mined  from  slender  body  theory  for  wing-body  combinations. 

Figure  (6-5)  is  a  plot  of  these  values.  The  wing  alone  lift 

curve  slope,  (C_  )  is  determined  from  thin  airfoil  theory  or 
L“  W 

experiment.  The  lift  curve  slope  of  delta  wings  with  super¬ 
sonic  leading  edges  is  given  by, 

CLa  =  |  where  S  =  V  MM2_1  k  <  £ 

For  subsonic  leading  edges  this  becomes. 


„  _  27Ttane/tany 

La  =  E8 


k  >  e 
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Where  E  is  the  elliptic  integral  of  the  second  kind  for 


vl  -  (tane/tany)  , 

Figure  (6)  shows  the  Mach  angle  y  and  sweep  angle,  e. 


Figure  (6-6) .  Wing  leading  edge. 


The  lift  of  the  nose  is  that  obtained  from  slender  body 
theory.  For  small  angles  of  attack 


“W 


N 


2/rad. 


The  tail-body  combination  lift  is  determined  in  the  same 
manner  as  the  wing.  The  above  equations  are  defined  for  the  case 
where  the  incidence  angle  is  zero  and  the  angle  of  attack,  (a) , 


160 


(CW  =  “V  +  (CNa)  S^"  +  (CNa) 

NaCM  NCt  N  W(B)  SrSf 


W 


B(W)  ref 


+  (CNa}  S  J  +  (CNcJ  sr:f 

T(b)  ref  B(t)  ref 


From  the  definition  of  the  interference  factors, 


W 


(C^J  =  <CmJ  +  {Kw  +  K  )  F- 

N^CM  N  N  N  5W  ”(B)  B(W)  Sref 


(CN  )  (K  +  K  ) 

Na  T  t(b)  b(t)  Sref 


(6) 


If  the  tail  is  the  control  surface,  a  similar  development  for 
the  control  effectiveness,  CN ^  yields  * 


=  ^CNa^  T  ^B  ^  S 

NO  CM  Na  T  T(b)  B(t)  Sref 


D.  MISSILE  PITCHING  MOMENT 

From  Figure  (6-7)  the  moment  about  the  center  of  gravity 
can  be  found.  The  centroid  of  the  wing  is  the  location  of  the 
center  of  pressure  of  a  wing  alone  in  supersonic  flow.  The 
effect  of  the  wing-body  combination  is  to  move  the  center  of 
pressure  aft  and  as  can  be  seen  in  Figure  (6-7)  the  center  of 
pressure  for  the  additional  lift  on  the  body  due  to  the  presence 
of  the  wing  is  aft  of  this  location.  Reference  [13]  has  an 
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■*}  ^(B) 
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f  XB(W) 

- 


T(B) 


*  B  (T) 

Figure  (6-7) .  Forces  Acting  on  the  Missile. 

excellent  discussion  on  how  to  find  these  centers  of  pressure. 

For  the  purpose  of  this  chapter,  which  is  initial  sizing  and 

placement  of  the  lifting  surfaces,  the  center  of  pressure  for 

both  of  these  forces,  N  and  N  is  taken  as  the  centroid 

(B)  B(W) 

of  the  wing  planform.  With  this  assumption,  the  moment  about 
the  center  of  gravity  is, 

(m)cm  =  nn  xn  +  (nw  +  nb  ]  ^  +  (nt  +  nb  }  xt 

CM  N  N  W(B)  B(w)  W  T(b)  B  (t)  T 

Following  the  same  development  as  for  CN 

CMa  •  “W„  3  +  <CNo.>w  <KW(B)  +  KB(W)>  (-5>  (s^7» 


+  (CNa)  (KT  +  KB  }  (“§)  (S^~) 
Na  T  T(B)  B(T)  D  Sref 


(7) 
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Also  if  the  tail  is  the  control  surface. 


'M<5 


(CNa\ 


S„ 

+  k  )  (~g)  (3— 

(B)  b(T)  d  sref 


Ome  must  be  careful  in  defining  the  moment  arms  in  the 
above  equations .  If  a  nose  up  pitching  moment  is  developed  in 
Figure  (6-7) ,  the  moment  arm  is  positive.  Conversely,  a  nega¬ 
tive  moment  arm  means  a  nose  down  pitching  moment  is  developed. 

E.  WEIGHT  AND  CENTER  OF  GRAVITY  VARIATIONS 

From  the  preceding  sections  of  this  chapter,  it  can  be  seen 
that  the  analysis  is  normally  for  one  point  in  the  flight  pro¬ 
file.  Since  the  launch  condition  has  been  defined,  the  tail 
sizing  can  be  accomplished  for  this  initial  condition.  The 
variation  of  missile  weight  and  center  of  gravity  is  due  to 
propellant  burning.  For  a  solid  propellant  this  variation  can 
be  quite  large.  Since  no  information  is  available  on  the  pro¬ 
pellant  at  this  point  in  the  design,  the  following  guidelines 
will  serve  for  initial  calculations.  For  air-to-air  missiles 
the  propellant  is  approximately  35  percent  of  the  launch  weight. 
The  center  of  gravity  travel  is  5  percent  of  the  length.  For 
surface-to-air  missiles  which  must  be  boosted  to  flight  speed, 
the  propellant  weight  can  be  taken  to  be  48  percent  of  the 
launch  weight  with  a  corresponding  center  of  gravity  travel  of 
approximately  8  percent  of  the  length.  These  are  initial  approxi¬ 
mations  taken  from  historical  data  and  can  be  refined  later  in 
the  design  process. 


VII.  LIFTING  SURFACE  DESIGN 


From  the  analysis  of  wing  lift  and  drag  in  Chapter  6,  it  can 
be  seen  that  the  performance  of  a  wing  will  vary  with  planform. 

Up  to  this  point  the  lifting  surfaces  have  been  considered  delta 
planforms  with  zero  taper  ratio.  This  is  not  necessarily  the 
optimum  planform  since  it  was  taken  from  a  historical  average 
and  does  not  apply  to  a  specific  missile.  This  chapter  deals  with 
the  sizing,  placement  and  planform  definition  of  the  missile. 

A.  WING  PLACEMENT 

The  wing  placement  on  the  missile  depends  upon  the  type  of 
control  used.  For  canard  control  the  wing  (aft  surface)  is  nor¬ 
mally  fixed  as  far  aft  as  possible  for  stability  purposes.  For 
a  wing  control  or  tail  control  missile,  where  the  wing  is  near 
the  center  of  gravity,  the  wing  placement  becomes  more  critical 
and  depends  upon  the  stability  margin  required.  Since  the  drag 
during  the  cruise  segment  includes  the  drag  due  to  lift,  one 
method  to  minimize  drag  would  be  to  place  the  wing  such  that 
zero  lift  is  produced  on  the  tail  during  cruise.  The  moment  about 
the  center  of  gravity  is  zero  for  trimmed  flight,  thus  equation 

(7)  yields,  for  (CM)  =  0. 

T 
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B.  MANEUVER  LOAD  FACTOR 

Regardless  of  the  type  of  missile  being  designed,  it  will  be 
required  to  maneuver  in  order  to  intercept  its  intended  target. 
For  an  air-to-air  or  surface-to-air  missile  the  level  of  this 
maneuver  may  be  quite  large.  The  maneuver  load  factor  required 
was  found  in  Chapter  3.  The  maximum  maneuver  the  missile  can 
sustain  depends  upon  the  maximum  trimmed  normal  force  the  missile 
can  develop. 


L  =  nW 


Figure  (7-1) .  Sustained  maneuver  of  a  missile. 


From  Figure  (7-1)  the  force  developed  by  a  missile  in  a  con¬ 
stant  acceleration  turn  is 


L  =  nW 

n  =  maneuver  load  factor 


(1) 
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»  normal  force  are  approximately  equal.  This  approximation  is 


good  for  angles  of  attack  of  up  to  10  degrees,  and  above  this 
value  of  maximum  trimmed  angle  of  attack  the  linear  theory 
becomes  inaccurate.  With  this  assumption  for  now,  equation  (1) 
becomes , 


L  =  N  =  nW 

(C„  )a  q  S  =  nW 
Na  ref 


(2) 


«W 


nW 


required 


^max  q  Sref 


Equation  (2)  gives  the  lift-curve  slope  required  to  develop 

the  required  normal  force  at  a  trimmed  angle  of  attack,  (a) 

^  max 

The  lift  curve  slope  developed  is  given  by  equation  (6)  in  Chapter 

6. 


'w 


(3) 


(CW  =  {Cnu)  +  {CNa!  (KW  +  KB  ]  S 
Na  CM  N  N  N  W  W(B)  B(W)  Bref 


+  (CNa5  (KT  +  KB  } 

wa  T  t(b)  b(t}  sref 


As  discussed  earlier  the  values  of  the  interference  factors 

depend  upon  the  wing  planform.  Typical  values  of  (K^  +  K  ) 

W(B)  B(W) 

and  (K  +  K  )  are  from  1.5  to  2.0.  As  a  conservative 
T<B)  B(T) 

estimate  a  value  of  1.5  is  assumed.  Equation  (3)  then  becomes, 


(CM  )  =  (C  )  +1.5  (C  ) 

Na  CM  Na  N  Na  W  Sref 
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+  1.5  (C„  ) 


Na  T  Sref 


(4) 


. — . . -  - 


Selecting  supersonic  leading  edges. 


"Ww  -  <«wT  - 


Therefore 
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CM 


lCNa> 
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8  S 


ref 


(Sw  ♦  sT) 


If  the  complete  missile  lift-curve  slope  is  set  equal  to  the 
required  lift-curve  slope,  the  lifting  surface  area  required  to 
maintain  the  maneuver  load  factor  can  be  found,  using  equation  (4) . 
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req 


=  [(CNa> 


req 


~  «W 
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(SW  +  V 


K  = 


req_ 


baseline 


If  the  same  ratio  of  wing  to  control  surface  area  is  kept  to 
minimize  stability  perturbations,  the  new  wing  and  tail  area  are. 
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K(V 

baseline 

K(S  ) 

baseline 


The  above  analysis  assumes  a  linear  variation  of  CNa  which  is 
good  for  small  angles  of  attack.  A  more  precise  analysis  will 
be  performed  in  Chapter  9.  The  above  analysis  should  be  performed 
at  the  expected  encounter  conditions. 
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1 .  Design  Example  (Maneuver  Load  Factor) 

Since  a  canard  control  was  chosen  for  the  design  example 
the  canard  and  wing  position  are  fixed  as  fqr  forward  and  aft 
as  possible.  The  analysis  is  done  for  a  conservative  Mach  number 
of  2.5  after  the  missile  has  slowed  from  its  cruise  velocity  due 
to  the  maneuver.  From  previous  results, 

n  =  31  g*s 

(S  +  S  )  *  9.02  ft2 

baseline 

W_,  =  794.83  lbs 

VJ 

If  the  missile  is  required  to  maneuver  at  one-half  its  powered 
range,  and  the  propellant  weight  is  35  percent  of  the  launch 
weight,  the  maneuver  weight  becomes 

W  *  655.73  lbs 

The  dynamic  pressure  for  MM  *  2.5  and  at  10,000  feet  altitude  is 


q  =  6369.84  lbs  ft2 

The  required  lift-curve  slope  becomes  for  (a  _„)  =  10  degrees 

IuclX 

=  .1745  rad. 


(C„  )  =  33.5311  per  rad. 

req 

The  required  lifting  surface  area  is; 


(SW  +  V 


[33.5311-2] 


6  S 


=6.57  ft2 


req 


168 


Therefore 


K  =  0.7284  . 


The  required  wing  and  canard  sizes  to  achieve  the  maneuver  are; 

Sw  =  (7.52)  (.7284)  =  5.48  ft2 
Sc  =  (1.50)  (.7284)  =  1.09  ft2 


C.  TAIL  SIZING 

For  a  wing-tail  combination,  the  primary  concern  in  the 
sizing  of  the  tail  is  the  static  stability  of  the  missile.  The 
missile  becomes  more  stable  as  the  mission  proceeds,  since  the 
center  of  gravity  moves  forward  as  the  propellant  burns.  As 
the  missile  becomes  more  stable,  control  of  the  missile  becomes 
sensitive.  If  the  tail  is  sized  for  zero  static  stability  at 
launch,  the  missile  control  will  remain  more  effective  during 
flight.  This  is  the  best  condition  possible  without  the  use  of 
some  form  of  stability  augmentation  system  at  launch.  Therefore 

at  launch  C„  =0. 

Na 


0  *  (CM  )  +  (C„  )  X..  (K.  +  K_  ) 

Na  N  "  Na  W  "  Sref  "(B)  B(W) 


(5) 


'Nc'  “T  S—  (KT  +  KB  ) 
Na  T  T  ref  T(B)  B(T) 


With  the  lifting  surface  area  fixed  due  to  the  maneuver  load 
factor,  the  tail  can  be  sized  to  satisfy  equation  (5) .  As  can 
be  seen  the  position  depends  highly  on  the  moment  arms  that  the 
lift  forces  act  through.  Since  the  missile  length  has  not  been 
fixed  at  this  point  in  the  design  and  may  vary  due  to  propulsion 
requirements,  this  analysis  will  be  completed  later. 


D.  WING  PLANFORM 


The  wing  planform  is  specified  by  the  leading  edge  sweep, 
Ale,  taper  ratio,  X,  aspect  ratio,  AR  and  planform  area  Sw-  The 
planform  area  was  fixed  due  to  maneuvering  requirements  in  a 
previous  section.  This  section  is  concerned  with  defining  the 
remaining  planform  parameters.  Figure  (7-2)  is  the  wing  plan- 
form  and  the  equations  used  to  define  these  parameters. 


2 

AR  =  b  /s 


X  *  Vcr 
2  s 

Cr  b (1+A) 


(i±2£i i) 


~  2  tan  ALE 


Figure  (7-2) .  Planform  geometry. 


The  value  of  the  lift-curve  slope  used  previously  was  derived 
from  linear  theory  and  is  applicable  only  to  simple  planforms. 

For  the  purpose  of  this  section,  which  will  include  low  aspect 
ratio  wings  with  supersonic  and  subsonic  leading  edges.  Figures 
(7-3)  and  (7-4)  are  used.  These  figures  have  been  corrected  for 
3-D  effects. 

1.  Effect  of  Taper  Ratio  and  Leading  Edge  Sweep 

The  lift  and  drag  characteristics  of  the  wing  are  the 
primary  parameters  of  interest.  The  objective  in wing  design  is 
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•  nr 


(pernd) 


(  ^ a )  theory 


is  to  obtain  maximum  lift  with  minimum  drag.  It  will  be  shown 
that  these  objectives  are  usually  conflicting;  therefore,  some 
compromise,  or  optimum,  planform  must  be  found.  From  Chapter  6 
the  equations  for  the  drag  of  a  wing  in  supersonic  flow  are  found. 
The  skin  friction  drag  depends  upon  the  mean  geometric  chord. 
Equation  (3)  in  Chapter  6  indicates  a  larger  mean  geometric 
chord  would  result  in  less  skin  friction  drag  for  a  fixed  area. 

The  mean  geometric  chord  can  be  expressed  in  terms  of  taper 
ratio  as  follows. 


-  _  2  2  s  /l+  X  +A2  \ 

c  "  3  bTI+TT  \  I+X  / 


i  /Zi  i+x+x2 

3  V  AR  (1+A)2 


This  equation  leads  to  a  zero  taper  ratio  to  maximize  c  and 
reduce  skin  friction  if  the  surface  area  and  aspect  ratio  were 
fixed.  The  wave  drag  is  constant  for  supersonic  leading  edges, 
and  decreases  when  the  leading  edge  goes  subsonic  or  the  leading 
edge  is  behind  the  mach  line.  The  lift  capabilities  also  de¬ 
crease  as  the  leading  edge  goes  subsonic. 

As  stated  earlier  this  leads  to  conflicting  performance  since 
the  objective  is  to  minimize  drag  while  maximizing  lift.  At 
this  point  an  example  best  illustrates  the  results  of  varying 
taper  ratio  and  leading  edge  sweep.  The  lift  -  curve  slope  is 
derived  using  the  methods  of  reference  (1).  Figures  (7-3)  and 
(7-4)  are  from  reference  [12]  and  are  used  to  find  the  supersonic 
linear  lift  corrected  for  3-D  flow  effects.  The  drag  methods 
of  Chapter  5  are  used  to  determine  the  drag  characteristics. 
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2.  Example 

2 

For  a  fixed  wing  area  of  4  ft  and  an  aspect  ratio  of 
2  the  lift  curve  slope  and  drag  were  determined  for  a  flight 
Mach  number  of  M=  2.  This  Mach  number  corresponds  to  a  sonic 
leading  edge  sweep  of  60°.  The  wing  leading  edge  sweep  was 
varied  from  zero  to  75°  for  taper  ratios  of  0,  1/2  and  1.  The 
results  are  plotted  in  Figures  (7-5)  and  (7-6) . 

From  this  example  certain  generalizations  can  be  made. 
From  Figure  (7-5)  it  can  be  seen  that  the  drag  is  relatively 
insensitive  to  taper  ratio.  There  is  a  reduction  in  drag  for 
increased  leading  edge  sweep.  For  this  example  there  is  approxi¬ 
mately  a  6  percent  drag  reduction  for  every  five  degree  increase 
beyond  the  sonic  value.  Figure  (7-6)  indicates  that  the  effect 
of  decreasing  the  taper  ratio  is  to  delay  the  drop  in  the  lift 
curve  slope  of  the  wing.  From  this  example  a  general  guideline 
would  be  to  fix  the  leading  edge  sweep  5  degrees  beyond  the  sonic 
value  and  the  taper  ratio  at  zero.  This  would  provide  a  6  per¬ 
cent  reduction  of  wave  drag  while  maintaining  the  maximum  lifting 
capabilities  of  the  wing. 

3.  Effect  of  Varying  Aspect  Ratio 

The  result  of  increasing  aspect  ratio  is  an  increase  in 
the  lift-curve  slope  of  the  wing  [12]  .  The  aspect  ratio  is 
given  by, 

AR  *  b2/s 

For  a  zero  taper  ratio  wing  this  becomes; 

AR  =  2b/cr  (6) 
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Effect  of  A  on  drag 


1 


From  equation  (6)  it  can  be  seen  that  increasing  aspect  ratio 
results  in  a  corresponding  increase  in  wing  span  for  a  given 
root  chord.  Missiles  are  normally  span  limited,  due  to  launcher 
constraints;  therefore,  there  is  a  maximum  span  which  can  be 
accomodated.  The  drag  of  the  wing  again  is  in  conflict  with 
the  lift  since  increasing  aspect  ratio  decreases  the  mean  geo¬ 
metric  chord  for  a  given  span  and  therefore  increases  skin  fric¬ 
tion  drag.  A  compromise  wing  AR  must  be  found  which  considers 
both  the  lift  and  the  drag  characteristics.  For  a  discussion 
of  optimum  aspect  ratio,  the  following  functions  are  defined 
for  convenience: 


F  =  F^  +  F2  “  Lift-drag  function 


(7) 


where 


F1  H  7c — J - ^Normalized  drag  function 

D0  max 


1/CNa 

F _  =  =7~.' ■  — s  Normalized  lift  function 

A/  . 

min 


From  equation  (7)  if  F  is  plotted  over  the  allowable 
range  of  aspect  ratio,  a  minimum  value  of  the  lift-drag  function 
fixes  the  desired  aspect  ratio.  This  aspect  ratio  corresponds 
to  minimizing  the  drag  function,  F^,  while  maximizing  the  lift 
function,  F2 .  For  convenience  the  abscissa  is  plotted  as  c/c^maxj 
as  shown  in  Figure  (7-7) . 

From  Figure  (7-7)  the  optimum  mean  geometric  chord  cor¬ 
responds  to  point  A. 
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Figure  (7-7) .  Lift-drag  function. 


c 

max 


The  aspect  ratio  is  then  given  by 


b  =  2s/cr 
AR  =  b2/s 


As  mentioned  earlier  a  missile  is  normally  span  limited.  The 
plot  of  the  lift-drag  function  F  is  normally  fairly  flat  on  the 
left,  or  for  c/c^maxj  approaching  zero.  For  this  reason  point 
B  may  be  chosen  as  the  optimum  since  F  varies  very  little  up 


m 


to  this  region.  Point  B  corresponds  to  increased  chord  and 
decreased  span. 

4 .  Design  Example  (Wing  Planform) 

From  the  previous  analysis  the  following  parameters  were 

defined. 


«M  =  3‘° 

S„  =  5. 48  ft2 

If  the  leading  edge  sweep  is  fixed  5  degrees  behind  the  mach 
line. 


y  =  sin 

'l  r  *  19-5° 

ale  -  95 

-  u  =  75.5° 

A  =  0 

planform  table 

becomes 

AR 

b 

c 

r 

c 

KW(B) 

kb(w) 

2.0 

3.31 

3.31 

2.21 

1.17 

.26 

in 

• 

r— 1 

2.87 

3.82 

2.55 

1.19 

.30 

1.0 

2.34 

4.68 

3.12 

1.22 

.35 

.75 

2.03 

5.40 

3.60 

1.24 

.40 

o 

in 

• 

1.66 

6.62 

4.41 

1.29 

.48 

.25 

1.17 

9.36 

6.24 

1.36 

.62 

The  wave  drag  is  constant  and  is  given  by. 


4  cot  ale 


C  ■  .0094 
DW 
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AR 

% 

C°0 

AR  tan  Ale 

6 

tan  ale 

(C„  )  (K__  +K„  1 

Nct  W  **(B)  B(W) 

2.0 

.0303 

.0397 

7.73 

.73 

2.1137 

1.5 

.0296 

.0390 

5.80 

.73 

1.9845 

1.0 

.0288 

.0382 

3.87 

.53 

1.8271 

.75 

.0282 

.0376 

2.90 

.73 

1.7813 

.50 

.0274 

.0368 

1.93 

.73 

1.4809 

.25 

.0260 

.0354 

.97 

.73 

.7937 

AR 

c/c 
'  max 

F1 

F2 

F 

2.0 

.35 

1.0000 

.3755 

1.3755 

1.5 

.41 

.9824 

.3999 

1.3823 

1.0 

.50 

.9622 

.4344 

1.3966 

.75 

.58 

.9471 

.4456 

1.3927 

o 

in 

• 

.70 

.9270 

.5360 

1.4630 

.25 

1.00 

.8917 

1.0000 

1.8917 

From  Figure  (7-8)  it  can  be  seen  that  the  lift-drag 
function  remains  relatively  constant  up  to  c/cmax  =  .58.  Since 
the  object  is  to  make  the  span  as  small  as  possible  this  point 
is  taken.  The  wing  planform  becomes. 
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c 


=  .58  c  =  3.62  ft  -  43.44  in 
max 

=  5.43  ft  =  65.16  in 
b  =  2.02  ft  =  24.24  in 
c  =  3.90  ft  =  46.80  in 
AR  =  0.74 


5.  Design  Example  (For  Canard  Planform) 

From  previous  analysis  and  development,  the  following 
parameters  have  been  defined: 


M  =  3.0 

S  =  1.09  ft2 
c 

(t/c)c  =0.03 


If  the  leading  edge  sweep  is  set  5  degrees  behind  the  sonic 
value  and  the  taper  ratio  is  set  to  zero,  the  following  planform 
table  results. 


ALE  =  75'5 


X  =  0 


AR 

b 

c 

r 

2.5 

1.63 

1.34 

2.0 

1.48 

1.47 

1.5 

1.28 

1.70 

1.0 

1.04 

2.10 

.5 

.74 

2.94 

.25 

.52 

4.19 

c 

(B) 

K 

B(W) 

.89 

1.29 

.52 

.98 

1.31 

.55 

1.14 

1.35 

.59 

1.40 

1.47 

.76 

1.96 

1.51 

.90 

2.79 

1.59 

1.02 
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Figure  (7-8) .  Wing  lift-drag  function. 


AR 

% 

C°0 

AR  tan 

6 

tan  ale 

Csa(k) 

2.5 

.0069 

.0087 

9.67 

.  73 

2.0 

.0067 

.0085 

7.73 

.73 

2.7900 

1.5 

.0067 

.0085 

5.80 

.73 

3.6339 

1.0 

.0064 

.0082 

3.87 

.73 

2.5952 

.5 

.0061 

.0079 

1.93 

.73 

2.0880 

.25 

.0058 

.0076 

.97 

.73 

1.0124 

AR 

c/c 
'  max 

F1 

F2 

F 

2.0 

.35 

1.0000 

.3628 

1.3628 

1.5 

.41 

1.000 

.3844 

1.3844 

1.0 

o 

in 

• 

.9647 

.3901 

1.3557 

.5 

.70 

.9294 

.4849 

1.4143 

.25 

1.00 

.8941 

1.0000 

1.8941 

From  Figure  (7-9)  it  is  seen  that  there  is  a  minimum  of 
the  lift-drag  function  at  c/c  =  0.5.  The  optimum  canard  is 

ITiaX 

then. 


c  *  0.5  c  =  1.40  ft 
max 

c  *  2.09  ft  =  25.08  in 
r 

b  =  1.04  ft  =  12.48  in 
AR  =  0.99 

c  =  2.01  ft  =  24.12  in 


=  16.80  in 


Figure  (7-9) .  Canard  lift-drag  function. 


E.  DESIGN  EXAMPLE  (REVISION  OF  ZERO  LIFT-DRAG  AND  LIFT-CURVE 


SLOPE) 

Wing:  c  =  3.62  ft 

cf  =  .0014 

(C_  )  =  .0750 

u0  W 

Canard:  c  =  1.40  ft 

cf  =  .0016 

(C_  )  =  .0166 

u0  c 

Body:  (C  )  =  .2292 

U0  B 

The  complete  missile  CD  including  the  interference  factor 

of  1.25,  is  now  0.4010.  From  the  previous  calculation  of  the 
baseline  zero  lift-drag  coefficient,  the  drag  has  been  reduced 
by  14.6  percent. 

From  equation  (6)  in  Chapter  6  the  lift-curve  slope  is  now 
CNa  =  25.17/rad 


Figure  (7-10)  is  the  missile  design  to  this  point. 


VIII.  PROPULSION  REQUIREMENTS 

The  following  discussion  presents  a  method  for  preliminary 
sizing  of  a  solid  rocket  motor  for  a  boost-sustain  trajectory 
of  an  air  launched  or  a  surface  launched  tactical  missile.  The 
analysis  consists  of  sizing  the  booster  from  incremental  velo¬ 
city  considerations  and  sizing,  the  sustainer  for  the  maximum 
range  required  at  the  operational  altitude.  The  method  assumes 
a  constant  acceleration  boost  and  a  constant  altitude  cruise. 


A.  BOOSTER  INITIAL  SIZING 

Since  the  control  system  cannot  respond  properly  while  in 
the  boost  phase  of  flight,  it  becomes  important  to  make  the  boost 
time  as  short  as  possible.  The  limiting  factor  for  boost  time 
is  the  maximum  axial  acceleration  the  airframe  and  components 
can  withstand.  This  maximum  acceleration  is  normally  around 
30  gfe.  The  boost  time  then  becomes. 


For  constant  acceleration. 


VV1 


AV. 


B 


(1) 


The  incremental  velocity  during  boost,  AVfi,  can  be  derived 
from  the  equations  of  motion. 
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Figure  (8-1) .  Forces  acting  on  the  missile. 


If  during  the  boost  phase  the  thrust,  drag  and  launch  angle  are 
considered  constant,  the  axial  acceleration  of  the  missile  is 
constant  and  may  be  written  as: 

dv 

mjj7=T-D-W  sm  yL 

dV  =  g  (^~.)  dt  -  g  sin  yl  dt 


If  the  launch  velocity  is  and  the  velocity  at  the  end  of 
boost  is  V2, 


(2) 


The  vehicle  weight  in  equation  (2)  is  a  function  of  time.  If  the 
propellant  weight  is  given  by  W  ,  and  the  propellant  grain  burn 

r 

is  linear  with  time,  the  missile  weight  can  be  expressed  as, 


w  =  W.  -  W 

L  p  fcb 


where  W_  =  Launch  weight. 
Equation  (2)  then  becomes, 


avb  =  v2  -  vx 


g(T  -  D) 


fh  tbdt 

J  (wLtb-wpt)  sin  yl  tb 


Integrating  the  first  term  yields. 


avb  = 


g(T  -  D)tt 
_ 


Jin 


Vb 


-  9  sin  yL  tb 


Since  the  empty  weight  is  given  by 


=  WT  -  W 
E  L  .  p 


The  incremental  velocity  can  be  written  as 
g (T  -  D)t,  W. 

- _ o  0  _  ii 


avb  = 


W 


inwZ'g  sin  yl  fcb 

£j 


(3) 


A  first  estimate  of  propellant  weight  can  be  made  if  the  drag  is 
assumed  zero  in  equation  (3)  ,  and  Igp  =  Ttb//W 

P 

iVB  ■  *  rsp  ln  wT  -  9  sin  Yl  tb  (4) 

£ 


In  the  above  expression  the  specific  impulse,  Igp,  is  charac¬ 
teristic  of  the  propellant  chosen  and  can  be  found  from  histori¬ 
cal  data,  or  a  specific  propellant  value  may  be  used. 
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WL 

Equation  (4)  can  then  be  solved  for  —  and  the  first  esti- 

WE 

mate  of  propellant  weight,  W ,  can  be  found  as  follows. 


MP  -  V1  -  iT> 

Li 


With  the  propellant  weight  and  specific  impulse  known,  the 
total  impulse  and  average  thrust  required  for  boost  can  be  found. 


XT  =  XSP  Wp 

T._ 


With  the  thrust  known  and  an  average  value  of  drag  assumed, 
equation  (3)  can  be  iterated  for  an  improved  estimate  of  the  pro¬ 
pellant  weight  needed. 

The  final  result  is  used  to  calculate  the  booster  combustor 
volume  and  length. 


vb  -  wP/ppnp 


where  np  =  volumetric  packing  factor 


L_  =  4V_  ,  _  2 
B  B/ttD,, 

c 


The  combustor  diameter  is  limited  by  the  missile  maximum 
diameter  and  is  a  design  choice.  The  propellant  density,  p  , 

r 

and  volumetric  packing  factor,  n  ,  are  determined  from  histori- 

r 

cal  data  or  given  for  a  specific  propellant. 
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B.  SUSTAINER  INITIAL  SIZING 

The  sustainer  thrust  required  to  maintain  cruise  is  the 

driving  factor  in  the  sustainer  motor  sizing.  For  initial  sizing 

purposes,  a  level,  constant  velocity  cruise  segment  is  assumed. 

In  this  analysis  the  thrust  required  is  equal  to  the  aerodynamic 

drag  developed  by  the  missile.  The  performance  requirements 

determine  the  maximum  range,  operating  altitude,  and  velocity  of 

the  missile.  From  these  requirements  the  sustainer  burn  time, 

t  ,  can  be  determined, 
s 


where  S  =  distance  traveled  during  boost. 

5 

Since  sustainer  thrust  is  equal  to  drag,  Tg  =  D,  the  total 
impulse  required  becomes, 


The  sustainer  nozzle  is  sized  for  the  operational  altitude 
of  the  missile.  The  thrust  coefficient,  C^,  can  be  expressed  as 


The  constant,  Cd,  is  the  nozzle  efficiency  and  has  been 
determined  from  historical  data  to  be  0.96.  The  nozzle  half 
angle  correction,  is  given  by? 
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NL 
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■ 

r 

_ 

u, 

; 

, 

J 

i 

i 

^  1 

_ 

A  =  (1  +  cos  a) /2 


oC 


where  a  is  the  nozzle  half  angle.  A  nominal  value  used  in  many 
designs  is  15°.  Larger  values  of  a  give  significant  non-axial 
flow  components  in  the  nozzle.  Values  less  than  15°  give  large 
nozzle  lengths  and  therefore  contribute  to  excessive  missile 
weights.  The  combustion  chamber  pressure  is  also  a  design  choice. 
A  lower  usable  value  of  is  200  psia  for  sustainers  and  500  psi 

for  boosters.  The  chamber  pressure  cannot  exceed  the  maximum 
expected  operating  pressure  of  the  missile  (MEOP) . 

The  thrust  of  the  sustainer  can  be  expressed  as  a  function 
of  the  thrust  coefficient,  chamber  pressure,  Pc,  and  throat  area, 
Afc.  The  throat  area  required  to  deliver  the  required  sustainer 
thrust  can  be  determined  from 


T 

s 


'Fd 


The  nozzle  area  ratio  is  a  function  of  the  pressure  ratio 
Pg/P^ ,  where  Pg  is  the  ambient  pressure  if  the  nozzle  expands 
the  flow  fully. 


Y+l 
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From  the  previous  equation  the  exit  area  can  be  determined. 
From  the  throat  and  exit  diameters,  the  nozzle  length  is  found 
by 

D-D 

L  =  _e _ t 

n  2  tan  a 


At  this  point  in  the  design  the  nozzle  exit  area  and  length 
should  be  checked  to  see  if  they  exceed  any  specifications  on 
the  missile;  i.e.,  the  nozzle  exit  area  must  be  less  than  the 
base  area  of  the  missile. 

Finally,  the  delivered  specific  impulse  is  typical  of  the 
propellant  chosen  for  the  sustainer.  The  propellant  weight  is 
then. 


Wp 


(Vs 

ISP 


The  sustainer  combustor  volume  and  length  are  then  determined 
in  the  same  manner  as  the  booster.  The  equations  concerning 
the  sustainer  nozzle  also  apply  to  the  booster  nozzle  with  the 
appropriate  booster  areas,  pressures  and  thrust  substituted. 


C.  ROCKET  MOTOR  CHAMBER  PRESSURE 

As  indicated  earlier  the  delivered  thrust  increases  with 

increasing  chamber  pressure,  P  ;  however,  the  wall  thickness 

c 

in  the  rocket  motor  depends  upon  the  expected  operating  pressure 
within  the  chamber.  If  the  wall  thickness  is  t  inches,  and  the 
yield  stress  of  the  casing  material  is  a ,  it  can  be  shown  that 
the  thickness  required  of  the  motor  casing  is  given  by. 
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t 


P  r 
max 


(5) 


where  r  is  the  chamber  radius,  and  P  is  the  maximum  chamber 

max 

pressure,  which  is  taken  to  be  1.2  P  to  allow  for  variations 

c 

within  the  propellant.  From  equation  (5)  the  weight  of  the  motor 
casing  can  be  determined  if  the  casing  is  cylindrical  and 
t  <  <  r 


W  »  p  V  =  27trJl  tp 
c  c  c  c  c 


Substituting  for  t 


W  = 
c 


2.4-rrr  £  p  P 
c  c  c 


(6) 


The  specific  impulse  of  the  rocket  motor  is  given  by. 

1/2 


V 


I  =  =  —  2g  (-2— )  R  T. 

sp  g„  g„  Y-l  0 


Y-l 

p0  y 
i  -  0 

C 


(7) 


From  equations  (6)  and  (7) ,  it  can  be  seen  that  a  compromise 
must  be  made  in  selecting  the  chamber  pressure;  since  increasing 
P  increases  the  specific  impulse  but  also  increases  the  casing 

C 

weight.  An  attempt  must  be  made  to  find  an  optimum  operating 
pressure. 

If  the  rocket  motor  weight  consists  of  the  nozzle,  propellant 
and  casing, 


W, 


N 


Wp  + 


w„ 


(8) 
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and  the  total  impulse,  which  is  constant  is  given  by, 


I  =  Tt,  =  I  w_ 
b  sp  p 


The  optimum  chamber  pressure  can  be  found  by  minimizing  the 
weight  and  maximizing  the  impulse. 


WM  j  W 

-4.-,  (-J1)  s  _d.  ( — ^ )  =  o 

dp  1  V  dp  W  V  u 

rc  c  sp  p 


differentiating  yields. 


jl  (!m 

^p  dpc  Wp 


WM  d(1/Isp) 
dp„ 


=  0 


(9) 


_1 _ d_ 

I  dp„  ' 


sp  ~*c  "p  W  I 

P  sp 


W„  d  I 
-S—  -5-SE-  0 
2  dP„ 


From  equation  ( 8 ) 


WM  Wm 

-M  =  -£  +  _N  +  i 
Vf  W„  W_ 

P  P  P 


(10) 


Substituting  equation  (10)  into  equation  (9) , 


_i 

«p  dPc 


W  dW 
_ C 

W_2  dpc 


WN  dW 


W. 


M 


d  I 


W_ 


dp. 


JR. 


W  I  dp_ 
p  sp  c 


=  0 
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Substituting  equations  (12)  and  (13)  into  equation  (11) 


Y— 1 

2.4p  W  W  RT.  p  y 

_ -£  .£ - e2_S _  (12.) 

0  CJ  CT  X  ^  D  'd  / 

p  y  sp  *c 


Substituting  for  gc  I  in  the  above  equation yie Ids , 


Y-l 

1 

2Y 

p  (l-y)/y 

(_i)  .1 

p„ 

Pc 

c 

0 


(14) 


If  a  propellant  and  casing  material  are  chosen,  equation  (14) 

can  be  solved  for  the  chamber  pressures,  p  ,  which  minimizes 

c 

WM/I  .  One  interesting  feature  of  equation  (14)  is  that  there 

a 

is  a  minimum  yield  strength  to  density  ration  that  will  yield 

pc 

a  usable  chamber  pressure,  furthermore  this  optimization  yields 
a  very  shallow  minimum.  Therefore,  if  the  thrust  or  exit  dia¬ 
meter  requires  an  increase  in  pc  the  penalty  paid  in  additional 
casing  weight  will  be  small. 

D.  TYPICAL  PROPELLANTS 

From  the  preceding  discussion  it  can  be  seen  that  some  know¬ 
ledge  of  propellants  to  be  used  in  the  rocket  motor  is  needed. 
The  thrust  developed  by  a  rocket  motor  depends  directly  upon  the 
pressure  in  the  chamber. 

T  -  CF  *T  PC 


196 


Since 


w 

P  Xsp  *sp 


dW 


I  dl 


W..  dl. 


c  ,  i  t.t  \  SP  S£  _ M  sp  _ 


-3 -  +  (W  +  W„) 

dp„  c  N 

c 


_  2  dp  I  dp 

I  _  sp  ^c 

sp 


From  equation  <8) 


W  +  W.  -  WM  -  W 
c  N  M  p 


Substituting  for  Wc  +  WN 


dW  W  dl 

c  ps? 


^  '  I.p  d*>c 


=  0 


From  equation  (6) 


dW  2.47rr  l  p_ 
c  _  _ c  c 

dp  “  a 
*c  y 


(11) 


Since 


W  =  irr  I  p 
P  c  p 


dW_  2. 4irp 

c  _  _  c  p 

dp  p  a 

Hp  y 


(12) 


from  equation  (3) 


dl. 


izl 

RTo  (foj  T 
gcpcIsp  pc 


(13) 
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Because  of  this  relationship,  high  pressures  are  needed  for  the 
boost  phase,  and  relatively  low  pressure  are  needed  for  the  sus¬ 
tain  phase  of  flight.  The  pressure  developed  in  a  rocket  motor 
chamber  is  a  function  of  the  burn  rate  and  burn  surface  area  of 
the  propellant.  In  order  to  provide  more  burn  area  for  the  boost 
phase,  the  grain  normally  has  an  internal  star  perforation;  while 
a  sustain  motor  is  normally  an  end  burning  grain  of  solid  cross- 
section.  The  volumetric  loading  of  a  rocket  motor  is  defined  as 
the  ratio  of  the  propellant  volume  to  the  rocket  motor  chamber 
volume . 

_  Grain  volume 
p  Chamber  volume 


Due  to  erosive  burning  effects  the  volumetric  loading  of  a  booster 
is  normally  limited  to  less  than  0.9.  For  efficient  packing 
the  ratio  varies  from  0.7  to  0.9.  The  volumetric  loading  of  an 
end  burning  sustainer  engine  is  normally  1.0.  The  range  of 
propellant  characteristics  is  shown  in  Table  (8-1). 

TABLE  (8-1).  Typical  Propellant  Properties  [14,15]. 


I  (sec) 
sp 

Pp  (lbm/in3) 
Y 


Sustainer 


180-210 
.059-. 062 
1.24  -  1.27 
1.00 


Booster 

210-260 


.062-. 065 
1.22  -  1.26 
0.85 
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E.  DESIGN  EXAMPLE  (BOOSTER) 

From  the  previous  analysis  the  launch  speed  is  M  =  1.5,  and 
the  cruise  speed  is  M  =  3.0  at  an  altitude  of  10,000  feet.  The 
boost  incremental  velocity  is  then, 


AVq  =  1616.10  ft/sec 

If  the  maximum  acceleration  during  boost  is  a  constant  30  gs 
(assume  sea  level  g  ) ,  the  time  of  burn  is, 

C 


s  -  4Va 


max 


1616.10  ft/sec 
30(32.2  ft/sec2) 


tfe  =  1.67  sec 


From  equation  (4)  with  yT  = 


AV. 


B 


'n  Wr 


g  I 


sp 


0 


Assume  I  =  250  sec  from  Table  (I) 
sp 


WL 

l  rr-  =  0.2008 
n  WE 

WE 

Wp  "  WL  (1  -  wf) 


For 


WT  =  794.83 

Li 


W  =  144.57  zbm 
P 
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The  thrust  provided  is 


(T)b 


I  W 

— 2-  =  21,642  Jib. 


Since  the  missile  flies  essentially  at  zero  angle  of 
attack  during  boost,  the  drag  at  the  end  of  boost  is. 


D 


(C„  ) 


q  S 


ref 


D  =  2006  Jlbf 


Therefore  if  an  average  drag  of  1000  lbs  is  assumed,  equation  (3) 
can  be  used  to  find  a  new  WL/WE 


AV„  W 

_ 3 _ EL_ 

g  (T  -  D) 


.2270 


W  =  161.40  Jib 
p  m 

This  gives  a  thrust  of  24,162  lbs.  One  more  iteration  of 
equation  (4)  gives, 

W  -  159-46  Jib 
p  m 

(T)  =  23,872  Jib. 

S  £ 
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Continued  iteration  of  equation  (4)  does  not  change  the  pro¬ 
pellant  weight. 

From  Table  (I);  p  »  .062  lb  /in3 

p  m 

np  =  *85 

The  booster  volume  then  becomes 
VB  =  3025.81  in3 


Allowing  one-half  inch  for  structure,  the  booster  length  becomes. 


L 


B 


4  V 


B 


42.69  inches 


F.  DESIGN  EXAMPLE  (SUSTAINER) 

The  maximum  range  of  the  missile  was  determined  to  be  35  nmiles. 

R  =  212,800  ft 
max 

The  distance  traveled  during  boost  is 

SB  "  *b 
-  V1  +  V2 

V  =  — -  =  2424.15  ft/sec 

S_  -  4048.33  ft 
3 

Therefore,  at  a  cruise  Mach  number  of  3.0,  the  sustainer  time 
of  burn  is 

t  =  64.59  sec. 
s 


201 


The  cruise  drag  consists  of  zero-lift  drag  and  induced  drag. 
The  lift  coefficient  of  the  missile  minus  the  booster  is; 


C 


L 


W 

^  ^ref 


CT  =  0.1270 


From  Chapter  6 


C 


D 


+  C. 


CD  •  .4016 

The  cruise  drag  then  becomes 
D  =  2003  Zbf 


The  total  impulse  required  is  then? 

(I  )  =  Dt  =  129,761  Zb*  -  sec 

T  s  s  f 


From  Table  (I)  the  specific  impulse  is  210  sec;  therefore, 
the  weight  of  propellant  needed  for  cruise  is 


(IJS 

(W  )  «  -=±-5. 

P  S  1  r, 
c  sp 


=  618  Zb. 


m 


The  volume  required  is 

V  =  9967.7  in3 
s 
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The  sustainer  length  becomes 


L  = 
s 


4V 


ttD 


=  140.62  inches 


The  total  rocket  motor  length  is,  L„„,; 

RM 


lrm  -  lb  +  Ls 


L_m  =  183.31  in 
RM 


The  rocket  motor  alone  exeeds  the  total  length  of  the  base¬ 
line  missile.  This  will  lead  to  large  length  to  diameter  ratios 
for  the  entire  missile;  therefore,  the  rocket  motor  length  must 
be  decreased.  That  is  the  subject  of  the  next  section. 


G.  REDUCING  ROCKET  MOTOR  LENGTH 

If  the  rocket  motor  length  of  the  last  section  leads  to 
excessive  length  to  diameter  ratios  for  the  missile,  the  motor 
length  must  be  decreased.  This  may  be  accomplished  in  one  of 
two  ways:  (1)  If  the  mission  profile  has  a  long  cruise  segment, 
the  cruise  altitude  may  be  increased  thereby  decreasing  the  drag 
and  the  total  impulse  required.  (2)  If  the  cruise  altitude 
cannot  be  varied,  the  missile  diameter  must  be  resized  for  pro¬ 
pulsion  considerations.  In  this  design  example  option  (2)  will 
be  selected. 

The  missile  length  can  be  expressed  as 


=  L  +  L_  +  L  +  L„„ 
n  G  c  WH 


+  L 


RM 


(15) 
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L  =  Nose  length 
n 

L_  =  Guidance  section  length 
Lc  =  Control  section  length 
LWH  =  Warhead  section  length 
=  Rocket  motor  length 


The  nose  and  warhead  sections  have  been  previously  defined. 
Since  the  propellant  weight  for  the  range  requirement  is  known, 

the  diameter  necessary  can  be  determined  if  a  maximum  missile 

/ 

length  is  specified.  The  drag  will  increase  slightly  due  to  the 
increased  diameter  but  the  total  impulse,  and  therefore  the  pro¬ 
pellant  weight,  will  change  only  slightly. 

The  rocket  motor  length  can  be  expressed  in  terms  of  the 
diameter  by 


L 


RM 


4W 

— ?? 

V° 


Equation  (15)  then  becomes 


(L) 


L  +  L_  +  L  +  L..„ 
n  G  c  WH 


4Wpc 
_  2 
VD 


The  guidance  and  the  control  sections  are  each  normally  about 

10  percent  of  the  missile  length. 

L  L  4W 

(L)  =  (— )D  +  .1L  +  .1L  +  (-§£)  D  +  -  £  ~ 

D  U  P 


(16) 
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As  mentioned  in  Chapter  2  the  missile  is  normally  designed  for 
an  existing  launcher.  This  launcher  will  have  a  maximum  length 


that  it  can  accomodate,  L  .  If  this  value  is  substituted  into 

IuclX 

equation  (16) ,  the  required  diameter  can  be  found. 

1 .  Design  Example  (Resizing  for  Propulsion) 

From  previous  analysis;  Ln/D  =  2.23 

Lwh/D  =2.50 

The  propellant  weight  from  before  was,  =  777.46  lbs.  If 

the  maximum  length  of  the  missile  is  taken  as  210  inches,  equa¬ 
tion  (16)  becomes, 

4.73D3  -  168D2  +  15,467.08  =  0.0 


The  required  diameter  is  then, 

(D)  -  11.75  inches 

req 

Allowing  0.25  in  for  structure,  the  missile  diameter  becomes 
12  inches.  From  the  equations  for  the  rocket  motor  lengths. 


26.62  in 
90.47  in 
117.09  in 


H.  DESIGN  EXAMPLE  (CHAMBER  PRESSURE) 

From  equation  (14) ,  the  strength  to  density  ratio  can  be 
solved  for  in  terms  of  p  . 


(» 


r  l_Y  i 

)  _  2y 

PQ  Y 

(f)  -1 

)  ~  (y-1) 

min 

p  (^1) 
PC  Pp> 


If  the  minimum  usable  chamber  pressure  is  Pc  =  500  psi,  and 
using  the  propellant  previously  selected. 

Y  =  1.24 

p  =  .062 

P 

pQ  =  10.11 

the  minimum  strength  to  density  required  can  be  found. 


(-J-)  =  225,542.87  in 

C  min 


Inconnel  718  is  a  common  metal  used  in  combustion  chambers 
and  has  the  following  properties  [4] : 


Oy  =  180,000  psi 


pc  = 
a 

JL  = 


.2662  b  /in' 
m 


676,183  >  {-*) 
pc 


mm 


The  q±imum  chamber  pressure  is  then  from  equation  (14) 
pc  =  1132.55  psi 

The  rocket  motor  chamber  wall  thickness  is  given  by  equation  (5) 
as 
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t 


t 


^max 


0.05 


r 


inches 


The  wall  weight  is  given  by  equation  (2) 


W  =  2tt  r£_._tp„ 
c  RM  c 


W  =  59.5  ib 
c  m 


I.  DESIGN  EXAMPLE  (SIZING  THE  NOZZLE) 

For  initial  analysis  it  is  assumed  that  the  flow  is  fully 
expanded,  in  a  nozzle  with  a  half  angle  of  15  degrees.  With 
the  chamber  pressure  of  1132.55  psi  the  thrust  coefficient  at 
altitude  is 

1=11 


C_  =  2.3023 
Fd 

The  throat  area  and  diameter  then  become 

*t*cpc’  9-1552  in2 

Fd 

dt  »  3.4142  in 
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mm 


From  isentropic  tables  for  Pg/Pc  =  .0089 


d 

e 


.  3342 


The  exit  diameter  becomes, 


d  =  10.22  in 
e 

The  nozzle  length  for  a  15  degree  half  angle  is, 
l  =  12.70  in 


The  nozzle  is  shown  in  Figure  (8-2)  .  Figure  (8-3)  is  the  design 
with  resized  diameter  and  length  for  propulsion  and  the  optimum 
wings  determined  in  Chapter  7. 


J.  TAIL  SIZING 

As  mentioned  in  Chapter  7  the  tail  (canard)  sizing  is  accom¬ 
plished  for  zero  static  stability  at  launch.  Now  that  the  missile 
length  has  been  fixed,  the  canard  can  be  sized.  Figure  (8-3) 
is  a  drawing  of  the  design  to  this  point.  The  total  lifting  sur¬ 
face  required  is 

S„  +  S„  -  6.57  ft2 
W  c 


if 


*  =■  sc/sw 


(16) 


From  equation  (5)  in  Chapter  7,  for  CMa  =  0 
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Figure  (8-2) .  Nozzl 


Equations  (16)  and  (17)  above  can  be  solved  for  Sw  and  k. 

Sw  =  4.70  ft2 

k  =  0.40 

Therefore,  S  =1.87 
c 

If  the  same  aspect  ratio  as  found  in  Chapter  7  is  used,  the  wing 
and  canard  planforms  become. 


Wing:  bw  =  1-86  ft 

(c  )  =  5.05  ft 

r  w 

(c)  ,  =  3.60  ft 
w 


Canard:  b  =  1.36  ft 

c 

(c  )  =  2.75  ft 

c 

(c)c  =  2.63  ft 
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IX.  NONLINEAR  AERODYNAMICS  AND  AEROl 


As  mentioned  in  Chpater  5  the  analysis  thus  far  has  been  for 
one  speed  and  small  angles  of  attack.  The  aerodynamics  change 
dramatically  from  subsonic  to  supersonic  flow  and  with  increasing 
angle  of  attack.  This  chapter  presents  the  methods  used  primarily 
by  USAF  Stability  and  Control  DATCOM  -and  covers  all  configura¬ 
tions  and  flight  speeds.  As  will  be  shown,  the  method  becomes 
very  involved  and  therefore  lends  itself  to  programming  on  a 
digital  computer. 

The  lifting  characteristics  of  both  wings  and  bodies  become 
nonlinear  as  the  missile  angle  of  attack  increases  above  10°. 

Up  to  this  point  in  the  thesis  development  only  the  linear 
contribution  has  been  considered.  For  large  angles  of  attack, 
the  nonlinear  term,  which  is  due  to  flow  separation,  must  be 
considered.  The  relative  effect  of  these  terms  on  the  normal 
force  and  pitching  moment  coefficients  is  shown  in  Figure  (9-1) . 


A.  VISCOUS  CROSS-FLOW 


As  can  be  seen  in  Figure  (9-1) ,  at  large  angles  of  attack 
the  forces  acting  are  primarily  nonlinear.  The  nonlinear  term 
is  normally  described  through  the  use  of  the  viscous  cross-flow 
coefficient,  C^ 


Sp  *  planform  area. 
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1 

] 

f 


Figure  (9-1)  .  Relative  contributions  of  linear  and  nonlinear 
terms  [ 19 ] . 

For  an  infinitely  long  circular  cylinder,  the  viscous  cross- 
flow  force  per  unit  length  is 

f  =  2r  C ,  g  (1) 

v  dc  n 

The  viscous  cross-flow  coefficient  is  determined  experimentally 
and  is  a  function  of  normal  Mach  number  and  Reynolds  number. 
Figure  (9-2)  gives  the  cross-flow  drag  coefficient  as  a  function 
of  cross- flow  Mach  number. 


L3 


Figure  (9-2).  Cross-flow  drag  coefficient  [18]. 


M  = 
n 


sin  a 


R  =  R  sin  a 
en  e 

vn  *  VM  sin  “ 


Since  in  real  flow  there  is  spillage  around  the  ends  of  a  finite 
length  cylinder,  the  drag  is  less  than  that  for  an  infinite,  2-D 
cylinder.  The  cross-flow  drag  proportionality  constant,  n,  is 
the  ratio  of  the  drag  coefficient  of  a  finite  cylinder  to  the 
drag  coefficient  of  an  infinite  cylinder. 


n=  (C7T 


c  2-D 


The  proportionality  constant  is  also  determined  experimentally 
and  is  given  in  Figure  (9-3)  .  Notice  that  the  cross-flow  drag 
proportionality  constant  approaches  one  as  the  length  to  dia¬ 
meter  ratio  becomes  large,  or  the  2-D  situation  is  approached. 

The  viscous  cross-flow  force  per  unit  length  for  a  finite 
cylinder  then  becomes 

pV  2 

£v  =  2rncd  — r-  (2> 

c 


Since 


.  2 
sin  a 


=  q 


.  2 
sm  a 


equation  (2)  becomes. 


f  =  2nC.  q  r  sin  a 
v  ac 
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B.  TOTAL  MISSILE  LIFT 

The  remainder  of  this  chapter  is  a  summary  of  references  [17] 
and  [18],  With  the  addition  of  the  force  due  to  viscous  cross- 
flow,  the  total  lift  on  the  body  becomes 


(CN} 


+  C, 


B 


N 


+  C 


B(W) 


N. 


+  C. 


B (T) -a 


N 


B(T)-6 


where , 

(C.T)  =  linear  and  nonlinear  lift  force  on  the  body 

N  B 

C  =  additional  lift  on  the  bodv  due  to  the  wing 

NB(W) 

C„  =  additional  lift  on  the  body  due  to  the  tail 

NB (T) -a 

(due  to  angle  of  attack) 

C.  =  additional  lift  on  the  body  due  to  the  tail 

B(T)-<5 

(due  to  control  deflection) 

The  additional  lift  on  the  body  in  the  presence  of  the  wing  and 
tail  can  be  found  as  described  in  Chapter  6. 

The  wing  lift  is  now  composed  of  a  linear  and  a  nonlinear 
component  and  becomes 


'N. 


=  C 


w 


N. 


+  C. 


W  (B) 


N 


w(v) 
s. 


cn„  ■  *W(B)  sina  +  cd  si”2a 

W  W  ref  n  ref 


The  tail  lift  is  found  in  a  similar  manner  only  with  the 
additional  term  for  any  control  surface  deflection,  6. 

NT  NT(B)-a  NT(B)-S  NT(V) 
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Figure  (9-4).  Drag  coefficient  for  a  flat  plate  normal  to  flow  [18] 


Notice  that  the  linear  terms,  in  the  above  equations  are  multiplied 
by  sina  instead  of  a  as  in  Chapter  6 .  The  linear  theory  from 
which  the  Chapter  6  equations  were  derived  assumes  small  angles 
of  attack  and  therefore  sinasa  .  This  is  not  true  at  higher 
angles  of  attack  and  is  therefore  included  here. 

The  total  normal  force  then  becomes ,  C^ 


'N 


*  C 


N 


+  C 


B 


N, 


W 


+  CN„  +  CI 


The  last  term, 
by 


is  the  lift-loss  due  to  downwash  and  is  given 


(C-VC^T 


(X'fB)31"01  +kW(B)Sin6]  1  (b  -  r)T  SW 
2n  ARt  <f„  -  rw)  Sre£ 


'  i  =  interference  factor 

f  =  votex  location 
w 

the  above  terms  are  found  by  the  methods  of  reference  [13]. 

C.  DRAG  CHARACTERISTICS  OF  A  MISSILE 

The  total  drag  acting  on  a  missile  is  the  sum  of  the  zero- 
lift  drag  and  the  induced  drag  (due  to  angle  of  attack  and/or 
control  surface  deflection) .  The  zero  lift  drag  of  bodies  and 
wings  is  highly  dependent  on  the  missile  speed. 

Three  speed  regimes  are  normally  considered. 
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Subsonic  - -  0  <  <  .8 

Transonic  -  . 8  <  <  1.2 

Supersonic -7—  M  >  1.2 

Because  of  the  empirical  nature  of  the  formulas  for  zero-lift 
drag,  the  procedures  followed  are  essentially  those  of  DATCOM  [16] 
1.  Zero  Lift  Drag 
a.  Subsonic 

At  subsonic  speeds  the  zero-lift  drag  consists  of 
skin  friction  (incompressible)  and  pressure  drag.  The  pressure 
drag  at  subsonic  speeds  is  usually  small  compared  to  the  drag 
due  to  skin  friction.  The  entire  zero-lift  drag  of  a  missile 
at  subsonic  speeds  is  given  by; 


+  CL 

D0  °0B  D0W  dot 


Where  the  components  are  the  body,  wing  and  tail  contributions 
(1)  Body  drag,  .  The  body  zero-lift  drag  is 

r  — 

given  by; 


OB 


=  1.02  C, 


OB 


'B 


1  +  1*5  +  7 
1  +  p72  +  p 


(Swet> 


B 


+  C. 


ref 


B 


where  C  is  found  by  the  methods  of  Chapter  5. 

AB 

The  first  term  is  the  skin  friction  contribution  and  the 
next  two  terms  are  the  pressure  contributions.  f  is  the  body 
fineness  ratio  and  is  given  by; 
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(2)  Wing  drag,  C  .  The  wing  zero-lift  drag  is 


given  by. 


OW 


'ow 


r  t  t  4 

*  C-  1  +  2(J)  +  100  (J) 

fw  L  c  c  J 


4  1  (SWet* 


W 


ref 


(Styfet^  ent;'-re  wetted  area  of  the  wings.  The  thickness 


W 


to  chord  ratio  is  given  by  (t/c) . 


(3)  Tail  drag,  C  .  The  tail  zero-lift  drag  is 


given  by. 


0T 


Cf 

uT  rT 


[l  +  iff)  *  100 (|)4] 


(SWet>, 

Sref 


The  above  analysis  assumes  fully  turbulent  flow  along  all  sur¬ 
faces.  C.  ,  C  _  and  C,.  are  the  local  average  skin  friction 
f B  fW  fT 

coefficients  based  on  the  local  Reynolds  number,  Rg. 

The  reference  lengths  are  the  body  length  and 
wing/tail  mean  aerodynamic  chords.  The  skin  friction  coefficient 
is  then  given  by 


Cf  = 


4.55 


UogiQ  Re) 


T758 


b.  Transonic 

The  total  transonic  zero-lift  drag  is  composed  of 
skin-friction  drag,  transonic  wave  drag,  pressure  drag  and 
base  pressure  drag.  The  compressible  skin  friction  drag  is 
found  by  using  a  correction  factor  on  the  incompressible  co¬ 
efficient,  Q,  ,  found  for  subsonic  flow.  The  compressibility 
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m 


correction  is  found  by  the  methods  of  Chapter  6.  The  skin 
friction  drag  then  becomes 


‘W 


C_  =  1.02  Cf  -5-  - 

DfB  £c  ref 


B 


The  subsonic  pressure  drag  is  the  same  as  before 


1.02  C 

rB 


1.5 

?7I 


+ 


(SWet> 


B 


ref 


The  above  equation  applies  for  Mach  numbers  in  the  range  of 

.8  to  1.0.  The  pressure  drag  then  decreases  linearly  to  zero 

at  a  Mach  number  of  1.2.  The  transonic  wave  drag  C  is  deter- 

rB 

mined  using  Figure  (9-5) ,  which  is  a  function  of  nose  fineness 

ratio  Vd- 

The  total  transonic  zero-lift  drag  of  the  body  then  becomes 

C°0B  =  C°f  +  C°PB  +  CDVB  +  \ 

where  SN  is  the  cross  sectional  area  at  the  nose-body  junction. 

The  total  transonic  zero-lift  drag  of  the  aerodynamic  surfaces 

is  composed  of  the  skin  friction  drag  and  a  drag  increment, 

AC  ,  which  is  the  transonic  wave  drag. 

D0 

Experimental  results  show  little  increase  in  the  viscous 
drag  of  aerodynamic  surfaces  from  the  subsonic  to  the  transonic 
regimes;  therefore  the  skin  friction  and  pressure  drag  is  the 
same  as  for  subsonic  flow  and  is  given  by 
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CD  =  Cf  1  +  2{|>  +  100(t)4]  T*1 

D0W  £w  1  c  c  J  Sref 


The  wave  drag  of  the  aerodynamic  surfaces  if  found  from 
Figure  (9-6)  and  is  a  function  of  AR  and  M.  For  swept  wings 
the  Mach  number  used  in  Figure  (9-6)  is  given  by. 


M'  =  M  [cos  A  c/4] 1//2 


The  transonic  wave  drag  increment  is  then  given  by 


2.5  SM 

AC  =  AC  '  [cos  A  c/4]  ^ - 

U0W  OW  aref 


where  AC  is  obtained  from  Figure  (9-6) .  The  tail  zero-lift 
U0W 

drag  is  found  in  the  same  manner  as  the  wing.  The  total  zero-lift 
drag  (transonic)  is  then  given  by; 


Cn  =  C_  +  C„  +  AC_  +  C„  +  AC^ 

Do  dob  dow  dow  dot  dot 


c.  Supersonic 

The  supersonic  zero-lift  drag  of  a  missile  is  deter¬ 
mined  empirically  by  assuming  a  parabolic  variation  of  Cn  with 

0 

Mach  number  between  1.2  and  3.0.  The  resulting  equation  is  for 
the  entire  missile  and  is  given  by 


C_  =  C  '  -  1.7209  (C  "  -  C  ’)  +  1.5708  (C  "  -  C  ’)  ^M 
D0  0  D0  D0  0  D0 
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Figure  (9-6) .  Transonic  C  for  unswept 


where 


c°; = v 


M  =  1. 2 


S  *  V 


M  =  3.0 


CD’  is  determined  by  using  the  methods  for  transonic  flow. 

The  magnitude  of  the  supersonic  wave  drag  is  highly 

dependent  upon  the  nose  shape.  For  this  reason  two  methods 

are  used  to  determine  C  " 

D0 

Method  I;  This  method  is  for  blunted  ogives,  pointed  ogives  and 
blunted  cones.  In  this  method 


For  Vd  4  -5 

For  ln/d  8'° 


C  "  =  C  ' 
\  \ 

\  ■  S’ 


+  AC, 


M  *  .8 


'OW 


+  ACD  +  CA 
0T  AB 


The  above  values  are  determined  from  the  methods  of  transonic 
flow. 


For  L 


N/D 


.5 


to  8.0 


K1  in  the  above  equation  is  derived  empirically  and  is  given  in 
Figure  (9-7) . 

Method  II:  This  method  is  for  pointed  conical  noses.  In  this 
method  C_"  is  determined  by 
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c  " 
Do 


=  (C„  ) 


+  AC, 


+  AC, 


M  =  .8 


OW 


OT 


+  (CD  } 
UVN 


+  C, 


M  =  3 


B 


The  first  three  terms  in  the  above  equation  are  found  by  transonic 

flow  methods.  The  forebody  wave  drag,  C  ,  for  M  =  3.0  is  found 

°VN 

using  Figure  (9-8).  The  nose  semi-vertex  angle,  8  ,  is  the  same 
as  ct  in  Chapter  6.  In  all  flow  regimes  the  base  pressure  drag 
is  found  as  in  Chapter  6 . 

2 .  Induced  Drag 

The  induced  drag  due  to  angle  of  attack  depends  upon  the 
flight  regime  the  missile  is  in.  For  M  <  .85  and  AR  >  3.0  the 
induced  drag  is , 

(CL)2 

C  =  — - - 

Di  irARe 

where  e  is  the  Oswald  efficiency  factor  a  0.7. 

The  tactical  missile  normally  has  an  aspect  ratio  of  less  than 
3.0.  For  all  flight  speeds  with  AR  <  3.0  the  induced  drag  is 

C_  =  C_  tan  a 
Di  L 

The  component  induced  drag  coefficients  are  found  in  the  same 
manner  as  above  using  component  lift  coefficients. 

3.  Total  Drag 

From  the  zero-lift  drag  and  induced  drag  the  missile 
total  drag  becomes, 


C 


D 


229 


Figure  (9-8)  .  Wave  drag  of  pointed  conical  nose 


D.  PITCHING  MOMENT  CHARACTERISTICS 


The  pitching  moment  of  a  missile  is  equal  to  the  sum  of 
the  pitching  moments  due  to  lift  and  drag  forces  acting  on  the 
body,  wings  and  tails.  If  only  small  angles  of  attack  are  of 
interest,  the  pitching  moment  is  due  primarily  to  the  lift 
forces.  The  methods  presented  here  are  for  all  angles  of  attack. 

1.  Body 

The  body  alone  pitching  moment  about  its  center  of  gravity 
is  given  from  slender  body  theory  and  viscous  cross- flow  theory 
as 


VB  -  SB(L  -  *CG> 
Sref  Lref 


sm  a  cos 


!  + 


n  C. 


S  x  _ 
P  \  /  cG 


ref 


-)  (- 


sm  a 


Jref 


S  =  Planform  area 
P 

Xp  =  Centroid  of  planform  area 
V_  =  Body  volume 
Sg  =  Base  area 


With  C 


:  given  the  center  of  pressure  of  the  body  becomes 

B 


"M. 


B 


“N, 


Jref 


2 .  Wing  (Fixed  Surface) 

The  center  of  pressure  location  for  the  wing  must  be 
found  before  the  pitching  moment  can  be  specified.  The 
center  of  pressure  of  the  various  lift  components  are  found 
by  the  method  of  Pitts,  Nielsen  and  l&atiari  [13].  The  lift 
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r»  ■  -  '  ^ 

components  are  the  same  as  in  the  lift  section.  The  center 

of  pressure  of  the  lift  on  the  wing  in  the  presence  of  the 

body,  (x)  ,  is  found  from  Figure  (9-10)  if  the  flow  is 

W(B) 

subsonic  and  Figure  (9-11)  if  the  flow  is  supersonic. 

The  center  of  pressure  of  the  additional  lift  on  the  body 

in  the  presence  of  the  wing,  (x)_  ,  is  found  using  Figure 

B  (W) 

(9-12)  if  the  flow  is  subsonic.  If  the  flow  is  supersonic 
the  center  of  pressure  is  found  from  Figure  (9-13)  or  Figure 
(9-14)  depending  on  the  parameters; 


BAR(1  +  X)  (1  +  j^|) 

X  =  taper  ratio 

m  =  co  +  A ■ 

Lib 

If  the  center  of  pressure  reference  is  moved  to  the  nose,  the 
following  expressions  result. 
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EFFECTIVE  ASPECT  RATIO,  /3AR 


Figure  (9-11} .  Wing  center  of  pressure 

(supersonic)  [18] . 
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EFFECTIVE  ASPECT  RATIO,  0AR 


Figure  (9-13) .  Body  center  of  pressure 


From  the  above  figure  the  entire  pitching  moment  about 


the  nose  of  the  missile  due  to  the  forces  acting  on  the  wing  is. 


;  «  [  (c 

“w  ^(B)  ^WCV)  i 


^°4<B) 


+  C_  +  C_  )cosct  +  AC_  sina]  (  _ 

L"'”'  L-‘  w-a  Lref 


(x  ) 

♦  tCT  cos.]  (  °P  Blw)) 
lB(W)  ref 


The  viscous  lift,  C  and  lift-loss  due  to  downwash, 

LW(V) 

CT  ,  are  not  shown  in  the  figure  but  are  the  same  forces  as 

lj . 

1 

found  in  the  lift  section.  These  forces  are  assumed  to  be 

acting  through  the  center  of  pressure,  (x  ) 

p  W(B) 

An  average  center  of  pressure  of  the  wing  due  to  all  forces 
acting  on  it  can  now  be  found  as 


C  L  , 

.  .  _  _ ”w  re£ _ _ 

'xcp„  Tel  +  C_  )  cosa  +  AC  sina 
w  LW  Li  DW-a 


where 

Lw  lw(b)  lb(w)  lw(v) 

3.  Tail  (Control  Surface) 

The  tail  pitching  moment  about  the  nose  of  the  missile 
is  found  in  the  same  way  except  now  a  control  surface  deflec¬ 
tion  must  be  included.  The  equations  now  become 
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C  =  C ( C_  +  Cl  +  C  ) cosa  +  C 

t  LT(B)  T  Li  LT ( B ) - 5 

(xcpJ 

(AC  +  AC  )  sina]  ( — = - •T— -  ■  ) 

DT-a  DT—  <S  ref 


(x  ) 
cp 


+  [C_  cosa  +  CT  ]  ( — = - 2ilL 

^B  (T)  ^B(T)-6  ref 


CMt  Lref 

\  =  _ _ ± _ _ 

cp  (CL  +  CT  )cosa  +  CT +  (AC„  +  AC„  Tsina 


T-a  i 


LT-6  °T-a  DT-6 


where. 


C  =  C  +  C  +  C 

LT-a  LT(B)  LB(T)  LT(V) 


Cr  -  C  +  C 

LT-6  LT(B)-<5  LB(T)-S 


The  wing  and  tail  pitching  moments  above  were  found  about 
the  nose  of  the  missile.  Transferring  the  pitching  moments 
about  the  center  of  gravity  the  complete  missile  pitching 
moment  becomes 


C  SB  C  +  C 

M  mb  *w 


’  x  _  —  (x  ) 

x  „  -  (x  ) 

<=G  cp  w 

4.  n 

cG  cp'T 

(x  ) 

cp  w  J 

«T 

(x  ) 

L  Cp  T  J 

E.  AER01  DESCRIPTION 

It  can  be  seen  from  the  preceding  description  of  the  aero¬ 
dynamics  of  a  missile,  that  the  process  of  obtaining  a  complete 


240 


description  of  the  aerodynamic  coefficients  becomes  an  involved 
task.  This  was  the  justification  for  initially  using  linear 
aerodynamics.  To  complete  an  analysis  a  fast  method  of  pre¬ 
dicting  the  static  aerodynamic  characteristics  is  needed. 

The  method  needs  to  include  both  linear  and  nonlinear  contri¬ 
butions  as  well  as  interference  factors  and  must  be  applicable 
to  all  speed  regimes.  This  analysis  lends  itself  to  program¬ 
ming  on  a  digital  computer.  AEROl  is  a  modification  of  the 
program  AEROCF  which  was  developed  at  the  Naval  Air  Development 
Center  by  Mr.  F.A.  Kuster,  Jr.  [17].  The  program  is  essentially 
as  written  except  for  calculation  of  planform  areas  and  centroids. 
The  program  was  also  modified  for  use  on  the  Naval  Postgraduate 
School  CP/CMS  system. 

AEROl  consists  of  a  main  program  and  three  subroutines.  The 
inputs  to  the  program  are  the  geometric  characteristics  of  the 
missile,  the  flight  conditions,  engine  and  inlet  type  and  the 
protuberance  drag.  The  output  consists  of  the  aerodynamic  co¬ 
efficients  for  lift,  drag  and  moment.  The  component  contribution 
to  these  coefficients  are  also  given.  The  component  and  overall 
center  of  pressure  are  also  determined. 

Subroutine  GEOSUB;  This  subroutine  calculates  the  missile 
wetted  area  and  the  Reynolds  number  per  foot  based  on  the  flight 
altitude. 

Subroutine  CLASUB;  This  subroutine  calculates  the  aerodynamic 
surface  lift-curve  slopes. 
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Subroutine  CATSUB;  This  subroutine  calculates  center  of 
pressure  locations,  cross-flow  drag  coefficients,  and  interference 
factors . 

Main  Program;  The  main  program  assumes  the  control  surface 
is  the  tail.  This  is  regardless  of  the  method  of  control  (Wing, 
Tail,  Canard).  Because  of  this,  care  must  be  taken  to  input  the 
right  data  for  the  tail.  For  instance,  if  the  missile  is  a  wing 
control  missile,  the  wing  data  is  input  as  the  tail  and  the 
tail  data  as  the  wing.  Figure  (9-15)  and  Table  (9-1)  give  a 
complete  listing  of  the  input  data.  Table  (9-II)  is  a  list  of 
the  output  data.  Appendix  A  is  a  listing  of  the  program  as  modi¬ 
fied  for  use  on  the  Naval  Postgraduate  School  IBM  360  computer. 

1.  Verification  of  AEROl 

Before  using  the  program  an  attempt  was  made  to  verify 
its  prediction  techniques  and  find  any  pitfalls  in  its  use.  To 
accomplish  this  the  program  was  run  for  various  configurations 
for  which  experimental  data  were  available,  and  the  results  were 
compared.  The  comparisons  are  shown  in  Figures  (9-16)  to  (9-25) 
from  references  [19]  -  [22]  which  are  NASA  technical  notes  and 
memoranda  which  report  results  of  wind  tunnel  tests  for  various 
body-wing- tail  combinations.  In  Figures  (9-16)  to  (9-25)  the 
solid  lines  are  AEROl  predictions, 
a.  NASA  TN  D-6996 

This  technical  note  presents  a  method  of  predicting 
aerodynamic  characteristics  for  bodies  alone  at  angles  of  attack 
from  0  to  180  degrees.  Several  nose-body canbinations  are  given. 
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TABLE  (9-1) 
Input  Data 


Variable 

Name 

Format 

1 

ICSC 

12 

2 

INOSE 

12 

3 

IDT 

12 

4 

IM 

12 

5 

IAL 

12 

6 

NBODY 

12 

7 

ISWPW 

12 

8 

IAFBW 

12 

9 

IWEPW 

12 

10 

NWING 

12 

11 

ISWPT 

12 

12 

IAFBT 

12 

13 

ISWEPT 

12 

14 

NTAIL 

12 

15 

XLAMW 

F10. 5 

16 

CLAMW 

F10. 5 

17 

BW 

F10. 5 

18 

CROOTW 

F10. 5 

19 

SW 

F10.5 

20 

XMACW 

F10. 5 

21 

XWING 

F10. 5 
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Meaning 

Type  of  control  (wing, 
tail,  canard) 

Nose  shape  (ellipsoidal, 
cone,  ogive) 

Number  of  control  surface 
deflections 

Numer  of  Mach  numbers 

Number  of  angles  of  attack 

Number  of  configurations 

Wing  shape  (delta, nondelta) 

Missile  body  after  wing 

Leading  edge  sweep  indicator 

Number  of  wings 

Tail  shape  (delta, nondelta) 

Missile  body  after  tail 

Leading  edge  sweep  indicator 

Number  of  tails 

Wing  taper  ratio 

Wing  leading  edge  sweep 

Wing  span 

Wing  root  chard 

Wing  exposed  area 

Wing  mean  geometric  chord 

Distance  to  wing  leading 
edge 


"1 


Variable 

Name 

Format 

Meaning 

22 

TOVCW 

F10.5 

Wing  thickness  to  chord  ratio 

23 

XLAMT 

F10. 5 

Tail  taper  ratio 

24 

CLAMT 

F10.5 

Tail  leading  edge  sweep 

25 

BT 

F10. 5 

Tail  span 

26 

CROOTT 

F10.5 

Tail  root  chord 

27 

ST 

F10. 5 

Tail  exposed  area 

28 

XMACT 

F10. 5 

Tail  mean  geometric  chord 

29 

XTAIL 

F10.5 

Distance  to  tail  leading 
edge 

30 

TOVCT 

F10.5 

Tail  thickness  to  chord  ratio 

31 

HT 

F10. 5 

Altitude 

32 

D 

F10. 5 

body  diameter 

33 

XL 

F10.5 

Body  length 

34 

XLNOSE 

F10.5 

Nose  length 

35 

XCG 

F10. 5 

Center  of  gravity  location 

36 

AREA 

F10.5 

Reference  area 

37 

XREF 

F10.5 

Reference  length 

38 

ENGINE 

F10. 5 

Engine  code 

39 

INLET 

F10. 5 

Inlet  code 

40 

BETA 

F10 .5 

Boattail  angle 

41 

DBASE 

F10.5 

Base  diameter 

42 

DJET 

F10.5 

Nozzle  exit  diameter 

43 

XLABOD 

F10. 5 

Boattail  length 

44 

CDPROT 

F10. 5 

Proturberance  drag 
coefficient 

L 
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TABLE  (9-II) 
Output  Data 


Variable  Name 
AL 

CLTOT 

CDTOT 

CLWP 

CLBW 

CLTP 

CLBT 

CLB 

CDI 

CNWP 

CNTP 

CLTD 

CDTD 

CN 

CA 

XCPW 

XCPT 

XCP 

CM 


Meaning 

Angle  of  attack 

Total  coefficient  of  lift 

Total  coefficient  of  drag 

Wing  panel  coefficient  of  lift 

Additional  lift  on  body  due  to 
wing  coefficient 

Tail  panel  coefficient  of  lift 

Additional  lift  on  body  due  to 
tail  coefficient 

Body  alone  lift  coefficient 

Induced  drag  coefficient 

Wing  panel  normal  force  coefficient 

Tail  panel  normal  force  coefficient 

Coefficient  of  lift  due  to  tail 
deflection 

Coefficient  of  drag  due  to  tail 
deflection 

Total  normal  force  coefficient 

Total  axial  force  coefficient 

Wing  center  of  pressure 

Tail  center  of  pressure 

Total  missile  center  of  pressure 

Total  pitching  moment  about  center 
of  gravity 
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Figures  ( 9 — 16 )  to  (9-18)  compare  the  coefficients  predicted  by 
AER01  with  those  obtained  for  body  number  9  in  the  NASA  report. 

The  normal  force  coefficient  is  predicted  well  throughout  the 
range.  The  moment  and  axial  force  trends  are  predicted  by  the 
program  but  large  errors  exist  throughout  the  range  of  angles 
of  attack. 

b.  NASA  TM  X-2367 

This  technical  memorandum  investigates  the  aero¬ 
dynamic  characteristics  of  various  cruciform  body-wing  combina¬ 
tions.  The  coefficients  for  this  configuration  agree  very  well 
with  experimental  values  up  to  10  degrees  angle  of  attack.  Beyond 
this  value  the  lift  and  moment  coefficients  predicted  by  AEROl 
exceed  the  experimental  values  by  as  much  as  10  percent.  Al¬ 
though  the  exact  cause  of  this  error  was  not  investigated,  it  may 
be  partially  explained  by  the  nose  shape  of  the  body.  The  nose 
is  a  combination  ogive  and  cone.  For  purposes  of  AEROl  it  was 
assumed  an  ogive.  The  greater  presented  area  of  the  ogive  would 
contribute  to  a  higher  and  Cm  through  the  cross-flow  terms  in 
these  coefficients.  Figures  (9-19)  through  (9-21)  compare  the 
predicted  with  the  experimentally  determined  coefficients  for 
the  wing-tail  configuration  of  this  reference  at  M  =  .9. 

c.  NASA  TM  X-2780  and  NASA  TM  X-2289 

These  technical  memoranda  investigate  the  aerodynamics 
of  a  delta  wing  missile  using  tail  control  and  a  tail-less  cruci¬ 
form  missile.  As  shown  in  Figures  (9-22)  to  (9-25)  there  is 
excellent  agreement  between  the  experimental  values  of  the  aero¬ 
dynamic  coefficients  and  those  predicted  by  AEROl. 
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Figure  (9- 


Figure  (9-20)  .  Tail-body  drag  force 


Figure  {9-21) .  Tail-body  moment 


Figure  (9-22) .  Wing-body  normal  force. 


Figure  (9-23) .  Wing-body  moment 


Figure  (9-24) .  Complete  missile  lift  force. 


Figure  (9-25) .  Complete  missile  drag 


F.  COMPONENT  WEIGHTS 


As  with  the  gross  weight,  the  component  weights  can  be 
determined  through  the  use  of  parametric  regression  equations. 
If  the  dimensions  of  the  components  are  known,  the  following 
formulas  can  be  used  to  determine  the  component  weights  [7] . 
Aero  Surface  Weight 

The  weight  of  one  wing  panel  is  given  by, 

WAS  =  6-77483  (Eas)1,02  (ARAS}‘56  (1) 

2 

Eas  Exposed  are  of  one  panel,  ft 
ARas  Aspect  ratio  of  one  panel 

Body  Structure  Weight 

«BS  *  -0604  (LBS>'64  (DBS>1'77  (2) 

L__  Length  of  body  to  be  covered  by 

structure.  This  does  not  include 
the  rocket  motor  unless  a  separate 
structure  surrounds  the  motor 

casing.  (inches) 

_  Diameter  of  body  structure 
BS 

(inches) 

Internal  Systems  Weight 

Wig  =  .00485  (WQ)-74  (Ljg)1*00  (DIS)*42  (3) 

W_  Gross  weight  of  the  missile 


L^g  Length  of  subsystem  (inches) 
DIS  Diameter  of  missile  (inches) 


In  the  above  equation,  the  internal  system  includes  guidance, 
avionics  and  control. 

1.  Design  Example  (Component  Sizing) 

If  the  guidance  and  the  control  sections  are  kept  at 
10  percent  of  the  total  length,  and  the  remaining  components 
are  as  sized  previously,  the  design  is  as  shown  in  Figure  (9- 
Since  these  components  are  considered  internal  systems,  their 
weights  can  be  determined  from  equation  (3) . 


W  =  .00485  (1809) '79  ( . 1L)  (12.5) *42 

cont 


W  t  =  W  .  .  =  83.64  lb 
cont  guid 

The  wing  and  canard  weights  are  found  using  equation  (1)  and 
Figure  (9-26) . 


Figure  (9-26) .  Aero  surface  weight. 


27)  . 
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Figure  (9-27) .  Final  design. 


The  wing  weight  is  four  times  the  weight  of  one  panel 


W„  =  43.42  lb 
W 

Similarly  the  canard  weight  is 

W  =  9.88  lb 
c 

If  the  entire  body  is  covered  by  structure,  the  body  structure 
weight  becomes 

W  =  .0604  (232)’64  (12. 5)1'77 
st 

W  =  172.38  lb 
st 


The  following  component  table  can  now  be  constructed. 


Component 

Length  (in) 

Weight  (lbs) 

Center  of  Gravity 

Control 

23.2 

83.64 

39.48 

Guidance 

23.2 

83.64 

62.68 

Warhead 

31.25 

28.26 

89.91 

Sustainer 

87.51 

618.92 

149.29 

Booster 

26.31 

159.46 

206.20 

Motor  Casing 

113.82 

59.50 

162.44 

Wings 

43.42 

211.40 

Canards 

9.88 

44.08 

Structure 

232.00 

172.38 

116.00 

(in) 
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The  complete  missile  weight  is  then  1264.62  lbs. 

The  center  of  gravity  from  the  above  table  is 

xcG  =  139.44  in 

The  center  of  gravity  at  the  end  of  boost  is 
(xcG)  =  129.85  in 

G.  DESIGN  EXAMPLE  (FINAL  ANALYSIS) 

From  the  design  parameters  so  far,  a  complete  description 
of  the  missile  can  be  determined.  Figure  (9-27)  is  a  drawing 
of  the  missile.  The  launch  conditions  for  the  missile  were 
specified  as. 


h  =  10,000  ft 
W_  =  1264.62 

Cj 

x  =  139.44  in 
cG 

From  these  conditions  the  input  data  for  AEROl  is  shown  in 
Table  (9-III) .  The  format  is  the  same  as  that  of  the  output  of 
AEROl  and  is  printed  as  a  check  to  ensure  that  the  input  data 
was  entered  properly.  The  output  is  shown  in  Table  (9-IV) . 

Figure  (9-28)  is  a  plot  of  the  coefficient  of  moment  versus  angle 
of  attack  for  the  launch  condition  and  for  the  beginning  of 
cruise.  The  center  of  gravity  has  moved  forward  approximately 
10  inches  for  the  beginning  of  cruise  so  that  the  stability  has 
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TABLE  (9-III) .  AEROl  Input 


r\  .D  n  N  n.  q  JN(T 
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o 
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Figure  (9-28).  Static  stability. 


increased  as  indicated  in  this  figure.  The  performance  objective 
from  the  beginning  of  the  design  was  to  provide  a  31  g  maneuver 
capability.  Table  (9-  V)  shows  the  output  of  AER01  for  Mach 
numbers  of  1.5  to  3.0  and  control  deflections  of  0.0  to  30  degrees. 
From  this  output  the  trimmed  normal  force  can  be  found  and  the 
corresponding  load  factor  is  then. 


'N 


q  s 


n  = 


TR 


ref 


w 


The  following  values  of  maneuver  load  factor  were  found  using 
AERO 1 . 


M 

q  (lb/ft2) 

n(g’s) 

6  req 

2.5 

18.80 

6369 

85.00 

10  deg. 

2.0 

24.50 

4076 

70.99 

10  deg. 

1.5 

12.00 

2293 

19.56 

10  deg. 

As  indicated  in  the  above  table,  the  missile  meets  the  maneuvering 
specifications  of  the  operational  requirements. 
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X.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  has  presented  the  methods  and  general  procedures 
for  the  conceptual  design  of  tactical  missiles.  As  mentioned, 
this  is  not  necessarily  the  best  design.  Now  that  the  design 
procedure  has  been  attempted  once,  the  reason  for  this  is  readily 
seen.  The  process  is  one  of  continuous  compromise.  As  was  shown 
in  Chapter  7,  the  optimum  wing  for  lift  is  not  the  best  wing 
for  minimizing  drag.  We  also  saw  in  Chapter  8  that  increasing 
chamber  pressure  increases  thrust  at  the  expense  of  increased 
weight.  This  thesis  has  tried  to  point  out  some  of  these  areas 
of  compromise  and  present  methods  to  deal  with  them. 

Areas  which  were  not  covered  which  need  to  be  investigated 
in  conceptual  design  are  structures ,  radar  cross-section  and  cost. 
With  the  increased  emphasis  on  survivability  and  the  decreasing 
budget,  it  becomes  increasingly  important  to  define  the  effects 
of  these  areas  on  design  early  in  the  process. 

The  complexity  of  the  design  process  and  the  need  to  obtain 
timely  and  accurate  information  have  made  it  ideally  suited  for 
the  digital  computer.  The  AEROCF  program  used  in  this  thesis 
is  part  of  a  large  scale  computer  program  (MISSYN)  which  con¬ 
sists  of  modules  for  each  section  of  the  design  analysis. 

With  a  good  understanding  of  the  theory  and  methods  used  in 
missile  design,  the  computer  aided  design  program  with  graphics 
capability  gives  the  designer  the  capacity  to  make  intelligent 


design  interations  and  almost  instantaneously  see  the  effect 
of  the  change  on  all  areas  of  design.  A  limited  example  of 
this  can  be  seen  in  the  use  of  AER01 .  A  change  in  the  per¬ 
formance  requirements  for  maneuver  capability  would  require 
a  redesign  of  the  lifting  surfaces.  A  change  in  the  lifting 
surface  design  would  change  the  drag  characteristics  of  the 
missile  and  therefore, the  propulsion  requirements.  The  AER01 
program  coupled  with  a  similar  propulsion  module  would  allow 
the  designer  to  make  the  changes  and  instantly  see  the  penalty 
or  savings  in  propellant  weight. 

One  pass  at  the  design  has  been  accomplished  in  this  thesis. 

As  was  seen  throughout  the  process,  decisions  in  one  area 

affect  the  design  in  others.  For  this  reason  the  design 

process  becomes  an  iterative  one.  The  final  design  of  the 

first  iteration  is  .the  baseline  missile  for  the  second  itera- 

« 

tion  and  the  design  is  started  again.  By  making  several 
passes  through  the  loop,  the  solution  converges  on  the  final 
design. 
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